2946: [Poi2000]公共串
Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 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
abcb
bca
acbc
Sample Output
2
二分答案+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
2
*/