POJ 3261 Milk Patterns - 后缀数组（论文题）

最新推荐文章于 2019-10-06 08:14:52 发布

x_1023

最新推荐文章于 2019-10-06 08:14:52 发布

阅读量217

点赞数

CC 4.0 BY-SA版权

分类专栏：字符串后缀数组

本文链接：https://blog.youkuaiyun.com/x_1023/article/details/77659486

字符串同时被 2 个专栏收录

5 篇文章

订阅专栏

后缀数组

3 篇文章

订阅专栏

本文介绍了一种解决特定论文题目中关于寻找字符串中可重叠且至少重复k次的子序列最大长度的方法。通过二分搜索和后缀数组等技术实现高效查找。

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

大概又是一道论文题。问可重叠的至少重复k次的子序列的最大长度。

毅然决然二分长度，然后搜一遍就好了。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>

using namespace std;

const int maxn=20005;
const int maxm=100005;

int n,k;
int s[maxn];
int c[maxm];
int height[maxn],sa[maxn],rank[maxn];
int temp[maxm],tmp[maxm];

int read()
{
	int x=0;char ch=getchar();
	while(ch<'0'||ch>'9')ch=getchar();
	while(ch<='9'&&ch>='0')x=x*10+ch-'0',ch=getchar();
	return x;
}

void getsa()
{
	int *x=tmp,*y=temp,cnt,m=maxm-2;
	for(int i=1;i<=m;i++)c[i]=0;
	for(int i=1;i<=n;i++)c[x[i]=s[i]]++;
	for(int i=2;i<=m;i++)c[i]+=c[i-1];
	for(int i=n;i>=1;i--)sa[c[x[i]]--]=i;
	for(int k=1;k<=n;k<<=1)
	{
		cnt=0;
		for(int i=n-k+1;i<=n;i++)y[++cnt]=i;
		for(int i=1;i<=n;i++)if(sa[i]>k)y[++cnt]=sa[i]-k;
		for(int i=1;i<=m;i++)c[i]=0;
		for(int i=1;i<=n;i++)c[x[y[i]]]++;
		for(int i=2;i<=m;i++)c[i]+=c[i-1];
		for(int i=n;i>=1;i--)sa[c[x[y[i]]]--]=y[i];
		swap(x,y);
		x[sa[1]]=cnt=1;
		for(int i=2;i<=n;i++)
			x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?cnt:++cnt;
		if(cnt>=n)break;
		m=cnt;
	}
	for(int i=1;i<=n;i++)rank[sa[i]]=i;
	cnt=0;
	for(int i=1;i<=n;i++)
	{
		if(cnt)cnt--;
		int j=sa[rank[i]-1];
		while(s[i+cnt]==s[j+cnt])cnt++;
		height[rank[i]]=cnt;
	}
}
bool judge(int mid)
{
	int num=0;
	for(int i=1;i<=n;i++)
	{
		if(height[i]>=mid)
		{
			num++;
			if(num>=k-1)return true;
		}
		else num=0;
	}
	return false;
}
int main()
{
	scanf("%d%d",&n,&k);
	for(int i=1;i<=n;i++)s[i]=read()+1;
	getsa();
	int l=0,r=n+1,mid;
	while(l<r)
	{
		mid=l+r+1>>1;
		if(!judge(mid))r=mid-1;
		else l=mid;
	}
	printf("%d",l);
}