BZOJ 2946: [Poi2000]公共串 后缀数组

2946: [Poi2000]公共串

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 1296  Solved: 576
[Submit][Status][Discuss]

Description

          给出几个由小写字母构成的单词，求它们最长的公共子串的长度。

任务：

l        读入单词

l        计算最长公共子串的长度

l        输出结果  

Input

  文件的第一行是整数 n，1<=n<=5，表示单词的数量。接下来n行每行一个单词，只由小写字母组成，单词的长度至少为1，最大为2000。  

Output

仅一行，一个整数，最长公共子串的长度。  

Sample Input

3
abcb
bca
acbc

Sample Output

２

二分答案+hash的做法很显然

不过BJ当然要选择逼格高的后缀数组

在单词中间加分隔符 求完ht之后

按rk从小到大扫一遍（具体看code吧）

之后就可以搞了

BJ这么搞 要特判只有一个串、所有串相同的情况

最开始就 WA 了

在 cog 关 O2 WA 一个点 没特判所有串相同

之后关 O2 A 了

一到有 O2 的bzoj 和 开 O2 cogs 就又 WA 了

疯狂认为是 O2 把我搞挂了 一顿狂调

最后发现后来的特判写的太不经意 没初始化

没开 O2 时因为数据范围小 随机的可能都比串长大 没影响

开了O2就挂了

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<set>
#include<map>
using namespace std;

typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=10100;

int n,tot;
char s[N];
int bel[N];

int sa1[N],sa2[N],rk1[N],rk2[N],buc[N],ht[N];
int *sa=sa1,*tp=sa2,*rk=rk1,*tmp=rk2;

void getsa()
{
	register int i,j,k,flag=0;
	for(i=1;i<=n;++i) buc[s[i]]++;
	for(i=1;i<=300;++i) buc[i]+=buc[i-1];
	for(i=n;i;i--) sa[buc[s[i]]--]=i;
	for(i=1;i<=n;++i) rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);
	
	for(k=1;k<=n;k<<=1)
	{
		for(i=1;i<=n;++i) buc[rk[sa[i]]]=i;
		for(i=n;i;i--)
			if(sa[i]>k)
				tp[buc[rk[sa[i]-k]]--]=sa[i]-k;
		for(i=n-k+1;i<=n;++i)tp[buc[rk[i]]--]=i;
		flag=0;
		for(i=1;i<=n;++i)
		{
			tmp[tp[i]]=tmp[tp[i-1]]+(rk[tp[i]]!=rk[tp[i-1]]||rk[tp[i]+k]!=rk[tp[i-1]+k]);
			if(tmp[tp[i]]==n)flag=1;
		}
		swap(rk,tmp);swap(sa,tp);
		if(flag)break;
	}
	
	k=0;
	for(i=1;i<n;++i)
	{
		j=rk[i]-1;
		while(s[sa[rk[i]]+k]==s[sa[j]+k])k++;
		ht[rk[i]]=k;
		if(k)k--;
	}
}

int mn[N][15],lg[N];

void initial()
{
	register int i,j;
	for(i=1;i<=n;++i)mn[i][0]=ht[i];
	for(j=1;(1<<j)<=n;++j)
		for(i=1;i+(1<<j)<=n+1;++i)
			mn[i][j]=min(mn[i][j-1],mn[i+(1<<(j-1))][j-1]);
	lg[1]=0;
	for(i=2;i<=n;++i)
		lg[i]= i>=(1<<(lg[i-1]+1)) ? lg[i-1]+1 : lg[i-1];
}

inline int query_mn(int l,int r)
{
	int k=lg[r-l+1];
	return min(mn[l][k],mn[r-(1<<k)+1][k]);
}

int ans;
int pos[6][N],sum[6][N];

void solve()
{
	register int i,j,tmp;
	register bool flag;
	for(i=1;i<=n;++i)
	{
		for(j=1;j<=tot;++j)
			sum[j][i]=sum[j][i-1];
		sum[bel[sa[i]]][i]++;
		pos[bel[sa[i]]][sum[bel[sa[i]]][i]]=i;
		flag=1;tmp=i;
		for(j=1;j<=tot;++j)
		{
			if(!sum[j][i]){flag=0;break;}
			tmp=min(tmp,pos[j][sum[j][i]]);
		}
		if(flag)ans=max(ans,query_mn(tmp+1,i));
	}
}

int main()
{
	scanf("%d",&tot);
	register int i,now,tmp=N;
	s[0]=126;
	for(i=1;i<=tot;++i)
	{
		scanf("%s",s+n+1);
		n=strlen(s+1);
		now=n;
		while(s[now]!=126)
			bel[now]=i,now--;
		tmp=min(tmp,n-now);
		s[n+1]=126;n++;
	}
	if(tot==1)
	{
		cout<<n-1<<endl;
		return 0;
	}
	getsa();
	initial();
	solve();
	cout<<min(tmp,ans)<<endl;
	return 0;
}
/*
3
abcb
bca
acbc

２
*/