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,花了不少时间,感觉很快就要忘了。

 

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstring>
 4 #include <cstdlib>
 5 #include <algorithm>
 6 #include <queue>
 7 #include <vector>
 8 #include <cmath> 
 9 #define min(a, b) ((a) < (b) ? (a) : (b))
10 #define max(a, b) ((a) > (b) ? (a) : (b))
11 #define abs(a) ((a) < 0 ? (-1 * (a)) : (a))
12 template <class T>
13 inline void swap(T &a, T &b)
14 {
15     T tmp = a;a = b;b = tmp;
16 }
17 inline void read(int &x)
18 {
19     x = 0;char ch = getchar(), c = ch;
20     while(ch < '0' || ch > '9') c = ch, ch = getchar();
21     while(ch <= '9' && ch >= '0') x = x * 10 + ch - '0', ch = getchar();
22     if(c == '-') x = -x;
23 }
24 
25 const int INF = 0x3f3f3f3f;
26 const int MAXN = 1000000 + 10;
27 
28 struct SuffixArray
29 {
30     int s[MAXN], sa[MAXN], rank[MAXN], height[MAXN];
31     int t[MAXN], t2[MAXN], c[MAXN];
32     int n;
33     void clear(){n = 0;memset(sa, 0, sizeof(sa));}
34     void build_sa(int m)
35     {
36         int i, *x = t, *y = t2;
37         for(i = 0;i < m;++ i) c[i] = 0;
38         for(i = 0;i < n;++ i) ++ c[x[i] = s[i]];
39         for(i = 1;i < m;++ i) c[i] += c[i - 1];
40         for(i = n - 1;i >= 0;-- i) sa[--c[x[i]]] = i;
41         for(int k = 1;k <= n;k <<= 1)
42         {
43             int p = 0;
44             for(i = n - k;i < n;++ i) y[p ++] = i;
45             for(i = 0;i < n;++ i) if(sa[i] >= k) y[p ++] = sa[i] - k;
46             for(i = 0;i < m;++ i) c[i] = 0;
47             for(i = 0;i < n;++ i) ++ c[x[y[i]]];
48             for(i = 1;i < m;++ i) c[i] += c[i - 1];
49             for(i = n - 1;i >= 0;-- i) sa[--c[x[y[i]]]] = y[i];
50             swap(x, y);
51             p = 1;x[sa[0]] = 0;
52             for(i = 1;i < n;++ i)
53                 x[sa[i]] = y[sa[i - 1]] == y[sa[i]] && y[sa[i - 1] + k] == y[sa[i] + k] ? p - 1 : p ++;
54             if(p >= n) break;
55             m = p;
56         }
57     }
58     
59     void build_height()
60     {
61         int i,j,k = 0;
62         for(i = 0;i < n;++ i) rank[sa[i]] = i;
63         for(i = 0;i < n;++ i)
64         {
65             if(k) -- k;
66             j = sa[rank[i] - 1];
67             while(s[i + k] == s[j + k] && i + k < n && j + k < n) ++ k;
68             height[rank[i]] = k;            
69         }
70     }
71 }A;
72 
73 int ma, n;
74 
75 int main()
76 {
77     read(A.n);n = A.n;
78     for(int i = 0;i < A.n;++ i)
79         scanf("%s", &A.s[i]), ma = max(ma, A.s[i]);
80     for(int i = 0;i < A.n;++ i)
81         A.s[i + A.n] = A.s[A.n - i - 1];
82     A.n += A.n;
83     A.build_sa(ma + 1);A.build_height();
84     int i = 0, j = n, k = 0;
85     while(i + j < n + n)
86     {
87         if(A.rank[i] < A.rank[j]) printf("%c", A.s[i]), ++ i;
88         else printf("%c", A.s[j]), ++ j;
89         ++ k;
90         if(k == 80) putchar('\n'), k = 0;
91     }
92     return 0;
93 }
BZOJ1692

 

posted @ 2018-01-22 21:00  嘒彼小星  阅读(184)  评论(0编辑  收藏  举报