UOJ35 后缀排序

http://www.elijahqi.win/2018/01/29/uoj35/
这是一道模板题。

读入一个长度为 nn 的由小写英文字母组成的字符串，请把这个字符串的所有非空后缀按字典序从小到大排序，然后按顺序输出后缀的第一个字符在原串中的位置。位置编号为 11 到 nn。

除此之外为了进一步证明你确实有给后缀排序的超能力，请另外输出 n−1n−1 个整数分别表示排序后相邻后缀的最长公共前缀的长度。
输入格式

一行一个长度为 nn 的仅包含小写英文字母的字符串。
输出格式

第一行 nn 个整数，第 ii 个整数表示排名为 ii 的后缀的第一个字符在原串中的位置。

第二行 n−1n−1 个整数，第 ii 个整数表示排名为 ii 和排名为 i+1i+1 的后缀的最长公共前缀的长度。
样例一
input

ababa

output

5 3 1 4 2
1 3 0 2

explanation

排序后结果为：

a
aba
ababa
ba
baba

限制与约定

1≤n≤1051≤n≤105

时间限制：1s1s

空间限制：256MB256MB256MB

qwq写一写板子 发现已经不会写了 好菜好菜啊

大概重新理解了一下这个求解height的过程即 为什么可以线性搞 因为假设我现在这个后缀为i-1号 并且假设按照排序之后在他前面的后缀是k开始的 那么我们可以知道 height[rk[i-1]]=某个数

那么可以知道k+1是k去掉一个数得到 i是i-1去掉一个数得到 那么显然k+1这个后缀必须排在i的位置之前 那么他们之间 也就是按照后缀顺序排好之后的顺序里 他们之间最小匹配值也是height[rk[i-1]-1这个-1即去掉了第i-1和第k位置上的匹配 那么搞定

还有一种情况 即我i-1&k从第一位开始就不同 他们在后缀排序中的位置也不得而知 那么可以想到 他们height[rk[i-1]]本身就是0 那么下面的height[rk[i]]显然至少会比-1大 这是显然的 所以 可以o(n)最欢o(2*n)线性时间内搞定

#include<cstdio>
#include<cstring> 
#define N 110000
int sa[N],height[N],rk[N<<1],count[N],m,k,tmp[N],rk1[N];
char s[N];
int main(){
    freopen("sa.in","r",stdin);
    scanf("%s",s+1);int n=strlen(s+1);
    for (int i=1;i<=n;++i) count[s[i]]=1;
    for (int i=1;i<=255;++i) count[i]+=count[i-1];
    for (int i=1;i<=n;++i) rk[i]=count[s[i]];m=255;
    for (int p=1;k!=n;p<<=1,m=k){
        for (int i=1;i<=m;++i) count[i]=0;
        for (int i=1;i<=n;++i) ++count[rk[i+p]];
        for (int i=1;i<=m;++i) count[i]+=count[i-1];
        for (int i=n;i;--i) tmp[count[rk[i+p]]--]=i;
        for (int i=1;i<=n;++i) count[i]=0;
        for (int i=1;i<=n;++i) ++count[rk[i]];
        for (int i=1;i<=m;++i) count[i]+=count[i-1];
        for (int i=n;i;--i) sa[count[rk[tmp[i]]]--]=tmp[i];
        memcpy(rk1,rk,sizeof(rk)>>1);rk[sa[1]]=k=1;
        for (int i=2;i<=n;++i){
            if(rk1[sa[i]]!=rk1[sa[i-1]]||rk1[sa[i]+p]!=rk1[sa[i-1]+p]) ++k;
            rk[sa[i]]=k; 
        } 
    }for (int i=1;i<=n;++i) printf("%d ",sa[i]);puts("");k=0;
    for (int i=1;i<=n;++i){
        if (rk[i]==1) continue;k=k==0?0:k-1;
        while(s[i+k]==s[sa[rk[i]-1]+k]) ++k;
        height[rk[i]]=k;
    }for (int i=2;i<=n;++i) printf("%d ",height[i]);
    return 0;
}