【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 0x7fffffff
using 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);
}