USACO06DEC 牛奶模式

题意：求最长的可重叠的 K重复子串 的长度

考虑二分长度s,转化为验证性问题。

对SA进行分组。保证组内Height最小为s。这样在组内RMQ就可以任意了，因为RMQ一定是大于S的。

只要组内元素个数大于等于K就是可行解。

#include <bits/stdc++.h>
using namespace std;

struct SA{
    int str[1000005];
    int x[1000005],y[1000005],u[1000005],v[1000005],r[1000005],o[1000005],hei[1000005],m=1000005,n;
    
    int a[5000005];
    
    void build(int p,int l,int r){
        if(l==r) a[p]=hei[l];
        else build(p*2,l,(l+r)/2), build(p*2+1,(l+r)/2+1,r), a[p]=min(a[p*2],a[p*2+1]);
    }
    
    int query(int p,int l,int r,int ql,int qr){
        if(l>qr||r<ql) return 1e+9;
        if(l>=ql&&r<=qr) return a[p];
        return min(query(p*2,l,(l+r)/2,ql,qr),query(p*2+1,(l+r)/2+1,r,ql,qr));
    }
    
    void calc(){
        int i,j,k=0;
        //for(i=1; i<=n; i++) r[x[i]]=i;
        for(i=1; i<=n; hei[r[i++]]=k)
            for(k?k--:0,j=x[r[i]-1]; str[i+k]==str[j+k]; k++);
    }
    
    int solve(){
        //n=strlen(str+1);
        //printf("len %d
",n);
        
        for(int i=1;i<=n;i++) u[str[i]]++;
        for(int i=1;i<=m;i++) u[i]+=u[i-1];
        for(int i=n;i>=1;i--) x[u[str[i]]--]=i;
        r[x[1]]=1;
        for(int i=2;i<=n;i++) r[x[i]]=r[x[i-1]]+((str[x[i]]-str[x[i-1]])?1:0);
        
        for(int l=1;r[x[n]]<n;l<<=1) {
            memset(u,0,sizeof u);
            memset(v,0,sizeof v);
            memcpy(o,r,sizeof r);
            for(int i=1;i<=n;i++) u[r[i]]++, v[(i+l<=n)?r[i+l]:0]++;
            for(int i=1;i<=n;i++) u[i]+=u[i-1], v[i]+=v[i-1];
            for(int i=n;i>=1;i--) y[v[(i+l<=n)?r[i+l]:0]--]=i;
            for(int i=n;i>=1;i--) x[u[r[y[i]]]--]=y[i];
            r[x[1]]=1;
            for(int i=2;i<=n;i++) 
                r[x[i]]=r[x[i-1]]+((o[x[i]] != o[x[i-1]])  ||  (((x[i]+l<=n)?o[x[i]+l]:0) != ((x[i-1]+l<=n)?o[x[i-1]+l]:0)));
        }
        calc();
        hei[1]=0;
        //build(1,1,n);
    }
    
    int lcp(int pos1,int pos2){
        return query(1,1,n,min(r[pos1],r[pos2])+1,max(r[pos1],r[pos2]));
    }
    
} sa;

int n,k;

int main(){
    cin>>n>>k;    
    sa.n=n;
    for(int i=1;i<=n;i++) cin>>sa.str[i];
    sa.solve();
    int l=1,r=1000000,ans=0;
    while(r-l) {
        int s=(l+r)/2,cnt=1,mx=0;
        for(int i=1;i<=n+1;i++)
            if(sa.hei[i]<s)
                mx=max(cnt,mx), cnt=1;
            else cnt++;
        if(mx>=k) ans=max(ans,s), l=s+1;
        else r=s;
    }
    cout<<ans<<endl;
            
        
    
}