【BZOJ】 3238: [Ahoi2013]差异
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3238
求:$${\sum _{i=1}^{n-1}\sum _{j=i+1}^{n}len(T_i)len(T_j)-2lcp(T_i,T_j)}$$
$T_i$表示从第$i$个字符开始的后缀。
原式${=n(n+1)(n-1)-\sum _{i=1}^{n-1}\sum _{j=i+1}^{n}2lcp(T_i,T_j)}$
这不就是NOI2015品酒大会的第一问啦,利用SAM构出后缀树,每个后缀对应一个后缀树上的点,${lcp(x,y)=deep[lca(x,y)]}$,树形统计一下即可。
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #include<vector> 5 #include<cstdlib> 6 #include<cmath> 7 #include<cstring> 8 using namespace std; 9 #define maxn 500010*2 10 #define llg long long 11 #define SIZE 26 12 #define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout); 13 llg n,m,ans,ans1[maxn],size[maxn]; 14 char s[maxn]; 15 16 struct SAM 17 { 18 struct 19 { 20 llg len,f,ch[SIZE]; 21 void init() 22 { 23 len=0,f=-1; 24 memset(ch,0xff,sizeof(ch)); 25 } 26 }e[maxn<<1]; 27 28 vector<llg>a[maxn<<1]; 29 llg idx,last,size[maxn<<1]; 30 31 void init(){idx=last=0; e[idx++].init(); memset(size,0,sizeof(size));} 32 33 llg newnode() {e[idx].init(); return idx++;} 34 35 void insert(llg c) 36 { 37 llg end=newnode(),tmp=last; 38 e[end].len=e[last].len+1; size[end]=1; 39 for (;tmp!=-1 && e[tmp].ch[c]==-1;tmp=e[tmp].f) e[tmp].ch[c]=end; 40 if (tmp==-1) e[end].f=0; 41 else 42 { 43 llg nxt=e[tmp].ch[c]; 44 if (e[tmp].len+1==e[nxt].len) e[end].f=nxt; 45 else 46 { 47 llg np=newnode(); 48 e[np]=e[nxt]; 49 e[np].len=e[tmp].len+1; 50 e[nxt].f=e[end].f=np; 51 for (;tmp!=-1 && e[tmp].ch[c]==nxt;tmp=e[tmp].f) e[tmp].ch[c]=np; 52 } 53 } 54 last=end; 55 } 56 57 void link_fa() {for (llg i=1;i<idx;i++) a[e[i].f].push_back(i);} 58 59 void dp(llg x) 60 { 61 llg w=a[x].size(),v; 62 for (llg i=0;i<w;i++) dp(a[x][i]); 63 for (llg i=0;i<w;i++) 64 { 65 v=a[x][i]; 66 ans1[e[x].len]+=size[x]*size[v]; 67 size[x]+=size[v]; 68 } 69 } 70 71 }sam; 72 73 int main() 74 { 75 yyj("bzoj3238"); 76 scanf("%s",s); 77 n=strlen(s); 78 sam.init(); 79 for (llg i=n-1;i>=0;i--) sam.insert(s[i]-'a'); 80 sam.link_fa(); 81 sam.dp(0); 82 ans=n*(n+1)*(n-1); 83 ans/=2; 84 for (llg i=1;i<=n;i++) ans-=i*ans1[i]*2; 85 cout<<ans<<endl; 86 return 0; 87 }
转载于:https://www.cnblogs.com/Dragon-Light/p/6385455.html
【BZOJ】 3238: [Ahoi2013]差异相关推荐
- BZOJ 3238: [Ahoi2013]差异 [后缀数组 单调栈]
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 2326 Solved: 1054 [Submit][Stat ...
- bzoj 3238: [Ahoi2013]差异(后缀数组+单调栈)
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 3443 Solved: 1562 [Submit][Stat ...
- @bzoj - 3238@ [Ahoi2013]差异
目录 @description@ @solution@ @accepted code@ @details@ @description@ 给定一个长度为 n 的字符串 S,令 Ti 表示它从第 i 个字 ...
- bzoj 3238: [Ahoi2013]差异
Description Input 一行,一个字符串S Output 一行,一个整数,表示所求值 Sample Input cacao Sample Output 54 HINT 2<=N< ...
- BZOJ 3238 [Ahoi2013]差异
后缀自动机+后缀树 有一种不难想到的后缀数组的做法,详见黄学长题解 然而我做这题的目的主要是熟悉一下后缀自动机和后缀树- - 一个有趣的性质:反串的后缀自动机的parent树就是正串的后缀树,其树上的 ...
- 3238: [Ahoi2013]差异
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 2106 Solved: 953 [Submit][Statu ...
- 【BZOJ】3238: [Ahoi2013]差异
[题意]给定长度为n的小写字母字符串,令Ti表示以i开头的后缀,求Σ[Ti+Tj-2*lcp(Ti,Tj)],1<=i<j<=n. [算法]后缀自动机 [题解]Σ(Ti+Tj)只与n ...
- 【bzoj 3238】差异
传送门~ 解题思路 求后缀数组以后搞个单调栈. 因为height数组的意义让这个单调栈很烦,注意要理清思路再写. 代码: #include<algorithm> #include<i ...
- BZOJ3238: [Ahoi2013]差异
3238: [Ahoi2013]差异 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 4840 Solved: 2298 [Submit][Stat ...
最新文章
- E8.Net 2005工作流平台版本发布
- qs.parse()、qs.stringify()使用方法
- 简述html 布局的原理,css布局原理与实现-2019年9月4日20时
- java图形界面猜字游戏,java程序,猜字游戏,希望大神帮忙
- react学习(67)--git 屏蔽文件不被追踪
- 网络编程C#篇(二):Socket无连接简单实例
- ApacheCN 未来发展方向(暂定)2019.8.2
- 蚂蚁集团回应 “ IPO 推迟”;iPhone12系列支持北斗系统;LLVM 11.0.0 发布 | 极客头条
- 中国AI人才仅美国1/17,网易要填这个坑?
- [转]文本编辑软件UltraEdit v16.20官方简/繁体中文版下载+注册码和破解方法
- matlab 蜗杆轮廓,基于MATLAB的直廓环面蜗杆的可视化
- 【IoT】产品设计:贴片电子电容封装类型:0202封装、0402封装、0603封装、0805封装、1206封装、1210封装、1812封装、2010封装、2512封装
- 机器学习 第一节 第一课
- CSDN 上传资源已经存在
- java实现qq邮箱发送附件和图片
- 【渝粤题库】陕西师范大学500010 量子力学 作业(专升本)
- 90后,第一批接棒白酒的“年轻人”
- 聊聊chatbot那些事
- 自由人NFT平台,全球首个以NFT+农业+元宇宙项目
- 百度、高德地图功能进一步探索之—电子围栏绘制(三)