BZOJ1692: [Usaco2007 Dec]队列变换

1692: [Usaco2007 Dec]队列变换

Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 1982 Solved: 829
[Submit][Status][Discuss]

Description

FJ
打算带他的N(1 <= N <=
30,000)头奶牛去参加一年一度的“全美农场主大奖赛”。在这场比赛中，每个参赛者都必须让他的奶牛排成一列，然后领她们从裁判席前依次走过。
今年，竞赛委员会在接受队伍报名时，采用了一种新的登记规则：他们把所有队伍中奶牛名字的首字母取出，按它们对应奶牛在队伍中的次序排成一列（比如说，如
果FJ带去的奶牛依次为Bessie、Sylvia、Dora，登记人员就把这支队伍登记为BSD）。登记结束后，组委会将所有队伍的登记名称按字典序升
序排列，就得到了他们的出场顺序。
FJ最近有一大堆事情，因此他不打算在这个比赛上浪费过多的时间，也就是说，他想尽可能早地出场。于是，他打算把奶牛们预先设计好的队型重新调整一下。
FJ的调整方法是这样的：每次，他在原来队列的首端或是尾端牵出一头奶牛，把她安排到新队列的尾部，然后对剩余的奶牛队列重复以上的操作，直到所有奶牛都
被插到了新的队列里。这样得到的队列，就是FJ拉去登记的最终的奶牛队列。
接下来的事情就交给你了：对于给定的奶牛们的初始位置，计算出按照FJ的调整规则所可能得到的字典序最小的队列。

Input

* 第1行: 一个整数：N

* 第2..N+1行: 第i+1行仅有1个'A'..'Z'中的字母，表示队列中从前往后数第i 头奶牛名字的首字母

Output

* 第1..??行: 输出FJ所能得到的字典序最小的队列。每行（除了最后一行）输 出恰好80个'A'..'Z'中的字母，表示新队列中每头奶牛姓名的首 字母

Sample Input

6

A

C

D

B

C

B

输入说明:

FJ有6头顺次排好队的奶牛：ACDBCB

Sample Output

ABCBCD

输出说明:

操作数 原队列 新队列

#1 ACDBCB

#2 CDBCB A

#3 CDBC AB

#4 CDB ABC

#5 CD ABCB

#6 D ABCBC

#7 ABCBCD

HINT

Source

【题解】

当年是暴力做的。。

不难发现优先选左右小的字母，相同就一直往里找，找到同一位置第一个不同的，比一比，先走小的那个，复杂度n^2

我们可以优化“相同就一直往里找，找到同一位置第一个不同的”

把原串和原串反向连起来，求后缀数组，用rank判断即可

背了一下SA,花了不少时间，感觉很快就要忘了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <vector>
#include <cmath>
#define min(a, b) ((a) < (b) ? (a) : (b))
#define max(a, b) ((a) > (b) ? (a) : (b))
#define abs(a) ((a) < 0 ? (-1 * (a)) : (a))
template <class T>
inline void swap(T &a, T &b)
{
T tmp = a;a = b;b = tmp;
}
inline void read(int &x)
{
x = 0;char ch = getchar(), c = ch;
while(ch < '0' || ch > '9') c = ch, ch = getchar();
while(ch <= '9' && ch >= '0') x = x * 10 + ch - '0', ch = getchar();
if(c == '-') x = -x;
}

const int INF = 0x3f3f3f3f;
const int MAXN = 1000000 + 10;

struct SuffixArray
{
int s[MAXN], sa[MAXN], rank[MAXN], height[MAXN];
int t[MAXN], t2[MAXN], c[MAXN];
int n;
void clear(){n = 0;memset(sa, 0, sizeof(sa));}
void build_sa(int m)
{
int i, *x = t, *y = t2;
for(i = 0;i < m;++ i) c[i] = 0;
for(i = 0;i < n;++ i) ++ c[x[i] = s[i]];
for(i = 1;i < m;++ i) c[i] += c[i - 1];
for(i = n - 1;i >= 0;-- i) sa[--c[x[i]]] = i;
for(int k = 1;k <= n;k <<= 1)
{
int p = 0;
for(i = n - k;i < n;++ i) y[p ++] = i;
for(i = 0;i < n;++ i) if(sa[i] >= k) y[p ++] = sa[i] - k;
for(i = 0;i < m;++ i) c[i] = 0;
for(i = 0;i < n;++ i) ++ c[x[y[i]]];
for(i = 1;i < m;++ i) c[i] += c[i - 1];
for(i = n - 1;i >= 0;-- i) sa[--c[x[y[i]]]] = y[i];
swap(x, y);
p = 1;x[sa[0]] = 0;
for(i = 1;i < n;++ i)
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 build_height()
{
int i,j,k = 0;
for(i = 0;i < n;++ i) rank[sa[i]] = i;
for(i = 0;i < n;++ i)
{
if(k) -- k;
j = sa[rank[i] - 1];
while(s[i + k] == s[j + k] && i + k < n && j + k < n) ++ k;
height[rank[i]] = k;
}
}
}A;

int ma, n;

int main()
{
read(A.n);n = A.n;
for(int i = 0;i < A.n;++ i)
scanf("%s", &A.s[i]), ma = max(ma, A.s[i]);
for(int i = 0;i < A.n;++ i)
A.s[i + A.n] = A.s[A.n - i - 1];
A.n += A.n;
A.build_sa(ma + 1);A.build_height();
int i = 0, j = n, k = 0;
while(i + j < n + n)
{
if(A.rank[i] < A.rank[j]) printf("%c", A.s[i]), ++ i;
else printf("%c", A.s[j]), ++ j;
++ k;
if(k == 80) putchar('\n'), k = 0;
}
return 0;
}

BZOJ1692