2946: [Poi2000]公共串 后缀自动机

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/Phenix_2015/article/details/50958582

二分答案+hash
二分答案+后缀数组
花3s感受一下:

说一下后缀自动机的做法。

我们对一个串建立后缀自动机,而由于后缀自动机中由S到每个节点的路径能表示出以此节点结束的后缀,所以我们可以把每个串放进去匹配,得到每个串在每个状态的最大匹配长度max,并且要沿parent树更新父节点的max,最后将每个状态各个串的最大匹配长度取最小值,最后求各最小值的最大值。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 10005
using namespace std;
char s[N];
int n,cnt,last,S;
int fa[N],ch[N][26],len[N],now[N],ans[N],q[N],cc[N];
inline void insert(int c)
{
    int p=last,np=++cnt; last=np;
    len[np]=len[p]+1;
    while (p&&!ch[p][c]) ch[p][c]=np,p=fa[p];
    if (!p) fa[np]=S;
    else
    {
        int q=ch[p][c];
        if (len[p]+1==len[q]) fa[np]=q;
        else
        {
            int nq=++cnt; len[nq]=len[p]+1;
            memcpy(ch[nq],ch[q],sizeof(ch[q]));
            fa[nq]=fa[q]; fa[np]=fa[q]=nq;
            while (ch[p][c]==q) ch[p][c]=nq,p=fa[p];
        }
    }
}
inline void solve()
{
    scanf("%s",s+1);
    int l=strlen(s+1);
    int p=1,tmp=0;
    memset(now,0,sizeof(now));
    for (int i=1;i<=l;i++)
    {
        int c=s[i]-'a';
        while (p&&!ch[p][c]) p=fa[p];
        if (!p) p=1,tmp=0;
        else tmp=min(tmp,len[p])+1,p=ch[p][c];
        now[p]=max(now[p],tmp);
    }
    for (int i=cnt;i;i--) now[fa[q[i]]]=max(now[fa[q[i]]],now[q[i]]);
    for (int i=1;i<=cnt;i++) ans[i]=min(ans[i],now[i]);
}
int main()
{
    scanf("%d",&n);
    scanf("%s",s+1);
    int l=strlen(s+1);
    S=last=++cnt;
    for (int i=1;i<=l;i++)
        insert(s[i]-'a');
    for (int i=1;i<=cnt;i++) ans[i]=len[i];
    for (int i=1;i<=cnt;i++) ++cc[len[i]];
    for (int i=1;i<=cnt;i++) cc[i]+=cc[i-1];
    for (int i=cnt;i;i--) q[cc[len[i]]--]=i;
    for (int i=1;i<n;i++) solve();
    int Ans=0;
    for (int i=1;i<=cnt;i++) 
        Ans=max(Ans,ans[i]);
    cout << Ans << endl;
}
阅读更多
换一批

没有更多推荐了,返回首页