题目链接

树状数组+单调栈

计算出每个后缀的前面、后面第一个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]差异相关推荐

  1. BZOJ3238[Ahoi2013]差异

    BZOJ3238[Ahoi2013]差异 题目描述 n<=500000n<=500000n,都是小写字母 输入 一行,一个字符串S 输出 一行,一个整数,表示所求值 Solution 公式 ...

  2. BZOJ3238: [Ahoi2013]差异

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 4840  Solved: 2298 [Submit][Stat ...

  3. [BZOJ3238] [AHOI2013] 差异 - 后缀自动机

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 1968  Solved: 896 [Submit][Statu ...

  4. [BZOJ3238][AHOI2013]差异 [后缀数组+单调栈]

    题目地址 - GO-> 题目大意: 给定一个长度为 nn 的字符串SS,令TiTi表示它从第ii个字符开始的后缀,求以下这个式子的值: ∑1≤i<j≤nlen(Ti)+len(Tj)−2× ...

  5. bzoj3238 [Ahoi2013]差异 后缀自动机

    这个题虽然看起来比较麻烦,但其实有一点搞懂了就好做了 每一个后缀的前缀的lcp一定是这个串里的子串, 对于全部子串,如果它不是后缀,那一定可以找到这个串对应的后缀, 而这样的子串和后缀子串是没有区别的 ...

  6. [BZOJ3238][Ahoi2013]差异

    Time Limit: 20 Sec Memory Limit: 512 MB Description Input 一行,一个字符串S Output 一行,一个整数,表示所求值 Sample Inpu ...

  7. bzoj3238 [Ahoi2013]差异 后缀自动机

    题意:给出一个串,求其中任意两个字串的lcp的总和. 我们可以对于这个串建一颗后缀自动机,实际上,他的parent边树就是一颗后缀树,我们在后缀树上统计答案,设f表示right集合的大小,可以理解为后 ...

  8. bzoj3238: [Ahoi2013]差异(后缀自动机)

    传送门 后缀自动机好题. 题意: 做法:samsamsam 废话 考虑翻转字串,这样后缀的最长公共前缀等于前缀的最长公共后缀. 然后想到parentparentparent树上面两个串的最长公共后缀跟 ...

  9. 【bzoj3238】[Ahoi2013]差异

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 3425  Solved: 1559 [Submit][Stat ...

最新文章

  1. linux el5 el6 区别,Oracle ACFS Support on Linux詳解
  2. JAVA窗帘_HomeControl 智能家具系统,包括灯光,窗帘的控制,设备,房间,情景模式的添加 Java Develop 240万源代码下载- www.pudn.com...
  3. mac 思科 链路聚合_EtherChannel Cisco 端口聚合详解
  4. 云原生生态周报 Vol. 13 | Forrester 发布企业级容器平台报告
  5. 洛谷 P1849 [USACO12MAR]拖拉机Tractor
  6. 推荐系统遇上深度学习(五)--DeepCross Network模型理论和实践
  7. 1.Java 面试题整理(基础篇一)
  8. CUBLAS变量解释(1)
  9. python flask web部署_webapp开发之使用apache部署python flask webapp
  10. 【十二】Jmeter:“CSV 数据文件设置”参数化请求出现中文乱码
  11. Security+ 学习笔记43 无线网络
  12. 转: 在CentOS 6.X 上面安装 Python 2.7.X
  13. 历年蓝桥杯日期类问题总结
  14. android 清除activity栈
  15. 三菱plc编程软件gx-developer在win10中无法使用及工程无法建立解决方法
  16. win10桌面计算机打不开,win10“此电脑”双击打不开,任何文件夹都打不开
  17. 【美学集】色彩之冷暖色
  18. SQL Server查询IP地址等信息
  19. LeetCode 二分查找
  20. VLAN的划分以及三层交换机理论,OMG,太详细了吧,看它!

热门文章

  1. Taro+react开发(26)taro基本骨架
  2. 前端学习(3050):vue+element今日头条管理-表格组件基本使用
  3. [vue] vue中是如何使用event对象的?
  4. [css] 如何取消页面中选中的文字?
  5. 工作177:表单重置项目处理
  6. 前端学习(2260)高校案例
  7. “约见”面试官系列之常见面试题之第六十五篇之ajax如何解析jsonp(建议收藏)
  8. 前端学习(95):ps基本操作与图片格式
  9. Java面试题18 牛客 假定Base b = new Derived();
  10. 第八十七期:爬了知乎“沙雕问题”,笑死个人!