BZOJ1692: [Usaco2007 Dec]队列变换

1692: [Usaco2007 Dec]队列变换

Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 594 Solved: 246
[Submit][Status]

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

HINT

Source

Gold

题解：

后缀数组第一题。。。

原理很简单：左右字母相同时比较正着和倒着的字典序即可，可以用后缀数组做预处理，O(1)比较

代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 70000

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }

 inline bool cmp(int *r,int a,int b,int l)

 {

     return r[a]==r[b]&&r[a+l]==r[b+l];

 }

 int wa[maxn],wb[maxn],wr[maxn],sa[maxn],rank[maxn],a[maxn],n;

 void da(int *r,int *sa,int n,int m)

 {

     int i,j,p,*x=wa,*y=wb,*t;

     for(i=;i<m;i++)wr[i]=;

     for(i=;i<n;i++)wr[x[i]=r[i]]++;

     for(i=;i<m;i++)wr[i]+=wr[i-];

     for(i=n-;i>=;i--)sa[--wr[x[i]]]=i;

     for(j=,p=;p<n;j<<=,m=p)

      {

          for(p=,i=n-j;i<n;i++)y[p++]=i;

          for(i=;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;

          for(i=;i<m;i++)wr[i]=;

          for(i=;i<n;i++)wr[x[y[i]]]++;

          for(i=;i<m;i++)wr[i]+=wr[i-];

          for(i=n-;i>=;i--)sa[--wr[x[y[i]]]]=y[i];

          for(t=x,x=y,y=t,p=,i=,x[sa[]]=;i<n;i++)

          x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;

      }

 }

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     n=read();

     for(int i=;i<n;i++)

      {

          char ch=' ';

          while(ch<'A'||ch>'Z')ch=getchar();

          a[i]=int(ch-'A'+);

      }

     int tot=*n+;

     a[n]=;a[tot]=;

     for(int i=;i<=n;i++)a[n+i]=a[n-i];

     da(a,sa,tot+,);

     //for(int i=0;i<=tot;i++)cout<<i<<' '<<sa[i]<<endl;

     for(int i=;i<=tot;i++)rank[sa[i]]=i;

     //for(int i=0;i<=tot;i++)cout<<i<<' '<<rank[i]<<endl;

     int l=,r=n+;

     while(l+r-n-<n)

     {

         if(rank[l]<rank[r])putchar(a[l++]+'A'-);else putchar(a[r++]+'A'-);

         if(!((l+r-n-)%))putchar('\n');

     }

     return ;

 }

秒客网