后缀数组（SA）

学习了LRJ神犇的代码。orz。

首先真心建议了解下基数排序！！且要有一定的c++程序经验，否则程序很难看懂。

然后对着下面的程序调试（假装你已经会了算法思想）

弄个一个礼拜一下午就能学会了。

该算法基于倍增，然后错位比较，得到二元对并排序。

具体待更。

代码如下：

 #include <bits/stdc++.h>

 using namespace std;

 #define re register
 #define rep(i, a, b) for (re int i = a; i <= b; ++i)
 #define repd(i, a, b) for (re int i = a; i >= b; --i)
 #define For(i, a, b, s) for (re int i = a; i <= b; s)
 #define maxx(a, b) a = max(a, b)
 #define minn(a, b) a = min(a, b)
 #define LL long long
 #define INF (1 << 30)

 inline int read() {
     , f = ; char c = getchar();
      : f, c = getchar();
     ) + (w << ) + (c ^ '), c = getchar();
     return w * f;
 }

 ;

 char s[maxn];
 int t[maxn], t2[maxn], sa[maxn], c[maxn], n;

 void build_sa(int m) { // m表示字符集大小
     int *x = t, *y = t2; // 这样写是个技巧，可以快速交换数组（实际上交换了数组地址）
     rep(i, , m-) c[i] = ;
     rep(i, , n-) c[x[i] = s[i]]++;
     rep(i, , m-) c[i] += c[i-];
     repd(i, n-, ) sa[--c[x[i]]] = i; // 到这完成了初始字符串的基数排序
     ; k <= n; k <<= ) {
         ;
         rep(i, n-k, n-) y[p++] = i;
         rep(i, , n-) if (sa[i] >= k) y[p++] = sa[i]-k; //这两句完成了第二关键字的排序，而第二关键字为x[i]+k。
         rep(i, , m-) c[i] = ;
         rep(i, , n-) c[x[y[i]]]++;
         rep(i, , m-) c[i] += c[i-];
         repd(i, n-, ) sa[--c[x[y[i]]]] = y[i]; // 再完成第一关键字的排序
         swap(x, y); // 交换数组，原来的y数组（对应后面的x）没有用了
         p = ; x[sa[]] = ; // 根据原来的x数组（对应为y数组）修改现在的x数组
         rep(i, , n-)
             x[sa[i]] = y[sa[i]] == y[sa[i-]] && y[sa[i]+k] == y[sa[i-]+k] ? p- : p++;
         if (p == n) return; // 如果所有数两两不同，到目前为止sa一定唯一不变了，退出
         m = p; // 修改字符集大小
     }
 }

 int main() {
     scanf("%s", s);
     n = strlen(s);
     build_sa(();
     rep(i, , n-) printf();
     ;
 }