POJ 3261 后缀数组

题意：

最长的出现k次的重复（可重叠）子串

 

题解：

和POJ 1743一样，稍作改变，求最大的话一般都是对后缀数组二分处理。。

http://www.cnblogs.com/proverbs/archive/2013/02/05/2893419.html

判定条件从位置相差>=m改为是否存在一个块内后缀的个数大于等于m

 

数比较大，快排搞定，好像基数排序也可以。

 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <algorithm>

#define N 50500

using namespace std;

int r[N],wa[N],wc[N],wv[N],sa[N],wb[N],height[N],rank[N];
int n,m;

inline bool cmp1(int a,int b)
{
    if(r[a]==r[b]) return a<b;
    return r[a]<r[b];
}

inline bool cmp2(int a,int b)
{
    if(wv[a]==wv[b]) return a<b;
    return wv[a]<wv[b];
}

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

inline void da(int *r,int *sa,int n,int m)
{
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0;i<n;i++) wc[i]=i,x[i]=r[i];
    sort(wc,wc+n,cmp1);
    for(i=0;i<n;i++) sa[i]=wc[i];
    for(j=1,p=1;p<n;j<<=1,m=p)
    {
        for(i=n-j,p=0;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<n;i++) wc[i]=i;
        sort(wc,wc+n,cmp2);
        for(i=0;i<n;i++) sa[i]=y[wc[i]];
        for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=0;i<n;i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
    }
}

inline void getheight(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++);
}

inline bool check(int x)
{
    int cnt=1;
    for(int i=1;i<=n;i++)
    {
        if(height[i]<x) cnt=1;
        else cnt++;
        if(cnt>=m) return true;
    }
    return false;
}

inline void getans()
{
    int l=1,r=n,mid,res;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(check(mid)) res=mid,l=mid+1;
        else r=mid-1;
    }
    printf("%d\n",res);
}

inline void go()
{
    for(int i=0;i<n;i++) scanf("%d",&r[i]);
    r[n]=0;
    da(r,sa,n+1,1000010);
    getheight(r,sa,n);
    getans();
}

int main()
{
    while(scanf("%d%d",&n,&m)!=EOF) go();
    return 0;
}