后缀数组学习笔记

现在来看倍增算法是非常好理解的。

直接放一篇blog写的挺好的：http://www.cnblogs.com/zinthos/p/3899725.html

虽然理论复杂度是$O(nlogn)$，但其中各种细节优化确实十分有必要的。

给自己放一个倍增的模板，有空填DC3的坑

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m;
char s[1000010];
int rk[1000010],sa[1000010],tmp[1000010],c[1000010];
void outint(int x){
    if(x>=10) outint(x/10);
    putchar(x%10+'0');
}
bool inline cmp(int *r,int a,int b,int l){
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void Getsa(){
    m=127;
    for(int i=0;i<=m;i++) c[i]=0;
    for(int i=1;i<=n;i++) c[rk[i]=s[i]]++;
    for(int i=1;i<=m;i++) c[i]+=c[i-1];
    for(int i=n;i>=1;i--) sa[c[rk[i]]--]=i;
    for(int i=1;i<=n;i<<=1){
        int p=0;
        for(int j=n-i+1;j<=n;j++) tmp[++p]=j;
        for(int j=1;j<=n;j++) if(sa[j]>i) tmp[++p]=sa[j]-i;
        for(int j=0;j<=m;j++) c[j]=0;
        for(int j=1;j<=n;j++) c[rk[tmp[j]]]++;
        for(int j=1;j<=m;j++) c[j]+=c[j-1];
        for(int j=n;j>=1;j--) sa[c[rk[tmp[j]]]--]=tmp[j];
        swap(rk,tmp);
        rk[sa[1]]=p=1;
        for(int j=2;j<=n;j++) rk[sa[j]]=cmp(tmp,sa[j],sa[j-1],i)?p:++p;
        if(p>=n) break;
        m=p;
    }
}
int main(){
    scanf("%s",s+1);
    n=strlen(s+1);
    Getsa();
    for(int i=1;i<=n;i++){
        outint(sa[i]);
        putchar(' ');
    }
    return 0;
}