【POI2000】【BZOJ2946】公共串

2946: [Poi2000]公共串
Time Limit: 3 Sec Memory Limit: 128 MB
Submit: 286 Solved: 121
[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

Source

1A就是爽OwO~(≧▽≦)/~
就是求多个串的最长公共子串.像我在后缀自动机的性质和应用里写的一样做就好了.我就不再写一遍了>_<反正很水…就是裸的SAM嘛= =
比起两个串的时候多开一个每个状态对某个串的最大匹配长度,每个状态对所有串的最大匹配长度的最小值.
然后把状态按其代表的最大子串长度用一下基数排序存在一个temp里面.

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#define MAXN 4010#define MAXINT 0x7fffffffusing namespace std;int Ans;int n;char ch[MAXN];struct sam{    int last,cnt,p,q,np,nq;    int a[MAXN][26],len[MAXN],fa[MAXN];//len每个状态代表的最长子串长度     int minn[MAXN],maxn[MAXN];//maxn在某串匹配过程中每个状态对该串的最长匹配长度    //minn每个状态对所有串的最长匹配长度的最小值     int count[MAXN],temp[MAXN];//把每个状态按其代表的最大子串长度排序,temp为排序后状态号.     sam()    {        last = ++cnt;    }    void insert(int c)    {        p=last;np=last=++cnt;len[np]=len[p]+1;        while (!a[p][c]&&p) a[p][c]=np,p=fa[p];        if (!p) fa[np]=1;        else        {            q=a[p][c];            if (len[q]==len[p]+1)   fa[np]=q;            else            {                nq=++cnt;len[nq]=len[p]+1;                memcpy(a[nq],a[q],sizeof(a[q]));                fa[nq]=fa[q];                fa[q]=fa[np]=nq;                while (a[p][c]==q)  a[p][c]=nq,p=fa[p];            }        }    }    void build()    {        scanf("%s",ch+1);        int l=strlen(ch+1);        for (int i=1;i<=l;i++)  insert(ch[i]-'a');    }    void init()    {        for (int i=1;i<=cnt;i++)    minn[i]=len[i];        for (int i=1;i<=cnt;i++)    count[len[i]]++;        for (int i=1;i<=cnt;i++)    count[i]+=count[i-1];        for (int i=cnt;i;i--)   temp[count[len[i]]--]=i;    }    void solve()    {        for (int T=1;T<n;T++)        {            scanf("%s",ch+1);            int l=strlen(ch+1),tmp=0,now=1;            memset(maxn,0,sizeof(maxn));            for (int i=1;i<=l;i++)            {                int c=ch[i]-'a';                while (!a[now][c]&&now) now=fa[now];                if (!now)   tmp=0,now=1;                else    tmp=min(tmp,len[now])+1,now=a[now][c];                maxn[now]=max(maxn[now],tmp);//能匹配到的最大串长             }            for (int i=cnt;i;i--)   maxn[fa[temp[i]]]=max(maxn[fa[temp[i]]],maxn[temp[i]]);            for (int i=1;i<=cnt;i++)    minn[i]=min(minn[i],maxn[i]);        }    }}sam;int main(){    scanf("%d",&n);    sam.build();    sam.init();    sam.solve();    for (int i=1;i<=sam.cnt;i++)    Ans=max(Ans,sam.minn[i]);    printf("%d",Ans);}