zoukankan      html  css  js  c++  java
  • 【模板】后缀排序

    传送门

    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!

  • 相关阅读:
    bzoj 1013: [JSOI2008]球形空间产生器sphere
    bzoj 1012: [JSOI2008]最大数maxnumber
    bzoj 1010: [HNOI2008]玩具装箱toy
    bzoj 1008: [HNOI2008]越狱
    bzoj 1007: [HNOI2008]水平可见直线
    UVa12105 越大越好
    POJ
    最优点配对问题(紫书)
    Ned 的难题
    UVA
  • 原文地址:https://www.cnblogs.com/PaperCloud/p/10197729.html
Copyright © 2011-2022 走看看