bzoj1692 队列变换

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

贪心+后缀数组。

肯定是让奶牛的顺序越往前越小。

以下摘抄http://blog.csdn.net/xueyifan1993/article/details/7773750

假设f[i]表示从第i个开始整个串的后缀 g[i]表示从i开始，整个串的前缀

那么如果遇到左右相等的情况只需要比较f[L]和g[R]的大小就可以选哪一个了。

显然如果f[L]>g[R]的话，那么假如 s[L + 0..p ] == s[R - 0..p] 相等，下一 ch[L + p + 1] > ch[R - p -1]一定会不相等，并且ch[R-p-1]更优，如果完全相等那么说明怎么取都一样。
因为这样后缀数组可以解决他

直接给串的后面加一个'0' 然后再加一个反过来的穿在后面

比如ABCD就变为

ABCD0DCBA

这样求出后缀数组的rank就有一些作用

因为ABCD

f[1]对应的串就是ABCD0DCBA 也就是rank[1]

g[4]对应的就是DCBA 也就是rank[6]

然后可以通过对比rank来得到正确解

#include<algorithm>
#include<iostream>
#include<fstream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<cmath>
#include<set>
#include<map>
#define INF 0x7fffffff
#define N 60005
using namespace std;
inline int Readint()
{
    int x=0,f=1;
	char ch = getchar();
    while(ch > '9' || ch < '0'){if(ch =='-')f = -1;ch = getchar();}
    while(ch >= '0' && ch <= '9'){x = x * 10 + ch - '0'; ch = getchar();}
    return x * f;
}
int n,tot,p,q,k;
int a[N],v[N],sa[2][N],rank[2][N];
void calsa(int sa[N],int rank[N],int SA[N],int Rank[N])
{
     for(int i = 1; i <= n; i ++)v[ rank[sa[i]] ] = i;
     for(int i = n; i >= 1; i --)
         if(sa[i] > k)
             SA[v[rank[sa[i] - k]] --] = sa[i] - k;
     for(int i = n - k + 1; i <= n; i ++)
         SA[ v[rank[i]] -- ] = i;
     for(int i = 1; i <= n; i ++)
         Rank[SA[i]] = Rank[ SA[i - 1] ] + (rank[ SA[i - 1] ] != rank[SA[i]] || rank[ SA[i - 1] + k ] != rank[ SA[i] + k ]);
}
void work()
{
    p = 0,q = 1;
	a[0] = -1;
    memset(v,0,sizeof(v));
    for(int i = 1; i <= n; i ++) v[ a[i] ] ++;
    for(int i = 1; i <= 256; i ++)v[i] += v[i - 1];
    for(int i = 1; i <= n; i ++)
       sa[p][ v[a[i]] -- ] = i;
    for(int i = 1; i <= n; i ++)
        rank[p][sa[p][i]] = rank[p][sa[p][i - 1]] + (a[sa[p][i]] != a[sa[p][i - 1]]);
    k = 1; 
    while(k < n)
    {
        calsa(sa[p],rank[p],sa[q],rank[q]);
        p ^= 1;
		q ^= 1;
		k <<= 1;
    }
}
void getans()
{
    int l = 1, r = (n >> 1) + 2;
    for(int i = 1; i <= (n >> 1); i ++)
    {
        if(rank[p][l] < rank[p][r])
            printf("%c",a[l ++] + 'A' - 1);
        else printf("%c",a[r ++] + 'A' - 1);
        if(!(i % 80))printf("\n"); 
    }
}
int main()
{
    n = Readint();
    for(int i = 1; i <= n; i ++)
    {
        char ch;
        cin>>ch;
        a[++ tot] = (int)ch - 'A' + 1;
    }
    a[++ tot] = 0;
    for(int i = n; i; i --)
        a[++ tot] = a[i];
    n = tot;
    work();
    getans();
    return 0;
}