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=,f=; char ch=getchar();
     while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
     while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
     return x*f;
 }
 char ch[N];
 int n,k,p,q;
 int a[N],rank[][N],sa[][N],height[N],v[N];
 void cal_sa(int sa[N], int rank[N], int Sa[N], int Rank[N])
 {
     for (int i=; i<=n; i++) v[rank[sa[i]]]=i;
     for (int i=n; i>=; i--)
         if (sa[i]>k) Sa[v[rank[sa[i]-k]]--]=sa[i]-k;
     for (int i=n-k+; i<=n; i++) Sa[v[rank[i]]--]=i;
     for (int i=; i<=n; i++) Rank[Sa[i]]=Rank[Sa[i-]]+(rank[Sa[i]]!=rank[Sa[i-]] || rank[Sa[i]+k]!=rank[Sa[i-]+k]);
 }
 void get_height()
 {
     k=;
     for (int i=; i<=n; i++)
     {
         if (rank[p][i]==) height[]=;
         else{
             int j=sa[p][rank[p][i]-];
             while (a[i+k]==a[j+k]) k++;
             height[rank[p][i]]=k;
             if (k>) k--;
         }
     }
 }
 void work()
 {
     p=,q=;
     for (int i=; i<=n; i++) v[a[i]]++;
     for (int i=; i<=; i++) v[i]+=v[i-];
     for (int i=; i<=n; i++) sa[p][v[a[i]]--]=i;
     for (int i=; i<=n; i++) rank[p][sa[p][i]]=rank[p][sa[p][i-]]+(a[sa[p][i]]!=a[sa[p][i-]]);
     for (k=; k<n; k<<=,p^=,q^=) cal_sa(sa[p],rank[p],sa[q],rank[q]);
     get_height();
 }
 int main()
 {
     scanf("%s",ch+);
     n=strlen(ch+);
     for (int i=; i<=n; i++) a[i]=ch[i]-'a'+;
     work();
     for (int i=; i<=n; i++)
         printf("%d ",sa[p][i]); cout<<endl;
     for (int i=; i<=n; i++)
         printf("%d ",height[i]);
     return ;
 }