[BZOJ1692][Usaco2007 Dec]队列变换（贪心+后缀数组）

题目描述

传送门

题解

首先贪心的思路是很好想的，因为队首和队尾字典序小的肯定要先放。

但是遇到队首和队尾相同的情况就要考虑一下，比较当前所剩的字符串的正序和反序的字典序。比如CBAC就要考虑比较CBAC和CABC，发现第二个字典序更小，显然从队尾取更优。

我们把字符串反转后接在原串的后面，然后求sa。最后比较的时候rank数组发挥了巨大的作用。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;

const int max_n=6e4+5;

char s[max_n];
int sa[max_n],c[max_n],x[max_n],y[max_n],rank[max_n],height[max_n];
int n,m=100,len,l,r,cnt;

inline void build_sa(){
for (int i=0;i<m;++i) c[i]=0;
for (int i=0;i<n;++i) c[x[i]=s[i]]++;
for (int i=1;i<m;++i) c[i]+=c[i-1];
for (int i=n-1;i>=0;--i) sa[--c[x[i]]]=i;

for (int k=1;k<=n;k<<=1){
int p=0;
for (int i=n-k;i<n;++i) y[p++]=i;
for (int i=0;i<n;++i) if (sa[i]>=k) y[p++]=sa[i]-k;

for (int i=0;i<m;++i) c[i]=0;
for (int i=0;i<n;++i) c[x[y[i]]]++;
for (int i=1;i<m;++i) c[i]+=c[i-1];
for (int i=n-1;i>=0;--i) sa[--c[x[y[i]]]]=y[i];

swap(x,y);
p=1; x[sa[0]]=0;
for (int i=1;i<n;++i) x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&
((sa[i-1]+k>=n?-1:y[sa[i-1]+k])==(sa[i]+k>=n?-1:y[sa[i]+k]))?p-1:p++;
if (p>n) break;
m=p;
}
}

inline void build_lcp(){
for (int i=0;i<n;++i) rank[sa[i]]=i;
height[0]=0;
int k=0;
for (int i=0;i<n;++i){
if (!rank[i]) continue;
if (k) --k;
int j=sa[rank[i]-1];
while (i+k<n&&j+k<n&&s[i+k]==s[j+k]) ++k;
height[rank[i]]=k;
}
}

int main(){
scanf("%d\n",&n);
len=n;
for (int i=0;i<n-1;++i) scanf("%c\n",&s[i]);
scanf("%c",&s[n-1]);
s
='a'-1;
for (int i=0;i<n;++i) s[n+n-i]=s[i];
n=n*2+1;

build_sa();
for (int i=0;i<n;++i) rank[sa[i]]=i;

l=0; r=len+1;
for (int i=0;i<len;++i){
if (rank[l]<rank[r])
putchar(s[l++]);
else putchar(s[r++]);
if (!((i+1)%80)) putchar('\n');
}
}