[BZOJ2946][POI2000]公共串（后缀自动机）

题目：

我是超链接

题解：

后缀自动机的应用—-公共子串问题
对第一个串建立后缀自动机，不过这次匹配的时候是多个串进行匹配。每次匹配的时候对于后缀自动机中的每个节点维护一个值h，表示的是到达该节点所能匹配上的最大长度。光匹配还不够，我们需要按照拓扑序倒序，用每个节点取更新他的parent节点，因为如果匹配到一个状态，那么实际上他parent链上的所有状态都匹配上了。然后对于每个串匹配后得到的每个位置的h数组取min得到数组g（这是对于所有串每个位置的横向比较，长的字符串可以将就短的，但是短的不能将就长的啊），g数组再取最大值就是一个纵向比较，这样一横一纵，对于所有位置的所有字符串就可以得到比较了。
PS：所谓的拓扑倒序就是按照step从大到小。

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int N=4005;
char st[N];
int p,np,q,nq,cnt,last,ch[N*2][30],fa[N*2],c[N],a[N],step[N],cl[N],mn[N];
void insert(int c)
{
    p=last; np=last=++cnt;
    step[np]=step[p]+1;
    while (p && !ch[p][c]) ch[p][c]=np,p=fa[p];
    if (p==0){fa[np]=1;return;}
    q=ch[p][c];
    if (step[q]==step[p]+1){fa[np]=q;return;}
    nq=++cnt; step[nq]=step[p]+1;
    memcpy(ch[nq],ch[q],sizeof(ch[q]));
    fa[nq]=fa[q]; fa[q]=fa[np]=nq; 
    while (ch[p][c]==q) ch[p][c]=nq,p=fa[p];
}
void solve()
{
    memset(cl,0,sizeof(cl));
    int l=strlen(st),tmp=0;p=1;这里的p一定要赋在起点
    for (int i=0;i<l;i++)
    {
        int c=st[i]-'a';
        if (ch[p][c]) p=ch[p][c],tmp++;
        else
        {
            while (p && !ch[p][c]) p=fa[p];
            if (!p) p=1,tmp=0;
            else tmp=step[p]+1,p=ch[p][c];
        }
        cl[p]=max(cl[p],tmp);
    }
    for (int i=cnt;i;i--)
    {
        int t=a[i];
        mn[t]=min(mn[t],cl[t]);
        if (fa[t] && cl[t]) cl[fa[t]]=step[fa[t]];//如果后面的可以匹配上那前面的一定可以匹配上啦，而且是最长方案哦 
    }
}
int main()
{
    int n;scanf("%d",&n);
    memset(mn,0x7f,sizeof(mn));
    for (int i=1;i<=n;i++)
    {
        scanf("%s",st);
        if (i==1)
        {
            last=cnt=1;
            int len=strlen(st);
            for (int j=0;j<len;j++) insert(st[j]-'a');

            for (int j=1;j<=cnt;j++) c[step[j]]++;
            for (int j=1;j<=cnt;j++) c[j]+=c[j-1];
            for (int j=1;j<=cnt;j++) a[c[step[j]]--]=j;
        }
        else solve();
    } 
    int ans=0;
    for (int i=1;i<=cnt;i++) ans=max(ans,mn[i]);
    printf("%d",ans);
}