*后缀自动机

最新推荐文章于 2024-10-05 19:23:25 发布

DCDCBigBig

最新推荐文章于 2024-10-05 19:23:25 发布

阅读量275

点赞数

分类专栏：算法-字符串文章标签：字符串后缀算法后缀自动机算法

本文链接：https://blog.csdn.net/DCDCBigBig/article/details/79814657

版权

算法-字符串专栏收录该内容

8 篇文章 0 订阅

订阅专栏

今天是2018/4/4，DCDCBigBig的第三十二篇博文

zrO WJMZBMR Orz
膜拜陈立杰……

后缀自动机

//SPOJ8222
//题目大意：
//给长度为 n 的字符串 S , 对任意的 L , 求长度为 L 的子串最多出现的次数 
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
int len,ans[500010];
char s[300001];
struct SAM{
    int now,tot,rt,last,son[500010][26],fa[500010],mx[500010],l1[500010],l2[500010],siz[500010];
    void init(){
        memset(son,0,sizeof(son));
        tot=rt=last=1;
        mx[1]=0;
    }
    void extend(int ch){
        int p=last,np=++tot;
        mx[np]=mx[p]+1;
        for(;p&&!son[p][ch];p=fa[p])son[p][ch]=np;
        if(!p)fa[np]=rt;
        else{
            int q=son[p][ch];
            if(mx[q]==mx[p]+1)fa[np]=q;
            else{
                int nq=++tot;
                mx[nq]=mx[p]+1;
                memcpy(son[nq],son[q],sizeof(son[q]));
                fa[nq]=fa[q];
                fa[q]=fa[np]=nq;
                for(;son[p][ch]==q;p=fa[p])son[p][ch]=nq;
            }
        }
        last=np;
        siz[np]=1;
    }
    void solve(){
        for(int i=1;i<=tot;i++)l1[mx[i]]++;
        for(int i=2;i<=len;i++)l1[i]+=l1[i-1];
        for(int i=1;i<=tot;i++)l2[l1[mx[i]]--]=i;
        for(int i=tot;i;i--){
            siz[fa[l2[i]]]+=siz[l2[i]];
            ans[mx[l2[i]]]=max(ans[mx[l2[i]]],siz[l2[i]]);
        }
    }
}sam;
int main(){
    memset(ans,0,sizeof(ans));
    scanf("%s",s);
    sam.init();
    len=strlen(s);
    for(int i=0;i<len;i++){
        sam.extend(s[i]-'a');
    }
    sam.solve();
    for(int i=1;i<=len;i++){
        printf("%d\n",ans[i]);
    }
    return 0;
}

纯模板：

struct SAM{
    int now,tot,rt,last,son[500010][26],fa[500010],mx[500010],l1[500010],l2[500010],siz[500010];
    void init(){
        memset(son,0,sizeof(son));
        tot=rt=last=1;
        mx[1]=0;
    }
    void extend(int ch){
        int p=last,np=++tot;
        mx[np]=mx[p]+1;
        for(;p&&!son[p][ch];p=fa[p])son[p][ch]=np;
        if(!p)fa[np]=rt;
        else{
            int q=son[p][ch];
            if(mx[q]==mx[p]+1)fa[np]=q;
            else{
                int nq=++tot;
                mx[nq]=mx[p]+1;
                memcpy(son[nq],son[q],sizeof(son[q]));
                fa[nq]=fa[q];
                fa[q]=fa[np]=nq;
                for(;son[p][ch]==q;p=fa[p])son[p][ch]=nq;
            }
        }
        last=np;
        siz[np]=1;
    }
}sam;