[后缀自动机 模板题] SPOJ 1812Longest Common Substring II & BZOJ 2946 [Poi2000]公共串

题目大意：多个串的最长公共子串

后缀自动机模板题啦
对一个串建后缀自动机，其他串在上面匹配
每个状态记录每个串最长匹配长度的最小值
然后找到所有状态中答案最长的一个即可
需要注意的是这里不仅要用到自动机的性质 还要用到树的性质
如果一个状态被走到了 那么他的link也一定出现了 因为link对应的就是这个东东的一个子串 而我们跑自动机的时候不一定跑到了link这个点

说来说去就是循环里的第三行

for (int i=ncnt+1;i;i--){      int x=tmp[i];      minl[x]=min(minl[x],maxl[x]);      if (x && maxl[x]) maxl[st[x].link]=st[st[x].link].len;    }

放spoj上的代码吧

#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;inline char nc(){  static char buf[100000],*p1=buf,*p2=buf;  if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }  return *p1++;}inline int read(char *s){  char c=nc(); int len=0;  for (;!(c>='a' && c<='z');c=nc()) if (c==EOF) return 0;  for (;c>='a' && c<='z';s[++len]=c,c=nc()); s[++len]=0; return len-1;}const int N=250005;struct state{  int len,link,next[26];}st[N<<1];int ncnt,last;inline void Extend(char c){  int cur=++ncnt,p; c-='a';  st[cur].len=st[last].len+1;  for (p=last;p!=-1 && !st[p].next[c];p=st[p].link)    st[p].next[c]=cur;  if (p==-1)    st[cur].link=0;  else{    int q=st[p].next[c];    if (st[q].len==st[p].len+1)      st[cur].link=q;    else{      int nq=++ncnt;      st[nq].len=st[p].len+1;      st[nq].link=st[q].link;      for (int i=0;i<26;i++) st[nq].next[i]=st[q].next[i];      for (;p!=-1 && st[p].next[c]==q;p=st[p].link)    st[p].next[c]=nq;      st[q].link=st[cur].link=nq;    }  }  last=cur;}int tmp[N<<1],sum[N<<1];inline void Sort(){  for (int i=0;i<=ncnt;i++) sum[st[i].len]++;  for (int i=1;i<=st[last].len;i++) sum[i]+=sum[i-1];  for (int i=0;i<=ncnt;i++) tmp[sum[st[i].len]--]=i;}char A[N],B[N];int na,nb;int maxl[N<<1];int minl[N<<1];int main(){  freopen("t.in","r",stdin);  freopen("t.out","w",stdout);  na=read(A);  st[0].link=-1;  for (int i=1;i<=na;i++) Extend(A[i]);  Sort();  for (int i=0;i<=ncnt;i++) minl[i]=st[i].len;  while (nb=read(B)){    for (int i=0;i<=ncnt;i++) maxl[i]=0;    int p=0,l=0;    for (int i=1;i<=nb;i++){      while (p && !st[p].next[B[i]-'a'])    p=st[p].link,l=st[p].len;      if (st[p].next[B[i]-'a'])    p=st[p].next[B[i]-'a'],l++;      maxl[p]=max(maxl[p],l);    }    for (int i=ncnt+1;i;i--){      int x=tmp[i];      minl[x]=min(minl[x],maxl[x]);      if (x && maxl[x]) maxl[st[x].link]=st[st[x].link].len;    }  }  int Ans=0;  for (int i=0;i<=ncnt;i++)    Ans=max(Ans,minl[i]);  printf("%d\n",Ans);  return 0;}