BZOJ_1692_[Usaco2007 Dec]队列变换_后缀数组
BZOJ_1692_[Usaco2007 Dec]队列变换_后缀数组
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
A
C
D
B
C
B
输入说明:
FJ有6头顺次排好队的奶牛:ACDBCB
Sample Output
ABCBCD
暴力。
考虑暴力怎么做,从两边找找到第一个不一样的,比较两边的字典序哪个小拿哪个。
但我们不能一次都拿走,比如AAABCAAA,拿AAAB显然不如拿AAAA,于是这个暴力决定了我们每次只能拿一个。
拿一次也可以,不妨降低比较字典序的时间,可以把原串的反串贴在后面,然后求后缀数组,比较两个串的rank即可。
代码:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 | #include <stdio.h> #include <string.h> #include <algorithm> using namespace std; #define N 600050 char s[N]; int r[N],sa[N],wa[N],wb[N],wv[N],rank[N],height[N],n,ws[N]; inline char nc() { static char buf[100000],*p1,*p2; return p1==p2&&(p2=(p1=buf)+ fread (buf,1,100000,stdin),p1==p2)?EOF:*p1++;} inline int rd() { int x=0; char c=nc(); while (c< '0' ||c> '9' ) c=nc(); while (c>= '0' &&c<= '9' ) x=(x<<3)+(x<<1)+c- '0' ,c=nc(); return x;} inline int rc() { char c=nc(); while (c< 'A' ||c> 'Z' )c=nc(); return c- 'A' +1;} void build_suffix_array() { int i,j,p,m=28,*x=wa,*y=wb,*t; for (i=0;i<m;i++) ws[i]=0; for (i=0;i<n;i++) ws[x[i]=r[i]]++; for (i=1;i<m;i++) ws[i]+=ws[i-1]; for (i=n-1;i>=0;i--) sa[--ws[x[i]]]=i; for (j=p=1;p<n;j<<=1,m=p) { for (p=0,i=n-j;i<n;i++) y[p++]=i; for (i=0;i<n;i++) if (sa[i]-j>=0) y[p++]=sa[i]-j; for (i=0;i<n;i++) wv[i]=x[y[i]]; for (i=0;i<m;i++) ws[i]=0; for (i=0;i<n;i++) ws[wv[i]]++; for (i=1;i<m;i++) ws[i]+=ws[i-1]; for (i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i]; for (t=x,x=y,y=t,i=p=1,x[sa[0]]=0;i<n;i++) { if (y[sa[i]]==y[sa[i-1]]&&y[sa[i]+j]==y[sa[i-1]+j]) x[sa[i]]=p-1; else x[sa[i]]=p++; } } for (i=1;i<n;i++) rank[sa[i]]=i; for (i=p=0;i<n-1;height[rank[i++]]=p) for (p?p++:0,j=sa[rank[i]-1];r[i+p]==r[j+p];p++); } int main() { n=rd(); int i,ln=n,k; for (i=0;i<n;i++) r[i]=rc(); r[n]=0; for (i=n+1;i<=n+n;i++) r[i]=r[n+n-i]; n<<=1; r[n++]=0; build_suffix_array(); int j=ln-1; i=0; for (k=1;k<=ln;k++) { if (rank[i]<rank[ln+ln-j]) printf ( "%c" ,r[i++]+ 'A' -1); else printf ( "%c" ,r[j--]+ 'A' -1); if (k%80==0) puts ( "" ); } return 0; } |
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步