codevs 3160(后缀自动机复习)

题目描述 Description

给出两个由小写字母组成的字符串，求它们的最长公共子串的长度。

输入描述 Input Description

读入两个字符串

输出描述 Output Description

输出最长公共子串的长度

样例输入 Sample Input

C++

1 2	yeshowmuchiloveyoumydearmotherreallyicannotbelieveit yeaphowmuchiloveyoumydearmother

样例输出 Sample Output

27

数据范围及提示 Data Size & Hint

单个字符串的长度不超过100000

解题思路：后缀自动机

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
char a[100001],b[100001];
int last,cnt,ans,len1,len2,now;
int mx[100010],fa[100010];
int ch[100010][27];


inline int read()
 {
  char y; int x=0,f=1; y=getchar();
  while (y<'0' || y>'9') {if (y=='-') f=-1; y=getchar();}
while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
return x*f;
 }


void build(int s)
 {
  int p=last; int np=last=++cnt; mx[np]=mx[p]+1;
  while (p!=0 && !ch[p][s]) ch[p][s]=np,p=fa[p];
  if (p==0) fa[np]=1;else
  {
  int q=ch[p][s];
  if (mx[q]==mx[p]+1)
  {
  fa[np]=q;
}else
 {
  int nq=++cnt; mx[nq]=mx[p]+1;
  memcpy(ch[nq],ch[q],sizeof(ch[nq]));
  fa[nq]=fa[q];
  fa[q]=fa[np]=nq;
  while (p!=0 && ch[p][s]==q) ch[p][s]=nq,p=fa[p];
 }
 }
 }


int main()
{
scanf("%s",a); scanf("%s",b);
len1=strlen(a); len2=strlen(b);
cnt=1; last=1;
for (int i=1;i<=len1;++i)
{
build(a[i-1]-'a');
}
now=1;
int ans=-1; int tmp=0; 
for (int i=1;i<=len2;++i)
{
  int s=b[i-1]-'a';
  if (ch[now][s]) ++tmp,now=ch[now][s];else
   {
    while (!ch[now][s] && now!=0) now=fa[now];
    if (now==0) now=1,tmp=0;else
    tmp=mx[now]+1,now=ch[now][s];
}
ans=max(ans,tmp);
}
printf("%d",ans);
}