bzoj1031[JSOI2007]字符加密Cipher——SuffixArray

传送门：点击打开链接

折腾半天，终于开始后缀数组了

个人感觉吧，后缀数组是比AC自动机难拍的，但思考难度不增加太多（或者说都不容易）

然而这题是比较简单的

这道题无非就是对n个排列进行排序，问题在于不能用n^2的空间存储下每一个串，自然可以想到用后缀数组把每一个后缀排序得到rank，以此代替要求的排列

【如果不能理解，请看这一行：因为n个排列是可排序的，在它们后面加上一些字符再进行排序并不会影响最终的结果（其实这句话不完备，因为两个排列可能相同，但因为只与尾字母有关，不妨交换两个排列的顺序，那添加字符后的结果就不影响了）】

代码如下：

#include <bits/stdc++.h>
using namespace std;
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
#define red(i, a, b) for(int i = (a); i >= (b); i--)
#define ll long long

const int maxn = 400100;
char s[maxn];
int sa[maxn], t[maxn], t2[maxn], c[maxn];
int rank[maxn], height[maxn], p[maxn>>1];
int n, len;

void build_sa(int m) {
	int i, *x = t, *y = t2;
	rep(i, 0, m-1) c[i] = 0;
	rep(i, 0, n-1) c[x[i] = s[i]]++;
	rep(i, 1, m-1) c[i] += c[i-1];
	red(i, n-1, 0) sa[--c[x[i]]] = i;
	for(int k = 1; k <= n; k <<= 1) {
		int p = 0;
		rep(i, n-k, n-1) y[p++] = i;
		rep(i, 0, n-1) if (sa[i] >= k) y[p++] = sa[i]-k;
		rep(i, 0, m-1) c[i] = 0;
		rep(i, 0, n-1) c[x[y[i]]]++;
		rep(i, 0, m-1) c[i] += c[i-1];
		red(i, n-1, 0) sa[--c[x[y[i]]]] = y[i];
		swap(x, y);
		p =  1; x[sa[0]] = 0;
		rep(i, 1, n-1) x[sa[i]] = y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k] ? p-1 : p++;
		if (p>=n) break;
		m = p;
	}
}

void getheight() {
	int j, k = 0;
	rep(i, 0, n-1) rank[sa[i]] = i;
	rep(i, 0, n-1) {
		if (k) k--;
		int j = sa[rank[i]-1];
		while(s[i+k] == s[j+k]) k++;
		height[rank[i]] = k;
	}
}

int main() {
	scanf("%s", s);
	n = strlen(s); len = n;
	rep(i, 1, n) s[n+i-1] = s[i-1];
	n = n * 2;
	build_sa(256);
	getheight();
	rep(i, 0, n / 2 - 1) p[rank[i]] = s[i+len-1];
	rep(i, 0, n) if (p[i]) printf("%c", p[i]);
	printf("\n");
	return 0;
}