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

任意选择一个串作为模式串，构建出后缀自动机。

然后用其他的串在后缀自动机上跑匹配。

然后就到了理解后缀自动机性质的时候。

在某一个节点的最大值是可以沿着parent树上传的。

然后用dp[i][j]表示字符串i在位置j上的匹配的最长长度。

然后上传即可。之后在每一个节点统计所有串的min，然后再取max即可。

写法很暴力，貌似已经接近n^2了，但是还是A掉了。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i<=k;--i)
#define inf 0x3f3f3f3f
#define ll long long
#define maxn 5005
 
struct sam{
    char s[maxn];
    int dp[maxn][10],n;
    int last,cnt,len,lth;
    int go[maxn][26],l[maxn],fa[maxn];
    void init()
    {
        last=cnt=1;
        memset(go,0,sizeof go);
        memset(dp,0,sizeof dp);
    }
    void add(int x)
    {
//      printf("add %d
",x);
        int p=last,np=last=++cnt; l[np]=l[p]+1;
//      printf("on point %d
",np);
        for (;p&&!go[p][x];p=fa[p]) go[p][x]=np;
        if (!p) fa[np]=1;
        else
        {
            int q=go[p][x];
            if (l[q]==l[p]+1) fa[np]=q;
            else
            {
                int nq=++cnt;
                l[nq]=l[p]+1;
                memcpy(go[nq],go[q],sizeof go[q]);
                fa[nq]=fa[q];
                fa[np]=fa[q]=nq;
                for (;p&&go[p][x]==q;p=fa[p]) go[p][x]=nq;
            }
        }
    }
    void solve()
    {
        scanf("%d",&n);
        scanf("%s",s+1);
        lth=len=strlen(s+1);
        F(i,1,lth) add(s[i]-'a');
        F(i,2,n)
        {
            scanf("%s",s+1);
            len=strlen(s+1);
            int now=1,t=0;
            F(j,1,len)
            {
                int x=s[j]-'a';
                if (go[now][x]) {t++;now=go[now][x];}
                else
                {
                    while (now&&!go[now][x]) now=fa[now];
                    if (!now) {t=0;now=1;}
                    else t=l[now]+1,now=go[now][x];
                }
                for (int k=now;k;k=fa[k]) dp[k][i]=max(dp[k][i],t);
            }
        }
        int ans=0;
        F(i,1,cnt)
        {
            int tmp=l[i];
            F(j,2,n)
                tmp=min(tmp,dp[i][j]);
            ans=max(ans,tmp);
        }
        printf("%d
",ans);
    }
}SAM;
 
int main()
{
    SAM.init();
    SAM.solve();
}