[CODEVS3160]最长公共子串|后缀数组|后缀自动机

  第一种做法是把两个字符串接起来，中间放一个奇怪的字符，然后建这个串的后缀数组，求出h数组，对于h[i]表示lcp（sa[i],sa[i-1]），如果sa[i]和sa[i-1]分布在奇怪的字符的两边就用h[i]更新答案。。

  然后补了一下后缀自动机，就是裸题了，对一个串建SAM，另一个串在上面匹配即可。。

SA:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<memory.h>
#define N 400005
#define clr(a) memset(a,0,sizeof(a))
using namespace std;
char s
,b
;
int n,i,d,sa
,rank
,h
,buc
,sec
,rk
;
bool cmp(int a,int b)
{
return rk[a]==rk[b]&&rk[a+d]==rk[b+d];
}
void getsa()
{
int i,j,p=1;
d=1;clr(buc);clr(rank);
for (i=1;i<=n;i++) buc[rank[i]=s[i]-'a'+1]++;
for (i=1;i<=256;i++) buc[i]+=buc[i-1];
for (i=n;i;i--) sa[buc[rank[i]]--]=i;
while (p<n)
{
p=0;clr(buc);
for (i=n-d+1;i<=n;i++) sec[++p]=i;
for (i=1;i<=n;i++) if (sa[i]>d) sec[++p]=sa[i]-d;
for (i=1;i<=n;i++) buc[rank[sec[i]]]++,rk[i]=rank[i];
for (i=1;i<=n;i++) buc[i]+=buc[i-1];
for (i=n;i;i--) sa[buc[rank[sec[i]]]--]=sec[i];
for (i=2,p=1,rank[sa[1]]=1;i<=n;i++)
rank[sa[i]]=cmp(sa[i-1],sa[i])? p:++p;
d<<=1;
}
}
void geth()
{
int i,j,k=0;
h[1]=0;
for (i=1;i<=n;i++)
if (rank[i]>1)
{
j=sa[rank[i]-1];
while (s[i+k]==s[j+k]) k++;
h[rank[i]]=k;if (k) k--;
}
}
int main()
{
scanf("%s%s",s+1,b+1);
n=strlen(s+1);
int n1=n+1;
s[n+1]='z'+1;
for (i=n+2;i<=n+1+strlen(b+1);i++) s[i]=b[i-n-1];
n+=strlen(b+1)+1;
getsa();geth();
int ans=0;
for (i=2;i<=n;i++)
if ((sa[i]<n1&&sa[i-1]>n1)||(sa[i]>n1&&sa[i-1]<n1)) ans=max(ans,h[i]);
printf("%d\n",ans);
}

SAM:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<memory.h>
#define N 200005
using namespace std;
char s1
,s2
;
int n,ans,i,now,nd,a
[26],fa
,l
;
void extend(int c)
{
int p=now,np=now=++nd;
l[np]=l[p]+1;
while (!a[p][c]&&p) a[p][c]=np,p=fa[p];
if (!p) fa[np]=1;
else
{
int q=a[p][c];
if (l[q]==l[p]+1) fa[np]=q;
else
{
int nq=++nd;l[nq]=l[p]+1;
fa[nq]=fa[q];
fa[q]=fa[np]=nq;
for (int i=0;i<26;i++) a[nq][i]=a[q][i];
while (p&&a[p][c]==q) a[p][c]=nq,p=fa[p];
}
}
}
void build()
{
for (int i=1;i<=strlen(s1+1);i++) extend(s1[i]-'a');
}
void solve()
{
ans=0;now=1;
int tmp=0;
for (i=1;i<=strlen(s2+1);i++)
{
if (a[now][s2[i]-'a']) tmp++,now=a[now][s2[i]-'a'];
else
{
while (!a[now][s2[i]-'a']) now=fa[now];
if (!now) tmp=0,now=1;else tmp=l[now]+1,now=a[now][s2[i]-'a'];
}
ans=max(ans,tmp);
}
}
int main()
{
scanf("%s",s1+1);
memset(a,0,sizeof(a));
l[1]=0;now=nd=1;fa[1]=0;
for (i=0;i<26;i++) a[0][i]=1;
build();
scanf("%s",s2+1);
solve();
printf("%d\n",ans);
}