USACO06DEC 牛奶模式

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

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

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

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

 #include <bits/stdc++.h>
 using namespace std;
 
 struct SA{
     int str[];
     int x[],y[],u[],v[],r[],o[],hei[],m=,n;
 
     int a[];
 
     void build(int p,int l,int r){
         if(l==r) a[p]=hei[l];
         else build(p*,l,(l+r)/), build(p*+,(l+r)/+,r), a[p]=min(a[p*],a[p*+]);
     }
 
     int query(int p,int l,int r,int ql,int qr){
         if(l>qr||r<ql) return 1e+;
         if(l>=ql&&r<=qr) return a[p];
         return min(query(p*,l,(l+r)/,ql,qr),query(p*+,(l+r)/+,r,ql,qr));
     }
 
     void calc(){
         int i,j,k=;
         //for(i=1; i<=n; i++) r[x[i]]=i;
         for(i=; i<=n; hei[r[i++]]=k)
             for(k?k--:,j=x[r[i]-]; str[i+k]==str[j+k]; k++);
     }
 
     int solve(){
         //n=strlen(str+1);
         //printf("len %d\n",n);
 
         for(int i=;i<=n;i++) u[str[i]]++;
         for(int i=;i<=m;i++) u[i]+=u[i-];
         for(int i=n;i>=;i--) x[u[str[i]]--]=i;
         r[x[]]=;
         for(int i=;i<=n;i++) r[x[i]]=r[x[i-]]+((str[x[i]]-str[x[i-]])?:);
 
         for(int l=;r[x[n]]<n;l<<=) {
             memset(u,,sizeof u);
             memset(v,,sizeof v);
             memcpy(o,r,sizeof r);
             for(int i=;i<=n;i++) u[r[i]]++, v[(i+l<=n)?r[i+l]:]++;
             for(int i=;i<=n;i++) u[i]+=u[i-], v[i]+=v[i-];
             for(int i=n;i>=;i--) y[v[(i+l<=n)?r[i+l]:]--]=i;
             for(int i=n;i>=;i--) x[u[r[y[i]]]--]=y[i];
             r[x[]]=;
             for(int i=;i<=n;i++)
                 r[x[i]]=r[x[i-]]+((o[x[i]] != o[x[i-]])  ||  (((x[i]+l<=n)?o[x[i]+l]:) != ((x[i-]+l<=n)?o[x[i-]+l]:)));
         }
         calc();
         hei[]=;
         //build(1,1,n);
     }
 
     int lcp(int pos1,int pos2){
         return query(,,n,min(r[pos1],r[pos2])+,max(r[pos1],r[pos2]));
     }
 
 } sa;
 
 int n,k;
 
 int main(){
     cin>>n>>k;
     sa.n=n;
     for(int i=;i<=n;i++) cin>>sa.str[i];
     sa.solve();
     int l=,r=,ans=;
     while(r-l) {
         int s=(l+r)/,cnt=,mx=;
         for(int i=;i<=n+;i++)
             if(sa.hei[i]<s)
                 mx=max(cnt,mx), cnt=;
             else cnt++;
         if(mx>=k) ans=max(ans,s), l=s+;
         else r=s;
     }
     cout<<ans<<endl;
 
 }