传送门:点击打开链接
折腾半天,终于开始后缀数组了
个人感觉吧,后缀数组是比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;
}