[bzoj2946]公共串

题目大意

求n个字符串的最长公共子串

SAM

对第一个字符串建SAM
每个节点为护f[i]表示所有串与该状态匹配的最长长度。
对于每一个字符串，都在SAM上跑，并每次都计算g[i]表示字符串与状态i匹配的最大长度。
做完后用g和f取min更新f。
每次跑完后，要逆推计算g。对于i，若g[i]>0，则g[pre[i]]=step[pre[i]]。因为i成功被匹配了，匹配长度>step[pre[i]]，pre[i]也被匹配了，一定是step[pre[i]]被匹配掉了。
最后在f中取最大值。

#include<cstdio>#include<algorithm>#include<cstring>#define fo(i,a,b) for(i=a;i<=b;i++)#define fd(i,a,b) for(i=a;i>=b;i--)using namespace std;const int maxn=2000+10;int pre[maxn*2],step[maxn*2],g[maxn*2][26],cnt[maxn*2],a[maxn*2];int f[maxn*2],b[maxn*2];char s[maxn];int i,j,k,l,t,n,m,p,x,tot,last,len,ans;void add(int x){    int np=++tot,p=last;    step[np]=step[p]+1;    while (p&&g[p][x]==0){        g[p][x]=np;        p=pre[p];    }    if (!p) pre[np]=1;    else{        int q=g[p][x];        if (step[q]==step[p]+1) pre[np]=q;        else{            int nq=++tot,i;            fo(i,0,25) g[nq][i]=g[q][i];            step[nq]=step[p]+1;            pre[nq]=pre[q];            pre[q]=nq;            pre[np]=nq;            while (p&&g[p][x]==q){                g[p][x]=nq;                p=pre[p];            }        }    }    last=np;}int main(){    scanf("%d",&n);    scanf("%s",s+1);    last=tot=1;    len=strlen(s+1);    fo(i,1,len) add(s[i]-'a');    fo(i,1,tot) f[i]=step[i];    fo(i,1,tot) cnt[step[i]]++;    fo(i,1,len) cnt[i]+=cnt[i-1];    fo(i,1,tot) a[cnt[step[i]]--]=i;    n--;    while (n--){        fo(i,1,tot) b[i]=0;        scanf("%s",s+1);        len=strlen(s+1);        p=1;t=0;        fo(i,1,len){            x=s[i]-'a';            while (p!=1&&!g[p][x]) p=pre[p];            if (!g[p][x]) t=0;else t=min(t,step[p])+1;            if (g[p][x]) p=g[p][x];            b[p]=max(b[p],t);        }        fd(i,tot,1)             if (b[a[i]]) b[pre[a[i]]]=step[pre[a[i]]];        fo(i,1,tot) f[i]=min(f[i],b[i]);    }    ans=0;    fo(i,1,tot)        ans=max(ans,f[i]);    printf("%d\n",ans);}