Suffix Array

struct SuffixArray
{
    int id;ll rk;
    bool operator<(const SuffixArray &b)const
    {
        return rk<b.rk;
    }
}sa[__];

//rk[i]: s[i…n]的排名
//sa[i].id: 第i名为s[sa[i].id…n]
template<class T>
void get_sa(T s[],int n,int rk[])
{
    for(int i=1;i<=n;i++)
        sa[i].id=i,sa[i].rk=s[i];
    for(int i=1;;i<<=1)
    {
        sort(sa+1,sa+1+n);
        for(int j=1;j<=n;j++)
        {
            rk[sa[j].id]=rk[sa[j-1].id];
            if(sa[j].rk!=sa[j-1].rk)
                rk[sa[j].id]++;
        }
        if(rk[sa[n].id]==n)break;
        for(int j=1;j<=n;j++)
        {
            sa[j].id=j;
            sa[j].rk=1ll*rk[j]*__;
            if(j+i<=n)
                sa[j].rk+=rk[j+i];
        }
    }
}

//he[i]: 第i名后缀与第i-1名后缀的LCP
template<class T>
void get_he(T s[],int n,int rk[],int he[])
{
    for(int i=1,cp=0;i<=n;i++)
    {
        if(cp)--cp;
        int j=sa[rk[i]-1].id,maxx=max(i,j);
        while(maxx+cp<=n && s[i+cp]==s[j+cp])
            cp++;
        he[rk[i]]=cp;//h[i]=cp;
    }
}
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容