CodeVS 3160 最长公共子串

最新推荐文章于 2020-12-19 09:20:56 发布

*ACoder*

最新推荐文章于 2020-12-19 09:20:56 发布

阅读量380

点赞数

CC 4.0 BY-SA版权

分类专栏： # 后缀数组

本文链接：https://blog.youkuaiyun.com/FSAHFGSADHSAKNDAS/article/details/52886949

后缀数组专栏收录该内容

21 篇文章

订阅专栏

本文介绍了一种使用后缀数组解决字符串匹配问题的算法实现。通过详细的代码示例，展示了如何构建后缀数组并利用其进行高效匹配。适用于对字符串处理和算法优化感兴趣的读者。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目

http://codevs.cn/problem/3160/
题解

刷水题辣

代码

//后缀数组水题 
#include <cstdio>
#include <algorithm>
#include <cstring>
#define maxn 200010
using namespace std;
int sa[maxn], rank[maxn], height[maxn], ws[maxn], wv[maxn], wa[maxn], wb[maxn], l1, l2,
	r[maxn];
char a[maxn], b[maxn];
bool cmp(int *y, int a, int b, int l)
{return y[a]==y[b] and y[a+l]==y[b+l];}
void init(int n, int m)
{
	int i, j, k=0, p, *x=wa, *y=wb, *t;
	for(i=0;i<m;i++)ws[i]=0;
	for(i=0;i<n;i++)ws[x[i]=r[i]]++;
	for(i=1;i<m;i++)ws[i]+=ws[i-1];
	for(i=n-1;i>=0;i--)sa[--ws[x[i]]]=i;
	for(j=p=1;p<n;m=p,j<<=1)
	{
		for(i=n-j,p=0;i<n;i++)y[p++]=i;
		for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
		for(i=0;i<n;i++)wv[i]=x[y[i]];
		for(i=0;i<m;i++)ws[i]=0;
		for(i=0;i<n;i++)ws[wv[i]]++;
		for(i=1;i<m;i++)ws[i]+=ws[i-1];
		for(i=n-1;i>=0;i--)sa[--ws[wv[i]]]=y[i];
		for(t=x,x=y,y=t,p=1,i=1,x[sa[0]]=0;i<n;i++)
			x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
	}
	for(i=0;i<n;i++)rank[sa[i]]=i;
	for(i=0;i<n-1;height[rank[i++]]=k)
		for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
}
void solve()
{
	int ans=0, i;
	for(i=2;i<=l1+l2;i++)
		if(sa[i-1]<l1 and sa[i]>l1 or sa[i-1]>l1 and sa[i]<l1)
			ans=max(ans,height[i]);
	printf("%d\n",ans);
}
int main()
{
	int i;
	scanf("%s",a);l1=strlen(a);
	scanf("%s",b);l2=strlen(b);
	for(i=0;i<l1;i++)r[i]=a[i];
	r[l1]='|';
	for(i=0;i<l2;i++)r[l1+1+i]=b[i];
	init(l1+1+l2,300);
	solve();
	return 0;
}