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