spoj 1811 LCS (后缀自动机 SAM)

http://www.spoj.com/problems/LCS/

题意：求两个字符串A,B的最长公共子串。字符串长度不超过250000。

 

思路：这应该算是后缀自动机的经典应用了吧，我们先构造A的SAM，然后用A的SAM一次读入B的每一个字符，初始时状态在root处，此时最大匹配数为tmp=0，（这里的最大匹配数是指以当前读入的字符结尾，往前能匹配的最大长度），设当前到达的状态为p，最大匹配数为tmp，读入的字符为x，若p->go[x]!=NULL，则说明可从当前状态读入一个字符x到达下一个状态，则tmp++,p=p->go[x],否则，找到p的第一个祖先s，s->go[x]!=NULL,若s不存在，则说明以x结尾的字符串无法和A串的任何位置匹配，则设tmp=0,p=root。否则，设tmp=s->tmp+1（因为我们不算x的话已经到达了状态p，这说明对于p的任意祖先已经匹配完毕）,p=s->go[x]。我们求tmp所达到的最大值即为所求。代码如下：

#include <iostream>
#include <string.h>
#include <algorithm>
#include <stdio.h>
#define maxn 500010
using namespace std;
char str[maxn>>1];
struct node
{
node *par,*go[26];
int val;
}*root,*tail,que[maxn],*top[maxn];
int tot,len,c[maxn];
void add(int c,int l)
{
node *np=&que[tot++],*p=tail;
np->val=l;
while(p&&p->go[c]==NULL)
{
p->go[c]=np;
p=p->par;
}
if(p==NULL)
np->par=root;
else
{
node *q=p->go[c];
if(q->val==p->val+1)
np->par=q;
else
{
node *nq=&que[tot++];
*nq=*q;
nq->val=p->val+1;
np->par=q->par=nq;
while(p&&p->go[c]==q)
{
p->go[c]=nq;
p=p->par;
}
}
}
tail=np;
}
void init(int n)
{
for(int i=0;i<=n;i++)
{
memset(que[i].go,0,sizeof(que[i].go));
que[i].val=0;
}
tot=0;
len=1;
root=tail=&que[tot++];
}
void solve()
{
int i;
scanf("%s",str);
int l=strlen(str),ans=0,tmp=0;
node *p=root;
for(i=0;i<l;i++)
{
int x=str[i]-'a';
if(p->go[x])
{
p=p->go[x];
tmp++;
}
else
{
while(p&&p->go[x]==NULL)
p=p->par;
if(p)
{
tmp=p->val+1;
p=p->go[x];
}
else
{
p=root;
tmp=0;
}
}
if(tmp>ans)
ans=tmp;
}
printf("%d\n",ans);
}
int main()
{
//freopen("dd.txt","r",stdin);
scanf("%s",str);
int l=strlen(str);
init(l*2);
for(int i=0;i<l;i++)
add(str[i]-'a',len++);
solve();
return 0;
}