bzoj 2251: [2010Beijing Wc]外星联络

黄学长说是模板题，，然而不会。。

大概是按后缀数组的顺序（也就是字典序）来枚举一下，然后再按长度枚举一下，各种各样的暴力，，，然用在height上搞2333，不会啊

#include<bits/stdc++.h>
#define N 100005<<2
#define LL long long
#define inf 0x3f3f3f3f
using namespace std;
inline int ra()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
char ch[N];
int n,k,p,q;
int a[N],rank[2][N],sa[2][N],h[N],v[N];
void cal_sa(int sa[N], int rank[N], int Sa[N], int Rank[N])
{
    for (int i=1; i<=n; i++) v[rank[sa[i]]]=i;
    for (int i=n; i>=1; i--)
        if (sa[i]>k) Sa[v[rank[sa[i]-k]]--]=sa[i]-k;
    for (int i=n-k+1; i<=n; i++) Sa[v[rank[i]]--]=i;
    for (int i=1; i<=n; i++) Rank[Sa[i]]=Rank[Sa[i-1]]+(rank[Sa[i]]!=rank[Sa[i-1]] || rank[Sa[i]+k]!=rank[Sa[i-1]+k]);
}
void get_height()
{
    k=0;
    for (int i=1; i<=n; i++)
    {
        if (rank[p][i]==1) h[1]=0;
        else{
            int j=sa[p][rank[p][i]-1];
            while (a[i+k]==a[j+k]) k++;
            h[rank[p][i]]=k; 
            if (k>0) k--;
        }
    }
}
void work()
{
    p=0,q=1;
    for (int i=1; i<=n; i++) v[a[i]]++;
    for (int i=1; i<=30; i++) v[i]+=v[i-1];
    for (int i=1; i<=n; i++) sa[p][v[a[i]]--]=i;
    for (int i=1; i<=n; i++) rank[p][sa[p][i]]=rank[p][sa[p][i-1]]+(a[sa[p][i]]!=a[sa[p][i-1]]);
    for (k=1; k<n; k<<=1,p^=1,q^=1) cal_sa(sa[p],rank[p],sa[q],rank[q]);
    get_height();
}
int main()
{
    n=ra(); scanf("%s",ch+1);
    for (int i=1; i<=n; i++) a[i]=ch[i]-'0'+1;
    work();
    for (int i=1; i<=n; i++)
        for (int j=h[i]+1; sa[p][i]+j-1<=n; j++)
        {
            int l=i,r;
        //    for (l=i; l>=1 && h[l]>=j; l--);
            for (r=i+1; r<=n && h[r]>=j; r++);
            if (r-l>1) printf("%d
",r-l);
        }
    return 0;
}