1692: [Usaco2007 Dec]队列变换
Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1982 Solved: 829[][][]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
Sample Output
HINT
Source
【题解】
当年是暴力做的。。
不难发现优先选左右小的字母，相同就一直往里找，找到同一位置第一个不同的，比一比，先走小的那个，复杂度n^2
我们可以优化“相同就一直往里找，找到同一位置第一个不同的”
把原串和原串反向连起来，求后缀数组，用rank判断即可
背了一下SA,花了不少时间，感觉很快就要忘了。
1 #include2 #include 3 #include 4 #include 5 #include 6 #include 7 #include 8 #include 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 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 SuffixArray29 {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 }