模板_后缀数组

#define N 200005
int t1
,t2
,c
;
int sa
,rk
,ht
;
int n,m=27,len;int s
;
void sort(int *x,int *y){
for(int i=0;i<m;i++) c[i]=0;
for(int i=0;i<n;i++) c[x[y[i]]]++;
for(int i=1;i<m;i++) c[i]+=c[i-1];
for(int i=n-1;~i;--i) sa[--c[x[y[i]]]]=y[i];
}
void DA(){
int *x=t1,*y=t2;
for(int i=0;i<m;i++) c[i]=0;
for(int i=0;i<n;i++) c[x[i]=s[i]]++;
for(int i=1;i<m;i++) c[i]+=c[i-1];
for(int i=n-1;~i;i--) sa[--c[x[i]]]=i;
for(int k=1,p=0;k<=n;k<<=1,p=0){
for(int i=n-k;i<n;i++) y[p++]=i;
for(int i=0;i<n;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
sort(x,y);swap(x,y);m=1,x[sa[0]]=0;
for(int i=1;i<n;i++) x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&
((sa[i]+k>=n?-1:y[sa[i]+k])==(sa[i-1]+k>=n?-1:y[sa[i-1]+k]))?m-1:m++;
if(m>n) break;
}
for(int i=0;i<n;i++) rk[sa[i]]=i;
for(int i=0,j,k=0;i<n;ht[rk[i++]]=k){
if(!rk[i]) continue;
j=sa[rk[i]-1],k=k?k-1:0;
while(i+k<n&&j+k<n&&s[i+k]==s[j+k]) k++;
}ht[0]=0;
}