思路:

论文题

后缀数组&RMQ

有一些题解写得很繁

  1. //By SiriusRen
  2. #include <cmath>
  3. #include <cstdio>
  4. #include <cstring>
  5. #include <algorithm>
  6. using namespace std;
  7. #define N 100050
  8. int n,cases,cntA[N],cntB[N],A[N],B[N],rk[N],sa[N],tsa[N],ht[N],f[N][18],Log[N],stk[N],top,maxx;
  9. char s[N];
  10. void SA(){
  11. memset(cntA,0,sizeof(cntA));
  12. for(int i=1;i<=n;i++)cntA[s[i]]++;
  13. for(int i=1;i<=256;i++)cntA[i]+=cntA[i-1];
  14. for(int i=1;i<=n;i++)sa[cntA[s[i]]--]=i;
  15. rk[sa[1]]=1;
  16. for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);
  17. for(int l=1;rk[sa[n]]<n;l<<=1){
  18. memset(cntA,0,sizeof(cntA));
  19. memset(cntB,0,sizeof(cntB));
  20. for(int i=1;i<=n;i++)
  21. cntA[A[i]=rk[i]]++,
  22. cntB[B[i]=i+l<=n?rk[i+l]:0]++;
  23. for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1];
  24. for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;
  25. for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];
  26. rk[sa[1]]=1;
  27. for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]);
  28. }
  29. for(int i=1,j=0;i<=n;i++){
  30. j=j?j-1:0;
  31. while(s[i+j]==s[sa[rk[i]-1]+j])j++;
  32. ht[rk[i]]=j;
  33. }
  34. }
  35. void RMQ(){
  36. for(int i=1;i<=n;i++)f[i][0]=ht[i];
  37. for(int i=1;i<=16;i++)
  38. for(int j=1;j<=n;j++)if(j+(1<<i)-1<=n)
  39. f[j][i]=min(f[j][i-1],f[j+(1<<(i-1))][i-1]);
  40. }
  41. int LCP(int x,int y){
  42. int tempx=rk[x],tempy=rk[y];
  43. if(tempx>tempy)swap(tempx,tempy);
  44. tempx++;int t=Log[tempy-tempx+1];
  45. return min(f[tempx][t],f[tempy-(1<<t)+1][t]);
  46. }
  47. int main(){
  48. Log[0]=-1;
  49. for(int i=1;i<N;i++)Log[i]=i&(i-1)?Log[i-1]:Log[i-1]+1;
  50. while(scanf("%s",s+1)&&s[1]!='#'){
  51. maxx=0,n=strlen(s+1),SA(),RMQ();
  52. for(int i=1;i<=n/2;i++){
  53. for(int j=1;i+j<=n;j+=i){
  54. int r=LCP(j,j+i),tot=r/i+1,st=j-(i-r%i);
  55. if(st>0&&r%i&&LCP(st,st+i)>=r)tot++;
  56. if(tot>maxx)maxx=tot,stk[top=1]=i;
  57. else if(tot==maxx)stk[++top]=i;
  58. }
  59. }
  60. int len=-1,st;
  61. for(int i=1;i<=n&&len==-1;i++)
  62. for(int j=1;j<=top;j++)
  63. if(LCP(sa[i],sa[i]+stk[j])>=(maxx-1)*stk[j])
  64. {len=stk[j],st=sa[i];break;}
  65. printf("Case %d: ",++cases);
  66. for(int i=0;i<len*maxx;i++)printf("%c",s[i+st]);puts("");
  67. }
  68. }

POJ 3693 后缀数组+RMQ的更多相关文章

  1. poj 3693 后缀数组 重复次数最多的连续重复子串

    Maximum repetition substring Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 8669   Acc ...

  2. POJ 3693 后缀数组

    题目链接:http://poj.org/problem?id=3693 题意:首先定义了一个字符串的重复度.即一个字符串由一个子串重复k次构成.那么最大的k即是该字符串的重复度.现在给定一个长度为n的 ...

  3. [poj 3693]后缀数组+出现次数最多的重复子串

    题目链接:http://poj.org/problem?id=3693 枚举长度L,看长度为L的子串最多能重复出现几次,首先,能出现1次是肯定的,然后看是否能出现两次及以上.由抽屉原理,这个子串出现次 ...

  4. POJ - 2406 ~SPOJ - REPEATS~POJ - 3693 后缀数组求解重复字串问题

    POJ - 2406 题意: 给出一个字符串,要把它写成(x)n的形式,问n的最大值. 这题是求整个串的重复次数,不是重复最多次数的字串 这题很容易想到用KMP求最小循环节就没了,但是后缀数组也能写 ...

  5. POJ 3693 (后缀数组) Maximum repetition substring

    找重复次数最多的字串,如果有多解,要求字典序最小. 我也是跟着罗穗骞菊苣的论文才刷这道题的. 首先还是枚举一个循环节的长度L,如果它出现两次的话,一定会包含s[0], s[L], s[2L]这些相邻两 ...

  6. 【uva10829-求形如UVU的串的个数】后缀数组+rmq or 直接for水过

    题意:UVU形式的串的个数,V的长度规定,U要一样,位置不同即为不同字串 https://uva.onlinejudge.org/index.php?option=com_onlinejudge&am ...

  7. spoj687 REPEATS - Repeats (后缀数组+rmq)

    A string s is called an (k,l)-repeat if s is obtained by concatenating k>=1 times some seed strin ...

  8. POJ - 2774~POJ - 3415 后缀数组求解公共字串问题

    POJ - 2774: 题意: 求解A,B串的最长公共字串 (摘自罗穗骞的国家集训队论文): 算法分析: 字符串的任何一个子串都是这个字符串的某个后缀的前缀. 求 A 和 B 的最长 公共子串等价于求 ...

  9. HDU2459 后缀数组+RMQ

    题目大意: 在原串中找到一个拥有连续相同子串最多的那个子串 比如dababababc中的abababab有4个连续的ab,是最多的 如果有同样多的输出字典序最小的那个 这里用后缀数组解决问题: 枚举连 ...

随机推荐

  1. Spring学习笔记(一) 简介

    版权声明 本文是摘自IBM上Naveen Balani的一篇文章,原文请点击此处:http://www.ibm.com/developerworks/cn/java/wa-spring1/ Sprin ...

  2. POJ 1466 最大独立点集

    思路:匈牙利 n-ans/2; // by SiriusRen #include <cstdio> #include <cstring> #define N 505 using ...

  3. Rsync 服务器搭建

    Rsync简介 rsync 是一个 Unix 系统下的文件同步和传输工具. 它具备以下特性: 1. 能更新整个目录和树和文件系统 2. 有选择性的保持符号链链.硬链接.文件属于.权限.设备以及时间 等 ...

  4. 格式化日期字符串 FormatSettings使用

    如果 你想要得到 YYYY-MM/DD 这样的字符串 你肯定说这太简单了  直接 ShowMessage(FormatDateTime('YYYY-MM/DD',now)); 运行结果 YYYY-MM ...

  5. BZOJ2440: [中山市选2011]完全平方数(莫比乌斯+容斥原理)

    2440: [中山市选2011]完全平方数 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 4920  Solved: 2389[Submit][Sta ...

  6. RecyclerView的刷新和加载更多

    1.RecyclerView :出现也不知道多久了,没怎么使用过,上次写的笔记乱七八糟的,再次仔细的整理下.   使用需加入依赖  compile 'com.android.support:recyc ...

  7. Sqlite简单操作

    1.student类 public class Student { int id; String name; String sex; String address; int money; public ...

  8. Unity 移动键Q的三种用法 For Mac,Windows类同

    拖动整个场景:三指 (任何模式下)ALT+三指:旋转当前镜头 (任何模式下)双指前后滑动:缩放镜头 ps1:Q键移动的游戏场景,W移动的是游戏对象 ps2:三指 = 左键拖动

  9. 【洛谷1131】 [ZJOI2007]时态同步

    树形结构的入门题,一遍DFS即可解决,注意答案开long long #include<cstdio> #include<iostream> #include<algori ...

  10. 斯特林公式--取N阶乘近似值

    斯特林公式(Stirling's approximation)是一条用来取n的阶乘的近似值的数学公式.一般来说,当n很大的时候,n阶乘的计算量十分大,所以斯特林公式十分好用,而且,即使在n很小的时候, ...