• 【模板】后缀排序


    传送门

    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!

  • 相关阅读:
    LeetCode 338. 比特位计数
    LeetCode 208. 实现 Trie (前缀树)
    初识restful api接口
    破解 Navicat Premium 12
    ES6 Reflect的认识
    ES6 WeakMap和WeakSet的使用场景
    sublime 注释模版插件DocBlockr的使用
    js call方法的使用
    ES6 Generator的应用场景
    ES6 Symbol的应用场景
  • 原文地址:https://www.cnblogs.com/PaperCloud/p/10197729.html
Copyright © 2020-2023  润新知