[后缀自动机 模板题] SPOJ 1811Longest Common Substring

来源:互联网 发布:中兴证券软件下载 编辑:程序博客网 时间:2024/05/17 14:15

题目大意:两个串的最长公共子串

后缀自动机模板题啦
对A建自动机 B在上面匹配
假设现在到B[i-1]的最长公共子串长度为l
如果当前节点有B[i]这个孩子,那么直接就l++即可。
如果没有就找一直向前找suffix_link,直到找到有B[i]这个孩子的节点
类似KMP中的next和AC自动机中的fail

#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());  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;}char A[N],B[N];int na,nb;int main(){  freopen("t.in","r",stdin);  freopen("t.out","w",stdout);  na=read(A);  nb=read(B);  st[0].link=-1;  for (int i=1;i<=na;i++)    extend(A[i]);  int p=0,l=0,maxv=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++;    maxv=max(maxv,l);  }  printf("%d\n",maxv);  return 0;}
0 0
原创粉丝点击