字符串 后缀数组（可重叠的k 次最长重复子串）poj 3261

开始用起了后缀数组这种高级的东西，虽然细节方面是勉强看懂了，但是能运用好就行

要是有做这套题的朋友实在是不得其解的话，去看看这篇论文吧。。

IOI2009 国家集训队论文后缀数组罗穗骞

height数组当中存了从第一名到最后一名的公共前缀的长度

height 数组：定义height[i]=suffix(sa[i-1])和suffix(sa[i])的最长公共前缀，也就是排名相邻的两个后缀的最长公共前缀。那么对于j 和k，不妨设
rank[j]<rank[k],则有以下性质：
suffix(j) 和suffix(k) 的最长公共前缀为height[rank[j]+1],height[rank[j]+2], height[rank[j]+3], … ,height[rank[k]]中的最小值。

那么此题我们可以二分枚举其中的答案。

代码如下:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using std::sort;
int const M = 1001000;
int const N = 20100;
int rank[N],wb[M],wv[M],ws[M],height[N],num[N],sa[N];
int n,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=rank,*y=wb,*t;
    for(i=0; i<m; i++) ws[i]=0;
    for(i=0; i<n; i++) ws[x[i]=r[i]]++;
    for(i=1; i<m; i++) ws[i]+=ws[i-1];
    for(i=n-1; i>=0; i--) sa[--ws[x[i]]]=i;
    for(j=1,p=1; p<n; j*=2,m=p)
    {
        for(p=0,i=n-j; i<n; i++) y[p++]=i;
        for(i=0; i<n; i++) if(sa[i]>=j) y[p++]=sa[i]-j;
        for(i=0; i<n; i++) wv[i]=x[y[i]];
        for(i=0; i<m; i++) ws[i]=0;
        for(i=0; i<n; i++) ws[wv[i]]++;
        for(i=1; i<m; i++) ws[i]+=ws[i-1];
        for(i=n-1; i>=0; i--) sa[--ws[wv[i]]]=y[i];
        for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1; i<n; i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
    }
    return;
}

void calheight(int *r,int *sa,int n)
{
    int i,j,k=0;
    for(i=1; i<=n; i++) rank[sa[i]]=i;
    for(i=0; i<n; height[rank[i++]]=k)
        for(k?k--:0,j=sa[rank[i]-1]; r[i+k]==r[j+k]; k++);
    return;
}
bool check(int len)
{
     int tot=0;
     for(int i=2;i<=n;i++)
     {
         if(height[i]<len)tot=0;
         else
         {
            tot++;
            if(tot==k-1)return true;
         }
     }
     return false;
}
int main()
{
    while(~scanf("%d %d",&n,&k))
    {
          for(int i=0;i<n;i++)
          {
              scanf("%d",&num[i]);
          }
          da(num,sa,n+1,M);
          calheight(num,sa,n);
          int ans=0;
          int l=0,r=n;
          while(l<=r)
          {
                int m=(l+r)>>1;
                if(check(m))
                   l=m+1,ans=m;
                else
                   r=m-1;
          }
          printf("%d\n",ans);
    }
    return 0;
}