模板：后缀数组

http://www.cnblogs.com/zwfymqz/p/8413523.html#_label5

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 1000005
using namespace std;
int n,m,sa[maxn],rk[maxn],tp[maxn],tax[maxn],p;
char s[maxn];
void Qsort(){
    for(int i=0;i<=m;i++)tax[i]=0;
    for(int i=1;i<=n;i++)tax[rk[i]]++;
    for(int i=1;i<=m;i++)tax[i]+=tax[i-1];
    for(int i=n;i>=1;i--)sa[ tax[rk[tp[i]]]-- ]=tp[i];
}

int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    for(int i=1;i<=n;i++)rk[i]=s[i]-'0',tp[i]=i;
    m=n;Qsort();
    for(int ws=1,p=0;p<n;m=p,ws<<=1){
        p=0;
        for(int i=1;i<=ws;i++)tp[++p]=n-ws+i;
        for(int i=1;i<=n;i++)if(sa[i]>ws)tp[++p]=sa[i]-ws;
        Qsort();
        swap(tp,rk);
        rk[sa[1]]=p=1;
        for(int i=2;i<=n;i++){
            rk[sa[i]]=(tp[sa[i-1]]==tp[sa[i]]&&tp[sa[i-1]+ws]==tp[sa[i]+ws])?p:++p;
        }
    }
    cout<<sa[1];
    for(int i=2;i<=n;i++)printf(" %d",sa[i]);
    cout<<endl;
    //cout<<(int)'a'<<' '<<(int)'z'<<' '<<(int)'0'<<' '<<(int)'9'<<'
';
    return 0;
}
 //aabaaaab