poj 3261 二分答案+后缀数组 求至少出现k次的最长重复子序列

 #include "stdio.h"
 #define maxn 20010

 int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
 int rank[maxn],height[maxn];
 int r[maxn],sa[maxn],ans[maxn];
 int n,res,k;

 int cmp(int *r,int a,int b,int l)
 {
     return r[a]==r[b]&&r[a+l]==r[b+l];
 }

 void da(int *r,int *sa,int n,int m)
 {
     int i,j,p,*x=wa,*y=wb,*t;
     for(i=; i<m; i++) ws[i]=;
     for(i=; i<n; i++) ws[x[i]=r[i]]++;
     for(i=; i<m; i++) ws[i]+=ws[i-];
     for(i=n-; i>=; i--) sa[--ws[x[i]]]=i;
     for(j=,p=; p<n; j*=,m=p)
     {

         for(p=,i=n-j; i<n; i++) y[p++]=i;
         for(i=; i<n; i++) if(sa[i]>=j) y[p++]=sa[i]-j;
         for(i=; i<n; i++) wv[i]=x[y[i]];
         for(i=; i<m; i++) ws[i]=;
         for(i=; i<n; i++) ws[wv[i]]++;
         for(i=; i<m; i++) ws[i]+=ws[i-];
         for(i=n-; i>=; i--) sa[--ws[wv[i]]]=y[i];
         for(t=x,x=y,y=t,p=,x[sa[]]=,i=; i<n; i++)
             x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
     }
     return;
 }

 void calheight(int *r,int *sa,int n)
 {
     int i,j,k=;
     for(i=; i<=n; i++) rank[sa[i]]=i;
     for(i=; i<n; height[rank[i++]]=k)
         for(k?k--:,j=sa[rank[i]-]; r[i+k]==r[j+k]; k++);
     return;
 }

 bool calc(int x)        //x:length
 {
     int rec[maxn];
     int nm=;
     for (int i=;i<=n;i++)
     {
         if (height[i]<x)
         {
             nm++;
             rec[nm]=i;
         }
     }
     nm++;   rec[nm]=n+;
     for (int i=;i<=nm;i++)
     {
         int tl=rec[i-],tr=rec[i];
         if (tr-tl>=k)      return true;
     }
     return false;
 }

 int main()
 {
     scanf("%d %d",&n,&k);
     for (int i=; i<n; i++)
         scanf("%d",&r[i]);
     r[n]=;

     da(r,sa,n+,);
     calheight(r,sa,n);

     //TODO: Please Code Here
     int l=,r=n,res=;
     while (r>=l)
     {
         int mid=(l+r)/;
         if (calc(mid))
         {
             if (mid>res)    res=mid;
             l=mid+;
         }
         else
             r=mid-;
     }
     printf("%d\n",res);
     return ;
 }