模板_后缀数组

倍增法

#define N 200005
int t1[N],t2[N],c[N];
int sa[N],rk[N],ht[N];
int n,m=27,len;int s[N];
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;
}