【模板】后缀排序

传送门

Solution

怎么写都要把自己绕进去的样子?


Code 

#include<bits/stdc++.h>
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x*f;
}
#define MN 1000005
int n,sa[2][MN],rk[2][MN],num[MN],mx;
char s[MN];
inline void Sort_sa()
{
    register int K,i,p,q;mx=122;
    for(i=1;i<=n;++i) num[s[i]]++;
    for(i=1;i<=mx;++i) num[i]+=num[i-1];
    for(i=n;i>=1;--i) sa[1][num[s[i]]--]=i;
    for(i=1;i<=n;++i) rk[1][sa[1][i]]=rk[1][sa[1][i-1]]+(s[sa[1][i]]!=s[sa[1][i-1]]);
    
    for(p=0,q=1,K=1;K<=n;K<<=1,p^=1,q^=1)
    {
        if(mx==n) break;
        for(i=1;i<=n;++i) num[rk[q][sa[q][i]]]=i;
    //	for(i=1;i<=2;++i) printf("%d ",num[i]);puts("");
        
        for(i=n;i;--i) sa[q][i]>K?sa[p][num[rk[q][sa[q][i]-K]]--]=sa[q][i]-K:0;
        for(i=n-K+1;i<=n;++i) sa[p][num[rk[q][i]]--]=i;
        
        for(i=1;i<=n;++i)
            rk[p][sa[p][i]]=rk[p][sa[p][i-1]]+(rk[q][sa[p][i]]!=rk[q][sa[p][i-1]]||rk[q][sa[p][i]+K]!=rk[q][sa[p][i-1]+K]);
        mx=rk[p][sa[p][n]];
    }
    for(i=1;i<=n;++i) printf("%d ",sa[q][i]);
}
int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);Sort_sa();
    return 0;
}


Blog来自PaperCloud,未经允许,请勿转载,TKS!

原文地址:https://www.cnblogs.com/PaperCloud/p/10197729.html