bzoj 2946(后缀自动机)

2946: [Poi2000]公共串

Time Limit: 3 Sec   Memory Limit: 128 MB
Submit: 681   Solved: 285
[ Submit][ Status][ Discuss]

Description

 
        给出几个由小写字母构成的单词,求它们最长的公共子串的长度。
任务:
l         读入单词
l         计算最长公共子串的长度
l         输出结果
 

Input

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

Output

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

Sample Input

3
abcb
bca
acbc

Sample Output


解题思路:后缀自动机。先对第一个单词建立后缀自动机,然后,后面的单词进行匹配,求得

每个单词对每个点的最长后缀的最小值(因为要公共),然后在求出最大的最小值。

加深对后缀自动机的理解。


#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int n,len,last,cnt,pa;
char c[21000];
int fa[21000],mx[21000],ch[21000][27];
int ans[5][21000];


inline int read()
{
char y; int x=0,f=1; y=getchar();
while (y<'0'||y>'9') {if (y=='-') f=-1; y=getchar();}
while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
return x*f;
}


void add(int s)
 {
  int p=last; int np=last=++cnt; mx[np]=mx[p]+1;
  while (p&&!ch[p][s]) ch[p][s]=np,p=fa[p];
  if (p==0) fa[np]=1;else
 {
  int q=ch[p][s];
  if (mx[q]==mx[p]+1)
  {
  fa[np]=q;
}else
 {
  int nq=++cnt; mx[nq]=mx[p]+1;
  fa[nq]=fa[q];
  fa[np]=fa[q]=nq;
  memcpy(ch[nq],ch[q],sizeof(ch[nq]));
  while (p!=0&&ch[p][s]==q)ch[p][s]=nq,p=fa[p];
 }
  } 
 }


int main()
{
n=read();
scanf("%s",c);
len=strlen(c);
last=cnt=1;
for (int i=1;i<=len;++i)
{
add(c[i-1]-'a');
}
int pa=0;
for (int i=1;i<=n-1;++i)
{
scanf("%s",c);
len=strlen(c);
int tmp=0,now=1;
for (int j=1;j<=len;++j)
{
  int s=c[j-1]-'a';
      if (ch[now][s])
       {
        ++tmp; now=ch[now][s]; ans[i][now]=max(ans[i][now],tmp); 
}else
{
while (now!=0 && ch[now][s]==0) now=fa[now];
if (now==0)
{
now=1; tmp=0;
}else
 {
  ans[i][now]=max(ans[i][now],mx[now]);
  tmp=mx[now]+1; now=ch[now][s]; ans[i][now]=max(ans[i][now],tmp); 
 }
}
}
}
for (int i=cnt;i>=2;--i)
{
for (int j=1;j<=n-1;++j)
 ans[j][fa[i]]=min(mx[fa[i]],max(ans[j][fa[i]],ans[j][i]));
}
for (int i=2;i<=cnt;++i)
{
int jo=0x7fffffff;
for (int j=1;j<=n-1;++j)
jo=min(jo,ans[j][i]);
pa=max(jo,pa);
}
printf("%d\n",pa);
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值