题目描述:给出两个字符串,求最长公共子串的长度。
题解:后缀数组。将两个字符串合并成一个,中间加一个字符'0',然后求出height数组。因为以最长公共子串为前缀的子串,他们的rank肯定相差为一,然后遍历sa数组,判断sa[i]和sa[i-1]的位置是否位于字符'0'的两侧,符合条件的最大的height值即为答案。
代码:
#include<bits/stdc++.h>
#define N 200010
using namespace std;
char s[N];
int sa[N],t[N],t2[N],c[N],n,rak[N],height[N];
void build_sa(int m,char *s)
{
int i,*x=t,*y=t2;
for (i=0;i<m;i++)c[i]=0;
for (i=0;i<n;i++)c[x[i]=s[i]]++;
for (i=1;i<m;i++)c[i]+=c[i-1];
for (i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
for (int k=1;k<=n;k<<=1)
{
int p=0;
for (i=n-k;i<n;i++)y[p++]=i;
for (i=0;i<n;i++)if (sa[i]>=k)y[p++]=sa[i]-k;
for (i=0;i<m;i++)c[i]=0;
for (i=0;i<n;i++)c[x[y[i]]]++;
for (i=0;i<m;i++)c[i]+=c[i-1];
for (i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
swap(x,y);
p=1; x[sa[0]]=0;
for (i=1;i<n;i++)
x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
if (p>=n) break;
m=p;
}
}
void getheight()
{
int i,j,k=0;
for (i=0;i<n;i++)rak[sa[i]]=i;
for (i=0;i<n;i++)
{
if (k)k--;
if (!rak[i])continue;
j=sa[rak[i]-1];
while (s[i+k]==s[j+k])k++;
height[rak[i]]=k;
}
}
int main()
{
while (~scanf("%s",s))
{
n=strlen(s);
s[n]='0'; int nn=n,ans=0;
scanf("%s",s+n+1);
n=strlen(s);
build_sa(200,s);
getheight();
for (int i=1;i<n;i++)
{
int x=sa[i-1],y=sa[i];
if (x>y) swap(x,y);
if (x<nn && nn<y) ans=max(ans,height[i]);
}
printf("%d\n",ans);
}
return 0;
}