bzoj2946 [Poi2000]公共串（SAM）

Description

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

任务：

读入单词

计算最长公共子串的长度
输出结果

 

Input

 

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

 

Output

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

 

Sample Input

3

abcb

bca

acbc

Sample Output

２

[ Submit][ Status][ Discuss] 

分析：
这道题和两字符串的最长公共子串大同小异
我们还是用第一个字符串建立$SAM$
之后我们要用剩下的$n-1$个字符串在$SAM$上进行匹配

还是用之前的思路：两个字符串匹配到的最远结点就是最长公共子串
匹配的时候，记录在每一个结点（每一个状态）匹配到的最大长度（一个结点可能是多个结点的后继）

考虑如果一个结点代表的状态匹配到的最大长度是$l[i]$
在$parant$树上，父结点代表的状态一定是子结点的子串，那么父结点一定都匹配上了
因此对于状态$i$，我们需要用ta的答案更新$parent$的答案

我们要逆序更新祖先答案：$l[fa[i]]=dis[fa[i]]$
（简单来说就是按照$dis$数组从大到小）

每个结点取所有串的 l l <script type="math/tex" id="MathJax-Element-365">l</script>最小值，得出最长公共子串长
所有公共子串长度的最大值就是答案

#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int N=2003;
int n,dis[N<<1],ch[N<<1][26],last=1,sz=1,root=1,fa[N<<1];
int mn[N<<1],v[N<<1],pos[N<<1],l[N<<1],len;
char s[N];

void insert(int x)
{
    int now=++sz,pre=last;
    last=now; 
    dis[now]=dis[pre]+1;
    for (;pre&&!ch[pre][x];pre=fa[pre]) ch[pre][x]=now;

    if (!pre) fa[now]=root;   //原parent树中没有该结点 
    else
    {
        int q=ch[pre][x];     //now的直系祖先 
        if (dis[q]==dis[pre]+1) fa[now]=q;
        else
        {
            int nows=++sz;    //new node
            dis[nows]=dis[pre]+1;
            memcpy(ch[nows],ch[q],sizeof(ch[q]));
            fa[nows]=fa[q]; fa[q]=fa[now]=nows; 
            for (;pre&&ch[pre][x]==q;pre=fa[pre]) ch[pre][x]=nows;
        }
    }
}

void solve()
{
    memset(l,0,sizeof(l));   //每个结点能匹配的最大长度 
    int now=root,tmp=0;
    for (int i=1;i<=len;i++)
    {
        int x=s[i]-'a';
        if (ch[now][x]) now=ch[now][x],tmp++;
        else
        {
            while (now&&!ch[now][x]) now=fa[now];
            if (!now) now=root,tmp=0;
            else tmp=dis[now]+1,now=ch[now][x];
        }
        l[now]=max(l[now],tmp);
    }
    for (int i=sz;i>=1;i--)  //倒序 
    {
        int t=pos[i];
        mn[t]=min(mn[t],l[t]);
        if (fa[t]&&l[t]) l[fa[t]]=dis[fa[t]];
    }
}

int main()
{
    memset(mn,127,sizeof(mn));

    scanf("%d",&n);
    scanf("%s",s+1);
    len=strlen(s+1);
    for (int i=1;i<=len;i++) insert(s[i]-'a');

    if (n==1) 
    {
        printf("%d\n",len);
        return 0;
    }

    for (int i=1;i<=sz;i++) v[dis[i]]++;          //prepare
    for (int i=1;i<=len;i++) v[i]+=v[i-1];
    for (int i=1;i<=sz;i++) pos[v[dis[i]]--]=i;   //pos[i]：按照长度排名，从小到大 
    for (int i=2;i<=n;i++)
    {
        scanf("%s",s+1); len=strlen(s+1);
        solve();
    }

    int ans=0;
    for (int i=1;i<=sz;i++) ans=max(ans,mn[i]);
    printf("%d",ans);
    return 0;
}