[后缀自动机 模板题] SPOJ 1812 Longest 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;
}