传送门
后缀数组好题。
考虑枚举循环节长度lenlenlen。
然后考虑枚举循环节的起点来更新答案。
但是直接枚举每次O(n)O(n)O(n)。
考虑枚举len∗k+1len*k+1len∗k+1作为起点。
然后用len∗(k−1)+1len*(k-1)+1len∗(k−1)+1和len∗k+1len*k+1len∗k+1之间的某一段有可能使循环节长度再加一。
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#define ri register int
using namespace std;
const int N=1e5+5;
char s[N];
int sa2[N],cnt[N],n,Log[N],m,sa[N],rk[N],ht[N],st[N][17],pos[N],top=0;
inline void Sort(){
for(ri i=1;i<=m;++i)cnt[i]=0;
for(ri i=1;i<=n;++i)++cnt[rk[i]];
for(ri i=2;i<=m;++i)cnt[i]+=cnt[i-1];
for(ri i=n;i;--i)sa[cnt[rk[sa2[i]]]--]=sa2[i];
}
inline void getsa(){
for(ri i=1;i<=n;++i)rk[i]=s[i]-'a'+1,sa2[i]=i;
m=130,Sort();
for(ri w=1,p=0;m!=n;p=0,w<<=1){
for(ri i=n-w+1;i<=n;++i)sa2[++p]=i;
for(ri i=1;i<=n;++i)if(sa[i]>w)sa2[++p]=sa[i]-w;
Sort(),swap(rk,sa2),rk[sa[1]]=p=1;
for(ri i=2;i<=n;++i)rk[sa[i]]=(sa2[sa[i]]==sa2[sa[i-1]]&&sa2[sa[i]+w]==sa2[sa[i-1]+w])?p:++p;
m=p;
}
for(ri i=1,j,k=0;i<=n;ht[rk[i++]]=k)for(k?--k:k,j=sa[rk[i]-1];s[i+k]==s[j+k];++k);
for(ri i=1;i<=n;++i)st[i][0]=ht[i];
for(ri j=1;j<=16;++j)for(ri i=1;i+(1<<j)<=n;++i)st[i][j]=min(st[i][j-1],st[i+(1<<(j-1))][j-1]);
}
inline int rmq(int l,int r){return min(st[l][Log[r-l+1]],st[r+1-(1<<Log[r-l+1])][Log[r-l+1]]);}
inline int Query(int l,int r){return rmq(min(rk[l],rk[r])+1,max(rk[l],rk[r]));}
inline void query(){
int tim=1,ans=1,len=0;
for(ri T=1;T<n;++T){
for(ri l=1,r=T+1;r<=n;l+=T,r+=T){
int L=Query(l,r),stp=L/T+1,pl=l-(T-L%T),pr=pl+T;
if(pl>0&&L%T)if(Query(pl,pr)>=L)++stp;
if(stp>tim)tim=stp,pos[top=1]=T;
else if(stp==tim)pos[++top]=T;
}
}
for(ri i=1;!len&&i<=n;++i)for(ri j=1,T=pos[j];j<=top;T=pos[++j])if(Query(sa[i],sa[i]+T)>=(tim-1)*T){ans=sa[i],len=T;break;}
for(ri i=ans,j=1;j<=len*tim;++j,++i)printf("%c",s[i]);
}
int main(){
for(ri i=2;i<=100000;++i)Log[i]=Log[i>>1]+1;
for(ri tt=1;;++tt,puts(""),top=0){
scanf("%s",s+1),n=strlen(s+1);
if(s[1]=='#')break;
printf("Case %d: ",tt),getsa(),query();
}
return 0;
}