http://poj.org/problem?id=3581

給一串數，將其分成三個區間並且顛倒這三個區間，使得新數列字典序最小。

參考：http://blog.csdn.net/libin56842/article/details/46417301

（不是論文題，但是還是借鑒了題解，我是真的菜）

首先我們分成的三個區間肯定是一個前綴+中間一段+一個後綴。

那麽我們先求前綴，方法很簡單，只需要將原數列倒著存，然後找字典序最小的後綴即可。可以用後綴數組實現。

（因為原數列前綴的顛倒肯定是原數列顛倒後的其中一個後綴）

再考慮如何求中間一段。我們先把前面求完的部分扣掉，然後求字典序最小的後綴……？

並不是這樣，比如3 5 0 5 0 5，按照這個方法我們得到的是0 5，但實際上我們要得到0 5 0 5。

所以為了避免這種情況，我們將這個數列倍長，這樣就能使得0505小於05了。

最後再把剩下的一點輸出即可。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cctype>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
using
 
 namespace std;
const int N=4e5+10;
int n,m,t,rank[N],sa[N],height[N],w[N],b[N],s[N];
inline bool pan(int *x,int i,int j,int k){
    int ti=i+k<n?x[i+k]:-1;
    int tj=j+k<n?x[j+k]:-1;
    return x[i]==x[j]&&ti==tj;
}
inline void SA_init(){
    int *x=rank,*y=height,r=N-10;
    
 
for(int i=0;i<r;i++)w[i]=0;
    for(int i=0;i<n;i++)w[s[i]]++;
    for(int i=1;i<r;i++)w[i]+=w[i-1];
    for(int i=n-1;i>=0;i--)sa[--w[s[i]]]=i;
    r=1;x[sa[0]]=0;
    for(int i=1;i<n;i++)
    x[sa[i]]=s[sa[i]]==s[sa[i-1]]?r-1:r++;
    for(int k=1;r<n;k<<=1){
    int yn=0;
    for(int i=n-k;i<n;i++)y[yn++]=i;
    for(int i=0;i<n;i++)
        if(sa[i]>=k)y[yn++]=sa[i]-k;
    for(int i=0;i<r;i++)w[i]=0;
    for(int i=0;i<n;i++)++w[x[y[i]]];
    for(int i=1;i<r;i++)w[i]+=w[i-1];
    for(int i=n-1;i>=0;i--)sa[--w[x[y[i]]]]=y[i];
    swap(x,y);r=1;x[sa[0]]=0;
    for(int i=1;i<n;i++)
        x[sa[i]]=pan(y,sa[i],sa[i-1],k)?r-1:r++;
    }
    for(int i=0;i<n;i++)rank[i]=x[i];
}
inline void LSH(){
    sort(b,b+m);
    m=unique(b,b+m)-b;
    for(int i=0;i<n;i++){
        s[i]=lower_bound(b,b+m,s[i])-b;
    }
    return;
}
int main(){
    scanf("%d",&n);
    for(int i=n-1;i>=0;i--){
    scanf("%d",&b[i]);
    s[i]=b[i];
    }
    m=n;
    LSH();
    SA_init();
    for(t=0;;t++){
    if(sa[t]>1)break;
    }
    for(int i=sa[t];i<n;i++)
    printf("%d\n",b[s[i]]);
    for(int i=0;i<sa[t];i++)
    s[i+sa[t]]=s[i];
    n=sa[t]*2;
    SA_init();
    for(t=0;;t++){
    if(sa[t]==0)continue;
    if(sa[t]<(n>>1))break;
    }
    for(int i=sa[t];i<(n>>1);i++)
    printf("%d\n",b[s[i]]);
    for(int i=0;i<sa[t];i++)
    printf("%d\n",b[s[i]]);
    return 0;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+歡迎訪問我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++

POJ3581：Sequence——題解