【BZOJ 2946】 2946: [Poi2000]公共串 （SAM）

2946: [Poi2000]公共串

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 1063  Solved: 469

Description

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

任务：

l        读入单词

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

l        输出结果

Input

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

Output

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

Sample Input

3
abcb
bca
acbc

Sample Output

２

HINT

Source

【分析】

　　重新学一次SAM，从刷水题开始。

　　同spoj1812。用第一个串建sam，然后用其他的串跑，ans存在那个点那里，做完一个串的时候还要根据parent边的拓扑序更新该点ans，最后取min即可。

　　当然后缀数组也是可以的。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define Maxn 2010

int mymax(int x,int y) {return x>y?x:y;}
int mymin(int x,int y) {return x<y?x:y;}

struct node
{
    int pre,son[30],step;
    // node() {pre=step=0;memset(son,0,sizeof(son));}
}t[Maxn*2];
int ans[Maxn*2],ad[Maxn*2];

struct sam
{
    int last,tot;
    /*void upd(int x)
    {
        memset(t[x].son,0,sizeof(t[x].son));
    }*/
    void extend(int k)
    {
        int np=++tot,p=last;
        t[np].step=t[last].step+1;
        while(p&&!t[p].son[k])
        {
            t[p].son[k]=np;
            p=t[p].pre;
        }
        if(!p) t[np].pre=1;
        else
        {
            int q=t[p].son[k];
            if(t[q].step==t[p].step+1) t[np].pre=q;
            else
            {
                int nq=++tot;//upd(tot);
                t[nq].step=t[p].step+1;
                memcpy(t[nq].son,t[q].son,sizeof(t[nq].son));
                t[nq].pre=t[q].pre;
                t[q].pre=t[np].pre=nq;
                while(p&&t[p].son[k]==q)
                {
                    t[p].son[k]=nq;
                    p=t[p].pre;
                }
            }
        }
        last=np;
    }
    
}sam;

char s[10],ss[Maxn];

int main()
{
    int n;scanf("%d",&n);
    scanf("%s",s);
    sam.last=sam.tot=1;
    int l=strlen(s);
    for(int i=0;i<l;i++) sam.extend(s[i]-'a'+1);
    memset(ans,63,sizeof(ans));
    for(int i=2;i<=n;i++)
    {
        scanf("%s",s);
        l=strlen(s);
        int nw=1,sp=0;
        for(int j=1;j<=sam.tot;j++) ad[j]=0;
        for(int j=0;j<l;j++)
        {
            int ind=s[j]-'a'+1;
            while(nw&&!t[nw].son[ind]) nw=t[nw].pre,sp=t[nw].step;
            if(t[nw].son[ind]) sp++,nw=t[nw].son[ind];
            else nw=1,sp=0;
            ad[nw]=mymax(ad[nw],sp);
        }
        for(int i=sam.tot;i>=1;i--) ad[t[i].pre]=mymax(ad[t[i].pre],mymin(ad[i],t[t[i].pre].step));
        for(int i=sam.tot;i>=1;i--) ans[i]=mymin(ans[i],ad[i]);
    }
    int mx=0;
    for(int i=1;i<=sam.tot;i++) if(ans[i]<=sam.tot) mx=mymax(mx,ans[i]);
    printf("%d
",mx);
    return 0;
}

2017-04-17 10:21:42