【BZOJ-1692&1640】队列变换后缀数组 + 贪心

1692: [Usaco2007 Dec]队列变换

Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 1153 Solved: 482
[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

Gold

Solution

后缀数组

前后都可以拿出来，所以其实需要前缀的，那么建反串好了，连在原串后面

然后思考一种贪心的策略：

题目要求字典序最小，即队首队尾取较小拿出来即可；那么如果出现类似ABBCDBBA的情况，取队首还是队尾就比较的尴尬了

只需要比较前缀后缀即可。那么用rank比较即可，非常好想好写

Code

#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
#define maxn 60010
char S[maxn];int SA[maxn],len;
int wa[maxn],wb[maxn],ws[maxn],wv[maxn];
int cmp(int *r,int a,int b,int l)
{
    return r[a]==r[b]&&r[a+l]==r[b+l];
}
void DA(char *r,int *sa,int n,int m)
{
    int p,*x=wa,*y=wb,*t;
    for (int i=; i<m; i++) ws[i]=;
    for (int i=; i<n; i++) ws[x[i]=r[i]]++;
    for (int i=; i<m; i++) ws[i]+=ws[i-];
    for (int i=n-; i>=; i--) sa[--ws[x[i]]]=i;
    p=; for (int j=; p<n; j*=,m=p)
        {
            p=; for (int i=n-j; i<n; i++) y[p++]=i;
            for (int i=; i<n; i++) if (sa[i]>=j) y[p++]=sa[i]-j;
            for (int i=; i<n; i++) wv[i]=x[y[i]];
            for (int i=; i<m; i++) ws[i]=;
            for (int i=; i<n; i++) ws[wv[i]]++;
            for (int i=; i<m; i++) ws[i]+=ws[i-];
            for (int i=n-; i>=; i--) sa[--ws[wv[i]]]=y[i];
            t=x,x=y,y=t;p=;x[sa[]]=;
            for (int i=; i<n; i++)
                x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
        }
}
int rank[maxn],height[maxn];
void calheight(char *r,int *sa,int n)
{
    int k=;
    for (int i=; i<=n; i++) rank[sa[i]]=i;
    for (int i=; i<n; height[rank[i++]]=k)
        {k?k--:;for (int j=sa[rank[i]-]; r[i+k]==r[j+k]; k++);}
}
int main()
{
    int n;
    scanf("%d\n",&n); char s[];
    for (int i=; i<=n; i++)
        scanf("%s",s),S[len++]=s[];
    S[len++]='@';
    for (int i=n-; i>=; i--) S[len++]=S[i];
    S[len]=;
    DA(S,SA,len+,); calheight(S,SA,len);
    int L=,R=n+,cnt=;
    while (cnt<n)
        {
            cnt++;
            if (rank[L]<rank[R]) putchar(S[L]),L++;
                else putchar(S[R]),R++;
            if (cnt%==) puts("");
        }
    return ;
}

轻松1A系列..话说数据有点弱..一开始手误少了条语句竟然也1A了...