tyvj 1860 后缀数组

真·模板题（然而还是TLE了，tyvj真是个毒瘤，输出double什么的就是 -0.00000000，这些就TLE2333）

简单的说一下本蒟蒻看了一天后缀数组的收获（这东西太神了，，，wcwc，，收获到我的膝盖233，还是不太理解这个神奇的东西）（鉴于本人太虚，只看了一下倍增，照着黄学长扒了一个板子23333）

一共有3个神奇的数组，sa(Sufiix Array)，rank，height

对于前两个数组的求法，用基数排序的方法，每次搞出长度为2^i（除了最后，，）的子串的rank，然后用这些子串（当前长度为2^i和下一个），来更新出长度为2^(i+1)的子串

所以，一共是进行了log(strlen(ch))次的。

里面有一些小细节自己YY一下就好了，（吐槽，各种各样的sa[v[rank[..]..].]太恶心了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],height[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) height[1]=0;
        else{
            int j=sa[p][rank[p][i]-1];
            while (a[i+k]==a[j+k]) k++;
            height[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()
{
    scanf("%s",ch+1);
    n=strlen(ch+1);
    for (int i=1; i<=n; i++) a[i]=ch[i]-'a'+1;
    work();
    for (int i=1; i<=n; i++)
        printf("%d ",sa[p][i]); cout<<endl;
    for (int i=1; i<=n; i++)
        printf("%d ",height[i]);
    return 0;
}