POJ 3693 Maximum repetitionsubstring(后缀数组:循环子串)
来源:互联网 发布:php读取txt文件前字符 编辑:程序博客网 时间:2024/06/17 00:34
POJ 3693 Maximum repetitionsubstring(后缀数组:循环子串)
http://poj.org/problem?id=3693
题意:给你一个长N的字符串,现在要你找出该字符串中具有最大循环节数目的连续子串,如果有多个答案输出字典序最小的.
分析: 罗穗骞《后缀数组——处理字符串的有力工具》论文例题,不过论文中讲解的比较难懂.
本题是SPOJ687的加强版,可以先看SPOJ687的解答:
http://blog.csdn.net/u013480600/article/details/23974119
假设我们现在已经记录了所有可以使得循环节最多的循环节长度(从小到大记录).
假设存在那个子串,该子串字典序最小且循环节最多.且该字串肯定是某一个后缀的前缀,我们只需要按sa数组去枚举所有后缀(因为sa[1]是字典序最小的后缀),然后对于每个sa[i],从小到大枚举可行循环节的长度(因为对于相同字典序时,长度短的字典序小),然后用LCP查询看看sa[i]与sa[i]+len的LCP是否>=len*ans.如果是,马上退出即可.(就算后面还有>=len*ans的情况,它们不是sa字典序大,就是len长度长).
AC代码:
#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int maxn=100000+10;struct SuffixArray{ char s[maxn]; int sa[maxn],rank[maxn],height[maxn]; int t1[maxn],t2[maxn],c[maxn],n; int mm[maxn];//mm[i]=k表示2^k是正好小于i的那个k. 而2^(k+1)>i了. int dmin[maxn][20]; void build_sa(int m) { int i,*x=t1,*y=t2; for(i=0;i<m;i++) c[i]=0; for(i=0;i<n;i++) c[x[i]=s[i]]++; for(i=1;i<m;i++) c[i]+=c[i-1]; for(i=n-1;i>=0;i--) sa[--c[x[i]]]=i; for(int k=1;k<=n;k<<=1) { int p=0; for(i=n-k;i<n;i++) y[p++]=i; for(i=0;i<n;i++)if(sa[i]>=k) y[p++]=sa[i]-k; for(i=0;i<m;i++) c[i]=0; for(i=0;i<n;i++) c[x[y[i]]]++; for(i=1;i<m;i++) c[i]+=c[i-1]; for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]] = y[i]; swap(x,y); p=1,x[sa[0]]=0; for(i=1;i<n;i++) x[sa[i]]= y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]? p-1:p++; if(p>=n) break; m=p; } } void build_height() { int i,j,k=0; for(i=0;i<n;i++)rank[sa[i]]=i; for(i=0;i<n;i++) { if(k)k--; j=sa[rank[i]-1]; while(s[i+k]==s[j+k])k++; height[rank[i]]=k; } } void initMin() { mm[0]=-1; for(int i=1;i<=n;i++) { dmin[i][0]=height[i]; mm[i]=((i&(i-1))==0)?mm[i-1]+1:mm[i-1]; } for(int j=1;j<=mm[n];j++) for(int i=1;i+(1<<j)-1<=n;i++) dmin[i][j]=min(dmin[i][j-1] , dmin[i+(1<<(j-1))][j-1]); } int RMQ(int L,int R)//取得范围最小值 { int k=mm[R-L+1]; return min(dmin[L][k] , dmin[R-(1<<k)+1][k]); } int LCP(int i,int j)//求后缀i和j的LCP最长公共前缀 { int L=rank[i],R=rank[j]; if(L>R) swap(L,R); L++; return RMQ(L,R); }}sa;int cnt;int res[maxn];int main(){ int kase=1; while(scanf("%s",sa.s)==1) { int n=strlen(sa.s); if(n==1&&strcmp(sa.s,"#")==0) break; sa.s[n]=0; sa.n=n+1; sa.build_sa(128); sa.build_height(); sa.initMin(); int ans=1; cnt=0; for(int len=1;len<n;len++) { for(int j=0;j+len<n;j+=len) { int k=sa.LCP(j,j+len); int now=k/len+1; int jj=j-(len-k%len);//往前推到jj位置 if(jj>=0) if(sa.LCP(jj,jj+len)/len+1 > now) now++; if(now>ans) { ans=now; cnt=0; res[cnt++]=len; } else if(now == ans) res[cnt++]=len; } } int st,best_len=-1; for(int i=1;i<n&&best_len==-1;i++) { for(int j=0;j<cnt;j++) { int len = res[j]; if(sa.sa[i]+len<n && sa.LCP(sa.sa[i],sa.sa[i]+len)>=(ans-1)*len ) { best_len=len; st=sa.sa[i]; break; } } } printf("Case %d: ",kase++); for(int i=0;i<best_len*ans;i++)printf("%c",sa.s[st+i]); printf("\n"); } return 0;}
0 0
- POJ 3693 Maximum repetitionsubstring(后缀数组:循环子串)
- POJ - 3693 Maximum repetition substring(后缀数组求重复次数最多的连续重复子串)
- POJ 3693 Maximum repetition substring(后缀数组求最长重复子串)
- POJ 3693 Maximum repetition substring (后缀数组+RMQ 求重复最多的连续子串)
- POJ 3693 Maximum repetition substring(后缀数组[重复次数最多的连续重复子串])
- POJ 3693 Maximum repetition substring (后缀数组+RMQ 求重复最多的连续子串)
- POJ 3693 Maximum repetition substring 后缀数组求重复次数最多子串
- poj 3693 Maximum repetition substring //后缀数组
- poj 3693 Maximum repetition substring (后缀数组)
- poj 3693 Maximum repetition substring(后缀数组)
- poj 3693 Maximum repetition substring(后缀数组)
- 【后缀数组】【poj 3693】Maximum repetition substring
- POJ 3693 Maximum Repetition Substring 后缀数组
- poj 3693/hdu 2459 Maximum repetition substring spoj 687. Repeats ( 后缀数组 重复次数最多的连续重复子串)
- Maximum repetition substring+POj+后缀数组之求重复次数最多的连续重复子串
- poj 2406 Power Strings 【KMP求最小循环节】【后缀数组求连续重复子串】
- Poj 3693 & Hdu 2459 Maximum repetition substring (08合肥Online 后缀数组+RMQ 重复次数最多的连续重复子串)
- POJ 题目 3693 Maximum repetition substring(后缀数组+RMQ+枚举求最小字典序的重复次数最多的子串)
- 以自己的亲身经历献给那些在沈阳找工作看见沈阳深讯招聘信息的程序员们
- Cracking The Coding Interview5.1
- 第7周作业1-循环大战
- c++ 调用 python
- MYSQL用户权限管理(Grant,Revoke)
- POJ 3693 Maximum repetitionsubstring(后缀数组:循环子串)
- win8系统打不来chm文档文件的解决方案
- ifndef define endif
- 初级单例模式
- 第7周作业2-极值计算
- extern "C"
- 移植u-boot-2009-11->S3c44b0公板
- LeetCode(Word Search)
- Mybatis 示例之 复杂(complex)属性(property)