记录编号 454407 评测结果 AAAAAAAAAAA
题目名称 [USACO Dec06]产奶的模式 最终得分 100
用户昵称 GravatarHzoi_Mafia 是否通过 通过
代码语言 C++ 运行时间 0.015 s
提交时间 2017-09-28 20:56:22 内存使用 2.61 MiB
显示代码纯文本
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
inline int read(){
	int sum(0);
	char ch(getchar());
	for(;ch<'0'||ch>'9';ch=getchar());
	for(;ch>='0'&&ch<='9';sum=sum*10+(ch^48),ch=getchar());
	return sum;
}
int n,m,k,size;
int num[20005],s[20005];
int t1[20005],t2[20005],t3[20005],buc[20005];
int sa[20005],rank[20005],height[20005];
inline void Suffix(){
	int i,j,k(0),p(0),*x(t1),*y(t2),*t;
	for(i=0;i<=m;++i)buc[i]=0;
	for(i=1;i<=n;++i)++buc[x[i]=s[i]];
	for(i=1;i<=m;++i)buc[i]+=buc[i-1];
	for(i=n;i>=1;--i)sa[buc[x[i]]--]=i;
	for(j=1;p<n;j<<=1,m=p){
		for(p=0,i=n-j+1;i<=n;++i)y[++p]=i;
		for(i=1;i<=n;++i)
			if(sa[i]>j)
				y[++p]=sa[i]-j;
		for(i=0;i<=m;++i)buc[i]=0;
		for(i=1;i<=n;++i)t3[i]=x[y[i]];
		for(i=1;i<=n;++i)++buc[t3[i]];
		for(i=1;i<=m;++i)buc[i]+=buc[i-1];
		for(i=n;i>=1;--i)sa[buc[t3[i]]--]=y[i];
		for(t=x,x=y,y=t,x[sa[1]]=1,p=1,i=2;i<=n;++i)
			x[sa[i]]=((y[sa[i]]==y[sa[i-1]])&&(y[sa[i]+j]==y[sa[i-1]+j]))?p:++p;
	}
	for(i=1;i<=n;++i)rank[sa[i]]=i;
	for(i=1;i<=n;height[rank[i++]]=k)
		for(k?--k:0,j=sa[rank[i]-1];s[i+k]==s[j+k];++k);
}
inline bool check(int x){
	int cnt(0),ans(0);
	for(int i=1;i<=n;++i){
		if(height[i]<x){
			ans=max(ans,cnt);
			cnt=0;
		}
		++cnt;
	}
	return max(ans,cnt)>=k;
}
inline int gg(){
	freopen("patterns.in","r",stdin);
	freopen("patterns.out","w",stdout);
	m=n=read(),k=read();
	for(int i=1;i<=n;++i)
		s[i]=num[i]=read();
	sort(num+1,num+n+1);
	size=unique(num+1,num+n+1)-num-1;
	for(int i=1;i<=n;++i)
		s[i]=lower_bound(num+1,num+size+1,s[i])-num;
	Suffix();
	int l(0),r(n),mid,ans(0);
	while(l<=r){
		mid=(l+r)>>1;
		if(check(mid))
			ans=mid,l=mid+1;
		else
			r=mid-1;
	}
	printf("%d",ans);
	return 0;
}
int K(gg());
int main(){;}