BZOJ 2946 [Poi2000]公共串

最新推荐文章于 2022-02-24 16:25:58 发布

magic_sheep

最新推荐文章于 2022-02-24 16:25:58 发布

阅读量588

点赞数

CC 4.0 BY-SA版权

分类专栏：后缀自动机文章标签：后缀自动机 bzoj POI

本文链接：https://blog.youkuaiyun.com/magic_sheep/article/details/62443134

后缀自动机专栏收录该内容

6 篇文章

订阅专栏

本文介绍了一种利用后缀自动机进行字符串匹配的方法。首先为首个字符串建立后缀自动机，随后在该自动机中匹配其他字符串。通过记录状态最小匹配值并在所有状态中寻找最大值来解决特定问题。文章提供了完整的C++代码实现。

第一个串建立后缀自动机，其余串在后缀自动机上进行匹配，每个状态记录所有串匹配最小值，最后对所有状态求最大值即可。

对mx[]排序可以得到后缀树和自动机的拓扑序。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=4010;
int n,m;
char ch[3005];
struct SAM
{
    int cnt,last;
    int fa[maxn],a[maxn][27],mx[maxn],len[maxn],ans[maxn];
    int v[maxn],q[maxn];
    SAM(){
        last=++cnt;
    }
    void extend(int c){
        int p=last,np=last=++cnt;mx[np]=mx[p]+1;
        while(!a[p][c]&&p)a[p][c]=np,p=fa[p];
        if(!p)fa[np]=1;
        else
        {
            int q=a[p][c];
            if(mx[p]+1==mx[q])fa[np]=q;
            else
            {
                int nq=++cnt;mx[nq]=mx[p]+1;
                memcpy(a[nq],a[q],sizeof(a[q]));
                fa[nq]=fa[q];
                fa[np]=fa[q]=nq;
                while(a[p][c]==q)a[p][c]=nq,p=fa[p];
            }
        } 
    }
    void pre(){
        for(int i=1;i<=cnt;i++)ans[i]=mx[i];
        for(int i=1;i<=cnt;i++)v[mx[i]]++;
        for(int i=1;i<=cnt;i++)v[i]+=v[i-1];
        for(int i=cnt;i;i--)q[v[mx[i]]--]=i; 
    }
    void solve(){
        scanf("%s",ch+1);
        memset(len,0,sizeof(len));
        int l=strlen(ch+1),p=1,tmp=0;
        for(int i=1;i<=l;i++)
        {
            int c=ch[i]-'a';
            while(!a[p][c]&&p)p=fa[p];
            if(p==0)p=1,tmp=0;
            else tmp=min(tmp,mx[p])+1,p=a[p][c];
            len[p]=max(len[p],tmp);
        }
        for(int i=cnt;i;i--)len[fa[q[i]]]=max(len[fa[q[i]]],len[q[i]]);
        for(int i=1;i<=cnt;i++)ans[i]=min(ans[i],len[i]);
    }
}SAM;
int main()
{
    scanf("%d",&n);
    scanf("%s",ch+1);
    int l=strlen(ch+1);
    for(int i=1;i<=l;i++)SAM.extend(ch[i]-'a');
    SAM.pre();
    for(int i=1;i<n;i++)SAM.solve();
    int ans=0;
    for(int i=1;i<=SAM.cnt;i++)
        ans=max(ans,SAM.ans[i]);
    printf("%d\n",ans);
    return 0;   
}