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

6
A
C
D
B
C
B

输入说明:

FJ有6头顺次排好队的奶牛：ACDBCB

Sample Output

ABCBCD

暴力。

考虑暴力怎么做，从两边找找到第一个不一样的，比较两边的字典序哪个小拿哪个。

但我们不能一次都拿走，比如AAABCAAA，拿AAAB显然不如拿AAAA，于是这个暴力决定了我们每次只能拿一个。

拿一次也可以，不妨降低比较字典序的时间，可以把原串的反串贴在后面，然后求后缀数组，比较两个串的rank即可。

代码：

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

}