这一定是我做的最恶心的一道后缀数组的题
简直恶心哭了
具体做法可以参见罗穗骞的论文
说一下字典序最小怎么找:
找出所有重复最多的子串长度
枚举sa数组最先找到的那个满足条件的 就是字典序最小的
ps:
1,st表查询的时候 左边的rank要加1!
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define T 100010
using namespace std;
char s[T];
int t1[T],t2[T],sa[T],cc[T],rank[T],height[T];
int st[T][17],z[T];
int len,ans,cnt,flag,num=0;
bool cmp(int *y,int a,int b,int k)
{
int a1=y[a],b1=y[b];
int a2=a+k>=len?-1:y[a+k];
int b2=b+k>=len?-1:y[b+k];
return a1==b1&&a2==b2;
}
void make_sa()
{
int *x=t1,*y=t2,m=333;
for(int i=0;i<m;i++)cc[i]=0;
for(int i=0;i<len;i++)++cc[x[i]=s[i]];
for(int i=1;i<m;i++)cc[i]+=cc[i-1];
for(int i=len-1;~i;i--)sa[--cc[x[i]]]=i;
for(int k=1;k<=len;k<<=1)
{
int p=0;
for(int i=len-k;i<len;i++)y[p++]=i;
for(int i=0;i<len;i++)if(sa[i]>=k)y[p++]=sa[i]-k;
for(int i=0;i<m;i++)cc[i]=0;
for(int i=0;i<len;i++)++cc[x[y[i]]];
for(int i=1;i<m;i++)cc[i]+=cc[i-1];
for(int i=len-1;~i;i--)sa[--cc[x[y[i]]]]=y[i];
swap(x,y);m=1;x[sa[0]]=0;
for(int i=1;i<len;i++)x[sa[i]]=cmp(y,sa[i],sa[i-1],k)?m-1:m++;
if(m>=len)return;
}
}
void make_height()
{
int k=0;
for(int i=0;i<len;i++)rank[sa[i]]=i;
for(int i=0;i<len;i++)
{
if(!rank[i])continue;
int j=sa[rank[i]-1];
if(k)k--;
while(s[i+k]==s[j+k])k++;
height[rank[i]]=k;
}
}
void make_st()
{
int w=log2(len);
st[0][0]=T;
for(int i=1;i<len;i++)st[i][0]=height[i];//cout << i <<" "<<height[i]<<endl;
for(int k=1;k<=w;k++)
for(int i=0;i<len;i++)
{
if(i+(1<<k)>=len)break;
st[i][k]=min(st[i][k-1],st[i+(1<<k-1)][k-1]);
}
}
int ask(int l,int r)
{
if(l>r)swap(l,r);l++;
int k=log2(r-l+1);
return min(st[l][k],st[r-(1<<k)+1][k]);
}
int main()
{
while(scanf("%s",s)&&s[0]!='#')
{
len=strlen(s);
make_sa();
make_height();
make_st();
ans=cnt=0;
for(int l=1;l<=len;l++)
for(int i=l;i<len;i+=l)
{
int k=ask(rank[i],rank[i-l]);
int m=l-k%l;
if(i>l&&k%l)k=max(k,ask(rank[i-l-m],rank[i-m]));
k=k/l+1;
if(k>ans)ans=k,cnt=1,z[1]=l;
else if(k==ans&&z[cnt]!=l)z[++cnt]=l;
}
flag=1;
for(int i=0;i<len&&flag;i++)
for(int j=1;j<=cnt&&flag;j++)
if(ask(rank[sa[i]],rank[sa[i]+z[j]]) >= (ans-1)*z[j])
{
//if(sa[i]+z[j]*ans>=len)continue;
flag=0;
printf("Case %d: ",++num);
s[sa[i]+z[j]*ans]='\0';
printf("%s\n",s+sa[i]);
}
}
}