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