后缀数组

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define maxn 100005
int n,m,tot,t1[maxn],t2[maxn],Rank[maxn],sa[maxn],hei[maxn],sum[maxn];
char st[maxn];
int main(){
    scanf("%s",st+1),n=strlen(st+1);
    int *x=t1,*y=t2;
    m=256;
    for (int i=1;i<=m;i++) sum[i]=0;
    for (int i=1;i<=n;i++) sum[x[i]=st[i]]++;
    for (int i=1;i<=m;i++) sum[i]+=sum[i-1];
    for (int i=n;i>=1;i--) sa[sum[x[i]]--]=i;
    for (int len=1,tot=0;tot<n;len<<=1,m=tot){
        tot=0;
        for (int i=n-len+1;i<=n;i++) y[++tot]=i;
        for (int i=1;i<=n;i++) if (sa[i]>len) y[++tot]=sa[i]-len;
        for (int i=1;i<=m;i++) sum[i]=0;
        for (int i=1;i<=n;i++) sum[x[i]]++;
        for (int i=1;i<=m;i++) sum[i]+=sum[i-1];
        for (int i=n;i>=1;i--) sa[sum[x[y[i]]]--]=y[i];
        swap(x,y),tot=0;
        x[sa[1]]=++tot;
        for (int i=2;i<=n;i++){
            if (y[sa[i]]!=y[sa[i-1]]||y[sa[i]+len]!=y[sa[i-1]+len]) ++tot;
            x[sa[i]]=tot;
        }
    }
    for (int i=1;i<=n;i++) Rank[i]=x[i];
    memset(hei,0,sizeof(hei));
    for (int i=1,j=0;i<=n;i++){
        if (Rank[i]==1) continue;
        while (st[i+j]==st[sa[Rank[i]-1]+j]) j++;
        hei[Rank[i]]=j;
        if (j) j--;
    }
    for (int i=1;i<=n;i++) printf("%d
",sa[i]);
    return 0;
}

后缀数组模板。