POJ2774 Long Long Message【后缀数组(最长公共子串)】

题目描述：

RT，n<=100000，小写字母，求两个串的最长公共子串。

题目分析：

两个串拼起来，答案就是两个在不同串中的后缀的最长公共前缀。

我傻了。。写了后缀数组还去写二分分组。。
实际上只需要取height对应的两个后缀在不同串中的最大值就可以了。

Code：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 200005
using namespace std;
int b[maxn],ary[4][maxn],h[maxn];
int *sa=ary[0],*rk=ary[1],*nsa=ary[2],*nrk=ary[3];
void build_sa(int n,int m,int *a){
	memset(b,0,(m+1)<<2),a[n+1]=rk[n+1]=0;
	for(int i=1;i<=n;i++) b[a[i]]++;
	for(int i=1;i<=m;i++) b[i]+=b[i-1];
	for(int i=1;i<=n;i++) sa[b[a[i]]--]=i;
	for(int i=1;i<=n;i++) rk[sa[i]]=rk[sa[i-1]]+(a[sa[i-1]]!=a[sa[i]]);
	for(int k=1;rk[sa[n]]<n;k<<=1){
		for(int i=1;i<=n;i++) b[rk[sa[i]]]=i;
		for(int i=n;i>=1;i--) if(sa[i]>k) nsa[b[rk[sa[i]-k]]--]=sa[i]-k;
		for(int i=n-k+1;i<=n;i++) nsa[b[rk[i]]--]=i;
		for(int i=1;i<=n;i++) nrk[nsa[i]]=nrk[nsa[i-1]]+(rk[nsa[i-1]]!=rk[nsa[i]]||rk[nsa[i-1]+k]!=rk[nsa[i]+k]);
		swap(sa,nsa),swap(rk,nrk);	
	}
	for(int i=1,k=0,j;i<=n;h[rk[i]]=k,i++)
		for(k&&(k--),j=sa[rk[i]-1];a[i+k]==a[j+k];k++);
}
int n,m,a[maxn],ans;
char s0[maxn],s1[maxn];
int main()
{
	scanf("%s%s",s0+1,s1+1);
	n=strlen(s0+1),m=strlen(s1+1),a[n+1]=27;
	for(int i=1;i<=n;i++) a[i]=s0[i]-'a'+1;
	for(int i=1;i<=m;i++) a[n+1+i]=s1[i]-'a'+1;
	build_sa(n+m+1,27,a);
	for(int i=2;i<=n+m+1;i++)
		if((sa[i-1]<=n)==(sa[i]>n)) ans=max(ans,h[i]);
	printf("%d\n",ans);
}