[POI2000 bzoj2946]公共串

题目描述：给出几个由小写字母构成的单词，求它们最长的公共子串的长度。
任务：
l 读入单词
l 计算最长公共子串的长度
l 输出结果

分析：
1.第一个串放进去建后缀自动机
2.以后的每个串依次匹配，看到达每个状态的最短长度
3.取所有状态的最大值

#include<cstdio>#include<algorithm>#include<cstring>#include<cmath>#include<iostream>using namespace std;const int maxn=4010;int n;char str[maxn];int last;int e;int len[maxn];void read(int& x){    x=0;    char c=getchar();    while(c<'0' || c>'9') c=getchar();    while(c>='0' && c<='9'){        x=x*10+c-'0';        c=getchar();    }}struct Suffix_Automaton{    int ch[maxn][26],pre[maxn],Max[maxn];    int cnt[maxn],rk[maxn];    int ans[maxn];    Suffix_Automaton(){        last=++e;    }    void extand(int c){        int p=last,np=++e;Max[np]=Max[p]+1;last=np;        for(;p && !ch[p][c];p=pre[p]) ch[p][c]=np;        if(!p) pre[np]=1;        else{            int q=ch[p][c];            if(Max[q]!=Max[p]+1){                int nq=++e;Max[nq]=Max[p]+1;                for(int i=0;i<26;i++) ch[nq][i]=ch[q][i];                pre[nq]=pre[q];                pre[q]=pre[np]=nq;                for(;ch[p][c]==q;p=pre[p]) ch[p][c]=nq;            }else pre[np]=q;        }    }    void Pre(){        for(int i=1;i<=e;i++) ans[i]=Max[i];        for(int i=1;i<=e;i++) cnt[Max[i]]++;        for(int i=1;i<=e;i++) cnt[i]+=cnt[i-1];        for(int i=e;i>=1;i--) rk[cnt[Max[i]]--]=i;    }    void solve(){        scanf("%s",str+1);        memset(len,0,sizeof(len));        int p=1,tmp=0;        for(int i=1;i<=strlen(str+1);i++){            int c=str[i]-'a';            for(;!ch[p][c] && p;p=pre[p]);            if(!p) p=1,tmp=0;            else tmp=min(tmp,Max[p])+1,p=ch[p][c];            len[p]=max(len[p],tmp);        }        for(int i=e;i>=1;i--) len[pre[rk[i]]]=max(len[pre[rk[i]]],len[rk[i]]);        for(int i=1;i<=e;i++) ans[i]=min(ans[i],len[i]);    }}sam;int ans;int main(){    read(n);    scanf("%s",str+1);    for(int i=1;i<=strlen(str+1);i++) sam.extand(str[i]-'a');    sam.Pre();    for(int i=1;i<n;i++) sam.solve();    for(int i=1;i<=e;i++) ans=max(ans,sam.ans[i]);    printf("%d",ans);    return 0;}

^_^