BZOJ 2946: [Poi2000]公共串【后缀自动机/后缀数组】

题面：

给出几个由小写字母构成的单词，求它们最长的公共子串的长度。
文件的第一行是整数 n，1<=n<=5，表示单词的数量。接下来n行每行一个单词，只由小写字母组成，单词的长度至少为1，最大为2000。

题目分析：

后缀数组做法：
将 n 个字符串连起来，中间用不相同的且没有出现在字符串中的字符隔开，
求后缀数组。然后二分答案，按照height值将后缀分成若干组，判断每组
的后缀是否包含所有的原串。

后缀自动机做法：摘自yyb的博客

后缀自动机的Code：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 4005
#define maxc 26
using namespace std;
int n,m,f[maxn],ans[maxn],b[maxn],id[maxn];
char s[maxn];
int fail[maxn]={-1},ch[maxn][maxc],len[maxn],sz,last;
void sa_extend(int c){
	int p=last,cur=++sz,q;
	len[last=cur]=len[p]+1;
	for(;p!=-1&&!ch[p][c];p=fail[p]) ch[p][c]=cur;
	if(p==-1) fail[cur]=0;
	else if(len[p]+1==len[q=ch[p][c]]) fail[cur]=q;
	else{
		int clone=++sz; len[clone]=len[p]+1,fail[clone]=fail[q];
		fail[q]=fail[cur]=clone;
		memcpy(ch[clone],ch[q],sizeof ch[q]);
		for(;p!=-1&&ch[p][c]==q;p=fail[p]) ch[p][c]=clone;
	}
}
int main()
{
	scanf("%d%s",&n,s+1),m=strlen(s+1);
	for(int i=1;i<=m;i++) sa_extend(s[i]-'a');
	for(int i=1;i<=sz;i++) b[ans[i]=len[i]]++;
	for(int i=1;i<=m;i++) b[i]+=b[i-1];
	for(int i=1;i<=sz;i++) id[b[len[i]]--]=i;
	while(--n)
	{
		scanf("%s",s+1),m=strlen(s+1);
		int p=0,l=0,v;
		memset(f,0,sizeof f);
		for(int i=1;i<=m;i++)
		{
			v=s[i]-'a';
			while(p&&!ch[p][v]) p=fail[p],l=len[p];
			if(ch[p][v]) p=ch[p][v],l++;
			f[p]=max(f[p],l);
		}
		for(int i=sz;i>=1;i--) f[fail[id[i]]]=max(f[fail[id[i]]],f[id[i]]);
		for(int i=1;i<=sz;i++) ans[i]=min(ans[i],f[i]);
	}
	printf("%d",*max_element(ans+1,ans+1+sz));
}

好像BZOJ上不更新fail也能过，但是很容易就能hack掉。。比如3 bbccbbcc bc cb