uoj#35 后缀排序（后缀数组模版）

 #include<bits/stdc++.h>
 #define N 100005
 using namespace std;
 char s[N];
 int a[N],c[N],t1[N],t2[N],sa[N],rk[N],ht[N];
 int m,n,p;
 void calcsa(int n,int m){
     int *x=t1,*y=t2,f=,p=;
     for(int i=;i<=m;i++)c[i]=;
     for(int i=;i<=n;i++)c[x[i]=a[i]]++;
     for(int i=;i<=m;i++)c[i]+=c[i-];
     for(int i=n;i>=;i--)sa[c[x[i]]--]=i;
     for(int i=;i<=n&&p<=n;i<<=){p=;
         for(int j=n-i+;j<=n;j++)y[++p]=j;
         for(int j=;j<=n;j++)if(sa[j]>i)y[++p]=sa[j]-i;
         for(int j=;j<=m;j++)c[j]=;
         for(int j=;j<=n;j++)c[x[y[j]]]++;
         for(int j=;j<=m;j++)c[j]+=c[j-];
         for(int j=n;j>=;j--)sa[c[x[y[j]]]--]=y[j];
         swap(x,y);x[sa[]]=;p=;
         for(int j=;j<=n;j++)
             x[sa[j]]=y[sa[j]]==y[sa[j-]]&&y[sa[j]+i]==y[sa[j-]+i]?p-:p++;
         m=p;
     }
     for(int i=;i<=n;i++)rk[sa[i]]=i;
     for(int i=;i<=n;i++){
         int j=sa[rk[i]-];
         if(f)f--;while(a[i+f]==a[j+f])f++;
         ht[rk[i]]=f;
     }
 }
 int main(){
     scanf("%s",s);int len=strlen(s);
     for(int i=;i<len;i++)a[++n]=s[i]-'a'+;
     calcsa(n,);
     for(int i=;i<=n;i++)printf("%d ",sa[i]);puts("");
     for(int i=;i<=n;i++)printf("%d ",ht[i]);puts("");
     return ;
 }