【洛谷P3809】 【模板】后缀排序（倍增）

题目背景
这是一道模板题。

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

输入输出格式
输入格式：
一行一个长度为 n 的仅包含大小写英文字母或数字的字符串。

输出格式：
一行，共n个整数，表示答案。

输入输出样例
输入样例#1：
ababa
输出样例#1：
5 3 1 4 2
说明
n <= 10^6

大概是搞明白后缀数组这个倍增的实现了
代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define ll long long
#define maxx 1000005
using namespace std;
int sa[maxx];
int tp[maxx];
int rak[maxx];
int tax[maxx];
char s[maxx];
int n,m;
void qSort()
{
    for(int i=0;i<=m;i++)tax[i]=0;
    for(int i=1;i<=n;i++)tax[rak[i]]++;
    for(int i=1;i<=m;i++)tax[i]+=tax[i-1];
    for(int i=n;i>=1;i--)sa[tax[rak[tp[i]]]--]=tp[i];
}
void suffix()
{
    m=75;
    for(int i=1;i<=n;i++)
        rak[i]=s[i]-'0',tp[i]=i;
    qSort();
    for(int w=1,p=0;p<n;m=p,w<<=1)
    {
        p=0;
        for(int i=1;i<=w;i++)tp[++p]=n-w+i;
        for(int i=1;i<=n;i++)if(sa[i]>w)tp[++p]=sa[i]-w;
        qSort();
        swap(tp,rak);
        rak[sa[1]]=p=1;
        for(int i=2;i<=n;i++)
            rak[sa[i]]=(tp[sa[i-1]]==tp[sa[i]]&&tp[sa[i-1]+w]==tp[sa[i]+w])?p:++p;
    }
}
int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    suffix();
    for(int i=1;i<n;i++)printf("%d ",sa[i]);
    printf("%d\n",sa[n]);
    return 0;
}