记录编号 476085 评测结果 AAAAAAAAAAA
题目名称 [USACO Dec06]产奶的模式 最终得分 100
用户昵称 Gravatarcstdio 是否通过 通过
代码语言 C++ 运行时间 0.100 s
提交时间 2017-11-21 14:02:56 内存使用 18.09 MiB
显示代码纯文本
//之前就rank[sa[i]]=i
//没有判断sliding window里够不够K-1个
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>
using namespace std;
#define rank rank2
const int SIZEC=1000010,SIZEN=20010;
void basesort(int a[],int b[],int c[],int n,int m){
    if(m==-1){
        for(int i=0;i<n;i++) m=max(m,c[i]);
    }
    static int sum[SIZEC];
    for(int i=0;i<=m;i++) sum[i]=0;
    for(int i=0;i<n;i++) sum[c[a[i]]]++;
    for(int i=1;i<=m;i++) sum[i]+=sum[i-1];
    for(int i=n-1;i>=0;i--) b[--sum[c[a[i]]]]=a[i];
}
void sort_suf(int S[],int rank[],int sa[],int N){
    static int A[SIZEN],B[SIZEN],x[SIZEN],y[SIZEN];
    for(int i=0;i<N;i++) x[i]=S[i],A[i]=i;
    basesort(A,sa,x,N,-1);
    rank[sa[0]]=1;
    for(int i=1;i<N;i++){
        if(x[sa[i]]==x[sa[i-1]]) rank[sa[i]]=rank[sa[i-1]];
        else rank[sa[i]]=rank[sa[i-1]]+1;
    }
    for(int k=1;k<=N;k<<=1){
        for(int i=0;i<N;i++){
            x[i]=rank[i];
            y[i]=i+k<N?rank[i+k]:0;
            A[i]=i;
        }
        basesort(A,B,y,N,N);
        basesort(B,sa,x,N,N);
        rank[sa[0]]=1;
        for(int i=1;i<N;i++){
            if(x[sa[i]]==x[sa[i-1]]&&y[sa[i]]==y[sa[i-1]]) rank[sa[i]]=rank[sa[i-1]];
            else rank[sa[i]]=rank[sa[i-1]]+1;
        }
    }
    for(int i=0;i<N;i++) rank[sa[i]]=i;
}
void calc_height(int S[],int rank[],int sa[],int height[],int N){
    int h=0;
    for(int i=0;i<N;i++){
        if(rank[i]==0) h=0;
        else{
            int k=sa[rank[i]-1];
            if(--h<0) h=0;
            while(S[i+h]==S[k+h]) h++;
        }
        height[rank[i]]=h;
    }
}
int N,K;
int S[SIZEN];
int rank[SIZEN],sa[SIZEN];
int height[SIZEN];
void work(void){
    if(K==1){
        printf("%d\n",N);
        return;
    }
    sort_suf(S,rank,sa,N);
    calc_height(S,rank,sa,height,N);
    //for(int i=0;i<N;i++){for(int j=sa[i];j<N;j++){cout<<S[j]<<" ";}cout<<endl;}
    //for(int i=0;i<N;i++) cout<<height[i]<<" ";cout<<endl;
    multiset<int> st;
    int ans=0;
    for(int i=0;i<N;i++){
        st.insert(height[i]);
        if(i>=K-1){
            st.erase(height[i-K+1]);
        }
        if(st.size()==K-1) ans=max(ans,*st.begin());
    }
    printf("%d\n",ans);
}
void read(void){
    scanf("%d%d",&N,&K);
    for(int i=0;i<N;i++) scanf("%d",&S[i]);
}
int main(){
    //freopen("input.in","r",stdin);
    freopen("patterns.in","r",stdin);
    freopen("patterns.out","w",stdout);
    read();
    work();
    return 0;
}