POJ 3693 后缀数组+RMQ
来源:互联网 发布:数据库系统概论电子书 编辑:程序博客网 时间:2024/06/11 16:23
思路:
论文题
后缀数组&RMQ
有一些题解写得很繁
//By SiriusRen#include <cmath>#include <cstdio>#include <cstring>#include <algorithm>using namespace std;#define N 100050int 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;char s[N];void SA(){ memset(cntA,0,sizeof(cntA)); for(int i=1;i<=n;i++)cntA[s[i]]++; for(int i=1;i<=256;i++)cntA[i]+=cntA[i-1]; for(int i=1;i<=n;i++)sa[cntA[s[i]]--]=i; rk[sa[1]]=1; for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]); for(int l=1;rk[sa[n]]<n;l<<=1){ memset(cntA,0,sizeof(cntA)); memset(cntB,0,sizeof(cntB)); for(int i=1;i<=n;i++) cntA[A[i]=rk[i]]++, cntB[B[i]=i+l<=n?rk[i+l]:0]++; for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1]; for(int i=n;i;i--)tsa[cntB[B[i]]--]=i; for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i]; rk[sa[1]]=1; 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]]); } for(int i=1,j=0;i<=n;i++){ j=j?j-1:0; while(s[i+j]==s[sa[rk[i]-1]+j])j++; ht[rk[i]]=j; }}void RMQ(){ for(int i=1;i<=n;i++)f[i][0]=ht[i]; for(int i=1;i<=16;i++) for(int j=1;j<=n;j++)if(j+(1<<i)-1<=n) f[j][i]=min(f[j][i-1],f[j+(1<<(i-1))][i-1]);}int LCP(int x,int y){ int tempx=rk[x],tempy=rk[y]; if(tempx>tempy)swap(tempx,tempy); tempx++;int t=Log[tempy-tempx+1]; return min(f[tempx][t],f[tempy-(1<<t)+1][t]);}int main(){ Log[0]=-1; for(int i=1;i<N;i++)Log[i]=i&(i-1)?Log[i-1]:Log[i-1]+1; while(scanf("%s",s+1)&&s[1]!='#'){ maxx=0,n=strlen(s+1),SA(),RMQ(); for(int i=1;i<=n/2;i++){ for(int j=1;i+j<=n;j+=i){ int r=LCP(j,j+i),tot=r/i+1,st=j-(i-r%i); if(st>0&&r%i&&LCP(st,st+i)>=r)tot++; if(tot>maxx)maxx=tot,stk[top=1]=i; else if(tot==maxx)stk[++top]=i; } } int len=-1,st; for(int i=1;i<=n&&len==-1;i++) for(int j=1;j<=top;j++) if(LCP(sa[i],sa[i]+stk[j])>=(maxx-1)*stk[j]) {len=stk[j],st=sa[i];break;} printf("Case %d: ",++cases); for(int i=0;i<len*maxx;i++)printf("%c",s[i+st]);puts(""); }}
0 0
- poj 3693 后缀数组+RMQ
- POJ 3693 后缀数组+RMQ
- POJ 3693 后缀数组+RMQ
- poj-3693(后缀数组+RMQ)
- POJ 3693 Maximum repetition substring(RMQ+后缀数组)
- poj 3693 Maximum repetition substring 后缀数组+RMQ
- POJ 3693 Maximum repetition substring 后缀数组 暴力 rmq
- POJ 3693 Maximum repetition substring 后缀数组 + RMQ预处理
- POJ 3693 Maximum repetition substring后缀数组加上RMQ
- poj 3693 Maximum repetition substring (后缀数组+RMQ)
- POJ 3693 Maximum repetition substring RMQ+后缀数组
- POJ 3693 Maximum Repetition Substring <后缀数组 + RMQ>
- 【后缀数组,RMQ】poj3693
- 后缀数组+RMQ
- POJ3693后缀数组+RMQ
- poj3693 后缀数组+RMQ
- HDU6194(后缀数组 + rmq)
- POJ 3693 Maximum repetition substring(08合肥 RMQ+后缀数组)
- could not get lock /var/lib/dpkg/lock -open等
- a label can only be part of a statement and a declaration is not a statement
- hadoop的小疑问:Map执行未结束便开始执行Reduce操作?
- leetcode--36. Valid Sudoku
- 计算单词出现的频率C语言实现的
- POJ 3693 后缀数组+RMQ
- javascript初探四
- android5.0沉浸状态栏
- List of devices getting official LineageOS Roms [ Full ]
- 读取assets目录下的文件
- Latex:TexStudio的使用
- 虚树
- 1
- linux 文本处理命令