bzoj3238 [Ahoi2013]差异
题目链接
树状数组+单调栈
计算出每个后缀的前面、后面第一个h[]比它小的(前闭后开),乘起来计算答案
1 #include<algorithm> 2 #include<iostream> 3 #include<cstdlib> 4 #include<cstring> 5 #include<cstdio> 6 #include<string> 7 #include<cmath> 8 #include<ctime> 9 #include<queue> 10 #include<stack> 11 #include<map> 12 #include<set> 13 #define rre(i,r,l) for(int i=(r);i>=(l);i--) 14 #define re(i,l,r) for(int i=(l);i<=(r);i++) 15 #define Clear(a,b) memset(a,b,sizeof(a)) 16 #define inout(x) printf("%d",(x)) 17 #define douin(x) scanf("%lf",&x) 18 #define strin(x) scanf("%s",(x)) 19 #define LLin(x) scanf("%lld",&x) 20 #define op operator 21 #define CSC main 22 typedef unsigned long long ULL; 23 typedef const int cint; 24 typedef long long LL; 25 using namespace std; 26 void inin(int &ret) 27 { 28 ret=0;int f=0;char ch=getchar(); 29 while(ch<'0'||ch>'9'){if(ch=='-')f=1;ch=getchar();} 30 while(ch>='0'&&ch<='9')ret*=10,ret+=ch-'0',ch=getchar(); 31 ret=f?-ret:ret; 32 } 33 char s[500050]; 34 int n,sa[500050],c[500050],t[500050],t2[500050],h[500050]; 35 void build_sa(int m) 36 { 37 int *x=t,*y=t2; 38 re(i,0,n-1)x[i]=s[i],c[x[i]]++; 39 re(i,1,m)c[i]+=c[i-1]; 40 rre(i,n-1,0)sa[--c[x[i]]]=i; 41 for(int k=1;k<=n;k<<=1) 42 { 43 int p=0; 44 rre(i,n-1,n-k)y[p++]=i; 45 re(i,0,n-1)if(sa[i]>=k)y[p++]=sa[i]-k; 46 re(i,0,m-1)c[i]=0; 47 re(i,0,n-1)c[x[y[i]]]++; 48 re(i,1,m-1)c[i]+=c[i-1]; 49 rre(i,n-1,0)sa[--c[x[y[i]]]]=y[i]; 50 swap(x,y); 51 x[sa[0]]=0,p=1; 52 re(i,0,n-1)x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?p-1:p++; 53 if(p>=n)return ; 54 m=p; 55 } 56 } 57 int rank[500050]; 58 void build_height() 59 { 60 int k=0; 61 re(i,0,n-1)rank[sa[i]]=i; 62 re(i,0,n-1) 63 { 64 if(k)k--; 65 if(!rank[i])continue; 66 int j=sa[rank[i]-1]; 67 while(s[i+k]==s[j+k])k++; 68 h[rank[i]]=k; 69 } 70 } 71 int l[500050],r[500050],sta[500050]; 72 LL ans; 73 int main() 74 { 75 strin(s); 76 n=strlen(s); 77 build_sa(256); 78 build_height(); 79 int top=0; 80 re(i,0,n-1) 81 { 82 while(top&&h[sta[top]]>=h[i])top--; 83 l[i]=(top?sta[top]+1:0); 84 sta[++top]=i; 85 } 86 top=0; 87 rre(i,n-1,0) 88 { 89 while(top&&h[sta[top]]>h[i])top--; 90 r[i]=(top?sta[top]-1:n-1); 91 sta[++top]=i; 92 } 93 ans+=(LL)(1+n)*n/2*(n-1); 94 re(i,0,n-1) 95 ans-=2LL*(i-l[i]+1)*(r[i]-i+1)*h[i]; 96 printf("%lld",ans); 97 return 0; 98 }
转载于:https://www.cnblogs.com/HugeGun/p/5225566.html
bzoj3238 [Ahoi2013]差异相关推荐
- BZOJ3238[Ahoi2013]差异
BZOJ3238[Ahoi2013]差异 题目描述 n<=500000n<=500000n,都是小写字母 输入 一行,一个字符串S 输出 一行,一个整数,表示所求值 Solution 公式 ...
- BZOJ3238: [Ahoi2013]差异
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 4840 Solved: 2298 [Submit][Stat ...
- [BZOJ3238] [AHOI2013] 差异 - 后缀自动机
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 1968 Solved: 896 [Submit][Statu ...
- [BZOJ3238][AHOI2013]差异 [后缀数组+单调栈]
题目地址 - GO-> 题目大意: 给定一个长度为 nn 的字符串SS,令TiTi表示它从第ii个字符开始的后缀,求以下这个式子的值: ∑1≤i<j≤nlen(Ti)+len(Tj)−2× ...
- bzoj3238 [Ahoi2013]差异 后缀自动机
这个题虽然看起来比较麻烦,但其实有一点搞懂了就好做了 每一个后缀的前缀的lcp一定是这个串里的子串, 对于全部子串,如果它不是后缀,那一定可以找到这个串对应的后缀, 而这样的子串和后缀子串是没有区别的 ...
- [BZOJ3238][Ahoi2013]差异
Time Limit: 20 Sec Memory Limit: 512 MB Description Input 一行,一个字符串S Output 一行,一个整数,表示所求值 Sample Inpu ...
- bzoj3238 [Ahoi2013]差异 后缀自动机
题意:给出一个串,求其中任意两个字串的lcp的总和. 我们可以对于这个串建一颗后缀自动机,实际上,他的parent边树就是一颗后缀树,我们在后缀树上统计答案,设f表示right集合的大小,可以理解为后 ...
- bzoj3238: [Ahoi2013]差异(后缀自动机)
传送门 后缀自动机好题. 题意: 做法:samsamsam 废话 考虑翻转字串,这样后缀的最长公共前缀等于前缀的最长公共后缀. 然后想到parentparentparent树上面两个串的最长公共后缀跟 ...
- 【bzoj3238】[Ahoi2013]差异
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 3425 Solved: 1559 [Submit][Stat ...
最新文章
- linux el5 el6 区别,Oracle ACFS Support on Linux詳解
- JAVA窗帘_HomeControl 智能家具系统,包括灯光,窗帘的控制,设备,房间,情景模式的添加 Java Develop 240万源代码下载- www.pudn.com...
- mac 思科 链路聚合_EtherChannel Cisco 端口聚合详解
- 云原生生态周报 Vol. 13 | Forrester 发布企业级容器平台报告
- 洛谷 P1849 [USACO12MAR]拖拉机Tractor
- 推荐系统遇上深度学习(五)--DeepCross Network模型理论和实践
- 1.Java 面试题整理(基础篇一)
- CUBLAS变量解释(1)
- python flask web部署_webapp开发之使用apache部署python flask webapp
- 【十二】Jmeter:“CSV 数据文件设置”参数化请求出现中文乱码
- Security+ 学习笔记43 无线网络
- 转: 在CentOS 6.X 上面安装 Python 2.7.X
- 历年蓝桥杯日期类问题总结
- android 清除activity栈
- 三菱plc编程软件gx-developer在win10中无法使用及工程无法建立解决方法
- win10桌面计算机打不开,win10“此电脑”双击打不开,任何文件夹都打不开
- 【美学集】色彩之冷暖色
- SQL Server查询IP地址等信息
- LeetCode 二分查找
- VLAN的划分以及三层交换机理论,OMG,太详细了吧,看它!
热门文章
- Taro+react开发(26)taro基本骨架
- 前端学习(3050):vue+element今日头条管理-表格组件基本使用
- [vue] vue中是如何使用event对象的?
- [css] 如何取消页面中选中的文字?
- 工作177:表单重置项目处理
- 前端学习(2260)高校案例
- “约见”面试官系列之常见面试题之第六十五篇之ajax如何解析jsonp(建议收藏)
- 前端学习(95):ps基本操作与图片格式
- Java面试题18 牛客 假定Base b = new Derived();
- 第八十七期:爬了知乎“沙雕问题”,笑死个人!