……另一道题用到，想看看部分代码对不对。于是又拿SAM交了一遍此题。
仅贴代码、
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 101000
#define T 30
using namespace std;
char s[N],t[N];
int pa[N<<1],son[N<<1][T];
int deep[N<<1],cnt,root,last;
inline int Newnode(int _deep){deep[++cnt]=_deep;return cnt;}
inline void SAM(int alp)
{
int u=last;
last=Newnode(deep[last]+1);
while(u&&!son[u][alp])son[u][alp]=last,u=pa[u];
if(!u)pa[last]=root;
else {
int v=son[u][alp];
if(deep[v]==deep[u]+1)pa[last]=v;
else {
int nv=Newnode(deep[u]+1);
memcpy(son[nv],son[v],sizeof son[v]);
pa[nv]=pa[v],pa[v]=pa[last]=nv;
while(u&&son[u][alp]==v)son[u][alp]=nv,u=pa[u];
}
}
}
int lena,lenb;
int solve(int n)
{
int i,x,v,ret=0,temp=0;
x=root=last=Newnode(0);
for(i=1;i<=lena;i++)SAM(s[i]-'a');
for(i=1;i<=lenb;i++)
{
if(son[x][t[i]-'a'])x=son[x][t[i]-'a'],temp++;
else {
while(x&&son[x][t[i]-'a']==0)x=pa[x];
if(!x)x=root,temp=0;
else temp=deep[x]+1,x=son[x][t[i]-'a'];
}
ret=max(ret,temp);
}
return ret;
}
int main()
{
// freopen("test.in","r",stdin);
int i,j,n,ans=0;
scanf("%s %s",s+1,t+1);
lena=strlen(s+1),lenb=strlen(t+1);
printf("%d\n",solve(n));
return 0;
}
【POJ2774】Long Long Message 后缀自动机
原文地址：http://blog.csdn.net/vmurder/article/details/42125337
