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;
}