后缀数组模板

最新推荐文章于 2021-08-14 21:36:30 发布

原创最新推荐文章于 2021-08-14 21:36:30 发布 · 223 阅读

0 ·

CC 4.0 BY-SA版权

模板同时被 2 个专栏收录

3 篇文章

订阅专栏

后缀数组

2 篇文章

订阅专栏

博客围绕后缀数组模板展开，虽无具体内容，但可推测会涉及后缀数组模板的相关知识，后缀数组在信息技术领域常用于字符串处理等方面。

#include<bits/stdc++.h>
const int maxn=100008;
int n,k;
int s[maxn];
int sa[maxn],rk[maxn],height[maxn],a[maxn],b[maxn],cnt[maxn];
inline void get_sa()
{
    int m=1000;
    int *rk_1=a,*sa_2=b,*tmp;
    for(int i=1;i<=n;++i)   ++cnt[rk_1[i]=s[i]];
    for(int i=1;i<=m;++i)   cnt[i]+=cnt[i-1];
    for(int i=n;i;--i)  sa[cnt[rk_1[i]]]=i,--cnt[rk_1[i]];
    
    for(int k=1,p;k<=n;k<<=1,m=p){
        p=0;
        for(int i=n-k+1;i<=n;++i)   sa_2[++p]=i;
        for(int i=1;i<=n;++i)   if(sa[i]>k) sa_2[++p]=sa[i]-k;
        
        for(int i=1;i<=m;++i)   cnt[i]=0;
        for(int i=1;i<=n;++i)   ++cnt[rk_1[sa_2[i]]];
        for(int i=1;i<=m;++i)   cnt[i]+=cnt[i-1];
        for(int i=n;i;--i)  sa[cnt[rk_1[sa_2[i]]]]=sa_2[i],--cnt[rk_1[sa_2[i]]];
        
        tmp=sa_2;
        p=1;
        tmp[sa[1]]=1;
        for(int i=2;i<=n;++i){
            if(sa[i]+k<=n&&sa[i-1]+k<=n&&rk_1[sa[i]]==rk_1[sa[i-1]]&&rk_1[sa[i]+k]==rk_1[sa[i-1]+k]){
                tmp[sa[i]]=p;
            }
            else{
                tmp[sa[i]]=++p;
            }
        }
        if(p==n)    break;
        swap(rk_1,sa_2);
    }
}
inline void get_height()
{
    int k=0;
    for(int i=1;i<=n;++i)   rk[sa[i]]=i;
    for(int i=1;i<=n;++i){
        if(rk[i]>1){
            if(k)   --k;
            int j=sa[rk[i]-1];
            while(i+k<=n&&j+k<=n&&s[i+k]==s[j+k]){
                ++k;
            }
            height[rk[i]]=k;
        }
    }
}