题目:
题解:
后缀自动机的应用—-公共子串问题
对第一个串建立后缀自动机,不过这次匹配的时候是多个串进行匹配。每次匹配的时候对于后缀自动机中的每个节点维护一个值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);
}