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;
    }
}
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容