poj2774 Long Long Message（后缀数组）

很简单的后缀数组，求两个串的最长公共子串。

把这两个串，用特殊字符隔开连在一起，更新合法的最大公共前缀即可。

#include <cstdio>
#include <cstring>
#define N 200005
char a[N];
int rank[N<<1],sa[N],st[256],rank1[N],count[N],tmp[N],ans=0,h[N];
inline int max(int x,int y){return x>y?x:y;}
int main(){
//	freopen("a.in","r",stdin);
	scanf("%s",a+1);
	int len=strlen(a+1);
	a[++len]='z'+1;
	scanf("%s",a+len+1);
	int n=strlen(a+1);
	memset(st,0,sizeof(st));
	memset(rank,0,sizeof(rank));
	memset(rank1,0,sizeof(rank1));
	for(int i=1;i<=n;++i) st[a[i]]=1;
	for(int i=1;i<=255;++i) st[i]+=st[i-1];
	for(int i=n;i>=1;--i) rank[i]=st[a[i]];
	int k=0;
	for(int l=1;k!=n;l<<=1){
		memset(count,0,sizeof(count));
		for(int i=1;i<=n;++i) count[rank[i+l]]++;
		for(int i=1;i<=n;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) tmp[count[rank[i+l]]--]=i;
		memset(count,0,sizeof(count));
		for(int i=1;i<=n;++i) count[rank[tmp[i]]]++;
		for(int i=1;i<=n;++i) count[i]+=count[i-1];
		for(int i=n;i>=1;--i) sa[count[rank[tmp[i]]]--]=tmp[i];
		memcpy(rank1,rank,sizeof(rank)>>1);
		rank[sa[1]]=k=1;
		for(int i=2;i<=n;++i){
			if(rank1[sa[i]]!=rank1[sa[i-1]]||rank1[sa[i]+l]!=rank1[sa[i-1]+l]) ++k;
			rank[sa[i]]=k;
		}
	}k=0;
	for(int i=1;i<=n;++i){//i--第i个后缀 
		if(rank[i]==1){h[1]=0;continue;}
		if(i==1||h[rank[i-1]]<=1) k=0;
		if(k) --k;
		while(a[i+k]==a[sa[rank[i]-1]+k]) ++k;
		h[rank[i]]=k;
	}
	for(int i=2;i<=n;++i){//i--排名第i 
		if(!(sa[i]>len^sa[i-1]<=len)) ans=max(ans,h[i]);
	}
	printf("%d\n",ans);
	return 0;
}