[poj 3261]后缀数组+滑窗最小值

题目链接：http://poj.org/problem?id=3261

这个是可以交叉的重复串，所以用height就可以了，但是题目说让重复k次以上，也就是直接做一个k-1长度的滑窗最小值，从这些最小值里取最大即可。

这里其实为了节省空间可以先给数字离散化一下，这样就只有20000了，不过不离散化空间也够用。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;

const int maxn=;

#define F(x) ((x)/3+((x)%3==1?0:tb))
#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)
int wa[maxn*],wb[maxn*],wv[maxn*],wss[maxn*];
int c0(int *r,int a,int b)
{
    return r[a]==r[b] && r[a+]==r[b+] && r[a+]==r[b+];
}
int c12(int k,int *r,int a,int b)
{
    if (k==) return r[a]<r[b] || (r[a]==r[b]&&c12(,r,a+,b+));
    else return r[a]<r[b] || (r[a]==r[b]&&wv[a+]<wv[b+]);
}
void sort(int *r,int *a,int *b,int n,int m)
{
    int i;
    for (i=;i<n;i++) wv[i]=r[a[i]];
    for (i=;i<m;i++) wss[i]=;
    for (i=;i<n;i++) wss[wv[i]]++;
    for (i=;i<m;i++) wss[i]+=wss[i-];
    for (i=n-;i>=;i--) b[--wss[wv[i]]]=a[i];
}
void dc3(int *r,int *sa,int n,int m)
{
    int i,j,*rn=r+n;
    int *san=sa+n,ta=,tb=(n+)/,tbc=,p;
    r[n]=r[n+]=;
    for (i=;i<n;i++) if (i%!=) wa[tbc++]=i;
    sort(r+,wa,wb,tbc,m);
    sort(r+,wb,wa,tbc,m);
    sort(r,wa,wb,tbc,m);
    for (p=,rn[F(wb[])]=,i=;i<tbc;i++)
        rn[F(wb[i])]=c0(r,wb[i-],wb[i])?p-:p++;
    if (p<tbc) dc3(rn,san,tbc,p);
    else for (i=;i<tbc;i++) san[rn[i]]=i;
    for (i=;i<tbc;i++) if (san[i]<tb) wb[ta++]=san[i]*;
    if (n%==) wb[ta++]=n-;
    sort(r,wb,wa,ta,m);
    for (i=;i<tbc;i++) wv[wb[i]=G(san[i])]=i;
    for (i=,j=,p=;i<ta&&j<tbc;p++)
        sa[p]=c12(wb[j]%,r,wa[i],wb[j])?wa[i++]:wb[j++];
    for (;i<ta;p++) sa[p]=wa[i++];
    for (;j<tbc;p++) sa[p]=wb[j++];
}
void da(int str[],int sa[],int rank[],int height[],int n,int m)
{
    for (int i=n;i<n*;i++)
        str[i]=;
    dc3(str,sa,n+,m);
    int i,j,k=;
    for (i=;i<=n;i++) rank[sa[i]]=i;
    for (i=;i<n;i++)
    {
        if (k) k--;
        j=sa[rank[i]-];
        while (str[i+k]==str[j+k]) k++;
        height[rank[i]]=k;
    }
}
int a[maxn*];
int ra[maxn*],height[maxn*],sa[maxn*];

deque< pair<int,int> > q;

int solve(int n,int k)
{
    // height[2..n]
    // length k-1 min
    int ma=;
    while (!q.empty()) q.pop_back();
    for (int i=;i<k-;i++)
    {
        pair<int,int> p=make_pair(height[i+],i);
        while (!q.empty()&&q.back()>p) q.pop_back();
        q.push_back(p);
    }
    for (int i=k-;i<=n-;i++)
    {
        while (!q.empty() && q.front().second<i-(k-)) q.pop_front();
        pair<int,int> p=make_pair(height[i+],i);
        while (!q.empty()&&q.back()>p) q.pop_back();
        q.push_back(p);
        ma=max(ma,q.front().first);
    }
    return ma;
}

int main()
{
    int n,k;
    while (~scanf("%d%d",&n,&k))
    {
        for (int i=;i<n;i++) scanf("%d",&a[i]),a[i]++;
        da(a,sa,ra,height,n,);
        printf("%d\n",solve(n,k));
    }
    return ;
}