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

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

BZOJ_1692_[Usaco2007 Dec]隊列變換_後綴數組