记录编号 |
454407 |
评测结果 |
AAAAAAAAAAA |
题目名称 |
[USACO Dec06]产奶的模式 |
最终得分 |
100 |
用户昵称 |
Hzoi_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(){;}