后缀树组 一个十分巧妙的东西，实际上是一种思想

参考博客：

https://www.cnblogs.com/jinkun113/p/4743694.html

https://www.cnblogs.com/victorique/p/8480093.html

https://www.cnblogs.com/jackdong/archive/2012/10/15/2724034.html

https://www.cnblogs.com/shanchuan04/p/5324009.html

https://www.cnblogs.com/zwfymqz/p/8413523.html

后缀树组的精髓在于倍增算法、基数排序、和一些神奇的优化，巧妙的利用各个后缀之间的关系。理解了很长时间。

基数排序在后缀数组中可以在O(n)的时间内对一个二元组(p,q)进行排序，其中p是第一关键字，q是第二关键字，比其他的排序算法都要优越。

附上洛谷3809模板代码：

#include<bits/stdc++.h>
using namespace std;
typedef unsigned int ui;
typedef long long ll;
typedef unsigned long long ull;
#define pf printf
#define mem(a,b) memset(a,b,sizeof(a))
#define prime1 1e9+7
#define prime2 1e9+9
#define pi 3.14159265
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define scand(x) scanf("%llf",&x) 
#define f(i,a,b) for(int i=a;i<=b;i++)
#define scan(a) scanf("%d",&a)
#define dbg(args) cout<<#args<<":"<<args<<endl;
#define inf 0x3f3f3f3f
const int maxn=1e6+5;
int n,m,t;
using namespace std;
char s[maxn];
int y[maxn],x[maxn],c[maxn],sa[maxn],rk[maxn],height[maxn],wt[30];
void getsa()
{
    for (int i=1;i<=n;++i) ++c[x[i]=s[i]];
    //c数组是桶 
    //x[i]是第i个元素的第一关键字 
    for (int i=2;i<=m;++i) c[i]+=c[i-1]; 
    //做c的前缀和，我们就可以得出每个关键字最多是在第几名 
    for (int i=n;i>=1;--i) sa[c[x[i]]--]=i; 
    for (int k=1;k<=n;k<<=1)
    {
        int num=0;
        for (int i=n-k+1;i<=n;++i) y[++num]=i;
        //y[i]表示第二关键字排名为i的数，第一关键字的位置 
        //第n-k+1到第n位是没有第二关键字的 所以排名在最前面 
        for (int i=1;i<=n;++i) if (sa[i]>k) y[++num]=sa[i]-k;
        //排名为i的数 在数组中是否在第k位以后
        //如果满足(sa[i]>k) 那么它可以作为别人的第二关键字，就把它的第一关键字的位置添加进y就行了
        //所以i枚举的是第二关键字的排名，第二关键字靠前的先入队 
        for (int i=1;i<=m;++i) c[i]=0;
        //初始化c桶 
        for (int i=1;i<=n;++i) ++c[x[i]];
        //因为上一次循环已经算出了这次的第一关键字 所以直接加就行了 
        for (int i=2;i<=m;++i) c[i]+=c[i-1];//第一关键字排名为1~i的数有多少个 
        for (int i=n;i>=1;--i) sa[c[x[y[i]]]--]=y[i],y[i]=0;
        //因为y的顺序是按照第二关键字的顺序来排的 
        //第二关键字靠后的，在同一个第一关键字桶中排名越靠后 
        //基数排序 
        swap(x,y);
        //这里不用想太多，因为要生成新的x时要用到旧的，就把旧的复制下来，没别的意思 
        x[sa[1]]=1;num=1;
        for (int i=2;i<=n;++i)
            x[sa[i]]=(y[sa[i]]==y[sa[i-1]] && y[sa[i]+k]==y[sa[i-1]+k]) ? num : ++num;
        //因为sa[i]已经排好序了，所以可以按排名枚举，生成下一次的第一关键字 
        if (num==n) break;
        m=num;
        //这里就不用那个122了，因为都有新的编号了 
    }
}
void geth()
{
    int k=0;
    for (int i=1;i<=n;++i) rk[sa[i]]=i;
    for (int i=1;i<=n;++i)  
    {
        if (rk[i]==1) continue;//第一名height为0 
        if (k) --k;//h[i]>=h[i-1]+1;
        int j=sa[rk[i]-1];
        while (j+k<=n && i+k<=n && s[i+k]==s[j+k]) ++k;
        height[rk[i]]=k;//h[i]=height[rk[i]];
    }
}
int main()
{
    gets(s+1);
    n=strlen(s+1);m=122;
    //n表示原字符串长度，m表示字符个数，ascll('z')=122 
    //我们第一次读入字符直接不用转化，按原来的ascll码来就可以了 
    //因为转化数字和大小写字母还得分类讨论，怪麻烦的 
    getsa();
    geth();
    for(int i=1;i<=n;i++)cout<<sa[i]<<" ";
}

另一份大佬详解：

#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;

const int N = 1000100;

char s[N];
int sa[N],t1[N],t2[N],c[N],rnk[N],height[N];
//sa[i]  排名为i的是谁 -下标 
//rnk[i] i的排名 -名次 
//height[i]  排名为i的后缀与排名为i-1的后缀的最长公共前缀长度
int n,m = 130;

void get_sa() {
    int *x = t1,*y = t2;
    //基数排序 
    for (int i=0; i<m; ++i) c[i] = 0;
    for (int i=0; i<n; ++i) c[x[i] = s[i]]++;
    for (int i=1; i<m; ++i) c[i] += c[i-1];
    for (int i=n-1; i>=0; --i) sa[--c[x[i]]] = i;
    /*
    每次循环
    都将两个长度k子串合并为一个长度为2k的串
    其中前k个字符构成的子串的排名为第一关键字，后k个字符为第二关键字
    并求出合并后的字符串的排名 
    */
    //每次排名得到一共有多少个排名p，由p优化m的值。 
    for (int k=1; k<=n; k<<=1) {
        int p = 0;
        /*
        在上一轮中，第一关键字已经排好了， 此时只需要排第二关键字即可。 
        y数组是第二关键字排序的结果，存储的是2k长度的字符串的第一关键字的下标 
        n-k到n-1中所有的元素第二关键字为0 
        */
        for (int i=n-k; i<n; ++i) y[p++] = i;
        for (int i=0; i<n; ++i) if (sa[i] >= k) y[p++] = sa[i] - k;
        for (int i=0; i<m; ++i) c[i] = 0;
        for (int i=0; i<n; ++i) c[x[y[i]]]++;
        for (int i=1; i<m; ++i) c[i] += c[i-1];
        for (int i=n-1; i>=0; --i) sa[--c[x[y[i]]]] = y[i];
        swap(x,y);
        p = 1;
        x[sa[0]] = 0;
        for (int i=1; i<n; ++i)
            x[sa[i]] = (y[sa[i-1]]==y[sa[i]] && sa[i-1]+k<n && sa[i]+k<n &&    
            y[sa[i-1]+k]==y[sa[i]+k]) ? p-1 : p++;
        if (p >= n) break;//退出条件：一共有n种排名，说明有序 
        m = p;
    }
}
void get_height() {
    for (int i=0; i<n; ++i) rnk[sa[i]] = i; // sa排名为i的下标，rnk[i]下标为i的后缀的名次。 
    int k = 0; 
    height[0] = 0;  
    /*
    性质：height[rnk[i]] >= height[rnk[i-1]]-1
    理解为起始下标为i的后缀和i-1的后缀除了开头第一个字符不同，其余的相同的
    即suffix[i,n]与suffix[i-1,n]只有第一个字符不同。 
    所以suffix[i-1,n]匹配上的子串，与suffix[i,n]匹配的子串，中k-1个是相等的。
    */ 
    for (int i=0; i<n; ++i) {
        if (!rnk[i]) continue; 
        if (k) k--;
        int j = sa[rnk[i]-1]; // 前一个排名的后缀的开始位置 
        while (i+k<n && j+k<n && s[i+k]==s[j+k]) k++;
        height[rnk[i]] = k;
    }
}
int main() {
    scanf("%s",s);
    n = strlen(s);
    get_sa();
    for (int i=0; i<n; ++i) 
        printf("%d ",sa[i]+1);
    return 0;    
}