[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[n]='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');
    }
}