BZOJ 1692: [Usaco2007 Dec]队列变换 后缀数组

1692: [Usaco2007 Dec]队列变换

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 1730  Solved: 722

[Submit][Status][Discuss]

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

和上一篇BZOJ 4278: [ONTAK2015]Tasowanie 后缀数组基本一样

把字符串倒过来贴一份，中间inf隔开

后缀数组求一下rank，就可以当成二路归并搞

通过这两个题 总结一下

二路归并字典序最小问题——》后缀数组

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<set>
#include<map>
using namespace std;

typedef long long ll;

inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
inline void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=60100;

int n;

char s
;

int buc
,sa1
,sa2
,rk1
,rk2
;
int *sa=sa1,*tp=sa2,*rk=rk1,*tmp=rk2;

void get_sa(int len)
{
register int i,k;register bool flag=0;
for(i=1;i<=len;++i)buc[s[i]]++;
for(i=1;i<=300;++i)buc[i]+=buc[i-1];
for(i=len;i;i--)sa[buc[s[i]]--]=i;
for(i=1;i<=len;++i)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);

for(k=1;k<=len;k<<=1)
{
for(i=1;i<=len;++i)buc[rk[sa[i]]]=i;
for(i=len;i;i--)if(sa[i]>k)tp[buc[rk[sa[i]-k]]--]=sa[i]-k;
for(i=len-k+1;i<=len;++i)tp[buc[rk[i]]--]=i;
for(i=1;i<=len;++i)
{
tmp[tp[i]]=tmp[tp[i-1]]+(rk[tp[i]]!=rk[tp[i-1]]||rk[tp[i]+k]!=rk[tp[i-1]+k]);
if(tmp[tp[i]]==len)flag=1;
}
swap(tmp,rk);swap(sa,tp);if(flag)break;
}
}

void solve()
{
register int i(1),j(n+2),k;//cout<<"asgF"<<endl;
while(i+j<=2*n+2)
{
for(k=1;k<=80;++k)
{rk[i]<rk[j]?(putchar(s[i]),i++):(putchar(s[j]),j++);if(i+j>2*n+2)break;}puts("");
}
}

int main()
{
n=read();
register char temp;
register int i,j,k;
for(i=1;i<=n;++i)
{
temp='0';
while(temp>'Z'||temp<'A')temp=getchar();
s[i]=temp;
}
s[1+n]=s[(n<<1)+2]=300;
for(i=n+2;i<=(n<<1)+1;++i)s[i]=s[(n<<1)+2-i];
get_sa((n<<1)+2);solve();
return 0;
}
/*
6
A
C
D
B
C
B

ABCBCD
*/