bzoj1692

http://www.elijahqi.win/2017/07/14/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

显然不考虑后缀树组，这个题目的做法显然是类似贪心的方法来搜索

从左边，右边选一个ASCII码最小的一个先丢过去。但是如果AABCAA这样的串就会出现不知道先选哪一个的情况了

假设f[i]表示从第i个开始整个串的后缀 (上面的例子：f[1] = AABCAA f[3]=BCAA)

g[i]表示从i开始，他的前缀(上面的例子 g[4] = CBAA g[3]=BAA)

那么我们在贪心的时候，碰到相同的情况下，只需要比较f[L]和g[R]的大小就可以选哪一个了。

为什么呢……

显然如果f[L]>g[R]的话，那么假如 s[L + 0..p ] == s[R - 0..p] 相等，下一s[L + p + 1] > s[R - p -1]一定会不相等，并且s[R-p-1]更优，如果完全相等那么说明怎么取都一样。

后缀数组可以解决他

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

比如ABCD就变为

ABCD0DCBA

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

因为ABCD

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

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

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

#include<cstdio>
#include<cstring>
#define NN 33000
char st,str1[2*NN];
int n,a[2*NN],count[2*NN],rank[4*NN],rank1[4*NN],s[2*NN],tmp[2*NN],sa[2*NN],k;
int main(){
    freopen("bzoj1692.in","r",stdin);
    freopen("bzoj1692.out","w",stdout);
    scanf("%d\n",&n);
    for (int i=1;i<=n;++i){
        scanf("%c\n",&st);
        str1[i]=st;
        a[i]=st;
    }
    a[n+1]=280;
    str1[n+1]='#';
    for (int i=n+2;i<=2*n+1;++i) a[i]=a[2*n+1-i+1];
    for (int i=n+2;i<=2*n+1;++i) str1[i]=str1[2*n+1-i+1];
    int nn=n;
    n=2*n+1;
    //for (int i=1;i<=2*n+1;++i) printf("%d ",a[i]);
    memset(rank,0,sizeof(rank));
    memset(rank1,0,sizeof(rank1));
    memset(s,0,sizeof(s));
    for (int i=1;i<=n;++i) s[a[i]]=1;
    for (int i=1;i<=280;++i) s[i]+=s[i-1];
    for (int i=1;i<=n;++i) rank[i]=s[a[i]];
    k=0;
    for (int j=1;k!=n;j+=j){
        memset(count,0,sizeof(count));
        for (int i=1;i<=n;++i) count[rank[i+j]]++;
        for (int i=1;i<=n;++i) count[i]+=count[i-1];
        for (int i=n;i>=1;--i) tmp[count[rank[i+j]]--]=i;
        memset(count,0,sizeof(count));
        for (int i=1;i<=n;++i) count[rank[i]]++;
        for (int i=1;i<=n;++i) count[i]+=count[i-1];
        for (int i=n;i>=1;--i) sa[count[rank[tmp[i]]]--]=tmp[i];

        memcpy(rank1,rank,sizeof(rank)/2);
        rank[sa[1]]=1;k=1;
        for (int i=2;i<=n;++i){
            if (rank1[sa[i]]!=rank1[sa[i-1]]||rank1[sa[i]+j]!=rank1[sa[i-1]+j])++k;
            rank[sa[i]]=k;
        }
    }
    //for (int i=1;i<=n;++i) printf("%d ",rank[i]);
    int i=1,j=nn+2;
    while (1){
        if (rank[i]<rank[j]) {
            printf("%c",str1[i]);++i;
        }else {
            printf("%c",str1[j]);++j;
        }
        if ((j-nn-2+i-1)%80==0) printf("\n");
        if ((j-nn-2+i-1)==nn) break;
    }
    return 0;
}