【bzoj2946】[Poi2000]公共串 后缀自动机

[Poi2000]公共串

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 1386  Solved: 620
[Submit][Status][Discuss]

Description

 

       给出几个由小写字母构成的单词，求它们最长的公共子串的长度。

任务：

l        读入单词

l        计算最长公共子串的长度

l        输出结果

 

Input

 

文件的第一行是整数 n，1<=n<=5，表示单词的数量。接下来n行每行一个单词，只由小写字母组成，单词的长度至少为1，最大为2000。

 

Output

仅一行，一个整数，最长公共子串的长度。

 

Sample Input

3
abcb
bca
acbc

Sample Output

２

HINT

后缀自动机

 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>

 #define N 2007
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
     return x*f;
 }

 int n;
 char ch[N];

 struct SAM
 {
     int cnt,last;
     int fa[],a[][],mx[],len[],ans[];
     int v[],q[];
     SAM(){last=++cnt;}
     void extend(int c)
     {
         int p=last,np=last=++cnt;mx[np]=mx[p]+;
         while(!a[p][c]&&p)a[p][c]=np,p=fa[p];
         if(!p)fa[np]=;
         else
         {
             int q=a[p][c];
             if(mx[p]+==mx[q])fa[np]=q;
             else
             {
                 int nq=++cnt;mx[nq]=mx[p]+;
                 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=;i<=cnt;i++)ans[i]=mx[i];
         for(int i=;i<=cnt;i++)v[mx[i]]++;
         for(int i=;i<=cnt;i++)v[i]+=v[i-];
         for(int i=cnt;i;i--)q[v[mx[i]]--]=i;
     }
     void solve()
     {
         scanf("%s",ch+);
         memset(len,,sizeof(len));
         int l=strlen(ch+),p=,tmp=;
         for(int i=;i<=l;i++)
         {
             int c=ch[i]-'a';
             while(!a[p][c]&&p)p=fa[p];
             if(p==)p=,tmp=;
             else tmp=min(tmp,mx[p])+,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=;i<=cnt;i++)ans[i]=min(ans[i],len[i]);
     }
 }sam;
 int main()
 {
     n=read();scanf("%s",ch+);
     int l=strlen(ch+);
     for(int i=;i<=l;i++)sam.extend(ch[i]-'a');
     sam.pre();
     for(int i=;i<n;i++)sam.solve();
     int ans=;
     for(int i=;i<=sam.cnt;i++)
         ans=max(ans,sam.ans[i]);
     printf("%d\n",ans);
 }