ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

【ybtoj】【kmp】公共子串

2021-09-17 15:01:51  阅读:143  来源: 互联网

标签:子串 前缀 ybtoj ll KMP long && kmp


题意

image

题解

本题的瓶颈在于不容易判断多个字符串中,两两匹配的公共子串是否是相同的。
模式串和匹配串的概念有点绕...在这里把要对自身求 \(nxt\) 的串记为 A,另一个串记为 B.
不要太着急,想到“定一移一”的方法,先单独取出任意一个字符串(本文取的是 \(s_1\))作为初始的 A 串,长度为 \(lenA\)。
那么每次枚举另外的 \(n-1\) 个串作为 B 串做 KMP,根据每次失配的位置 \(j\) 可以确定出 A 串的前缀 \([1,j-1]\) 是 B 串的一个子串,而且由于是前缀所以这部分子串是所有 B 串公共的。
那么所有 B 串中的(子串=A 串前缀)的最大长度,取最小值,就是一个答案。
再枚举 A 串从什么位置开始,每次新的 A 串为 \([i,lenA]\) 的部分,重新做 KMP。

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int INF = 0x3f3f3f3f,N = 2021;
inline ll read()
{
	ll ret=0;char ch=' ',c=getchar();
	while(!(c>='0'&&c<='9')) ch=c,c=getchar();
	while(c>='0'&&c<='9') ret=(ret<<1)+(ret<<3)+c-'0',c=getchar();
	return ch=='-'?-ret:ret;
}
char s[6][N],t[N];
int n,nxt[N],lt,lenT;
inline void pretreat()
{
	int j=0;
	nxt[1]=0;
	for(int i=1;i<lenT;i++) 
	{
		while(j&&t[i+1]!=t[j+1]) j=nxt[j];
		if(t[i+1]==t[j+1]) j++;
		nxt[i+1]=j;
	}
}
inline int kmp(int op)
{
	int lenn=strlen(s[op]+1),j=0,ret=0;
	for(int i=0;i<lenn;i++)
	{
		while(j&&s[op][i+1]!=t[j+1]) ret=max(ret,j-1),j=nxt[j];
		if(t[j+1]==s[op][i+1]) j++,ret=max(ret,j);
		if(j==lenT) return lenT;
	}
	return ret;
}
int main()
{
	n=read();
	for(int i=1;i<=n;i++) scanf("%s",s[i]+1);
	lt=strlen(s[1]+1);
	int ans=0;
	for(int i=1;i<=lt;i++)
	{
		lenT=lt-i+1;
		for(int j=1;j<=lenT;j++)
			t[j]=s[1][j+i-1];
		pretreat();
		int res=INF;
		for(int j=2;j<=n;j++) res=min(res,kmp(j));
		ans=max(ans,res);
	}
	printf("%d",ans);
	return 0;
}

标签:子串,前缀,ybtoj,ll,KMP,long,&&,kmp
来源: https://www.cnblogs.com/conprour/p/15304674.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有