poj 3581
题意：给你一个长度为n的数组A，第一个最大，要求你把它切成三段然后分别翻转，问你翻转完后的字典序最小的数组
分析：切成三段需要确定两个分割点，对于第一个分割点由于第一个数是最大的，那么只要求一下第一段翻转后字典序最小的就是 答案，
求这个字典序最小第一段的方法就是对翻转后的A建立后缀数组，取第一个符合要求的就是第一个分割点（第一段长度至少1，后两段至少为2）
对于第二个分割点，不能简单的直接去最小的字典序，第二段会影响到第三段。将一个序列分两段然后分别翻转的结果可以看做，
将原序列拼接在一起后在翻转的字串。
这样就可以继续像前面一样用后缀数组搞了。
但是这题有两个坑点：
1、有多余数据，不能使用while(scanf("%d",&N))多case读入，而是直接scanf("%d",&N) 单case
2、A的数据范围没有给，是由负数的，所以在计算后缀数组时，要进行一定的修改，在找rank值的时候。
我的做法是在初始化rank的时候对A做一次离散化的处理，然后初始化rank。计算后缀数组的算法是倍增。具体的可以看代码。
1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 #include <algorithm> 5 #include <queue> 6 #include <string> 7 #include <vector> 8 #include <set> 9 #define maxn 200000*2+10 10 #define maxm 2010 11 #define maxt 110 12 #define INF 0x3f3f3f3f 13 #define mod 1009 14 #define MAX_STATE 2500 15 using namespace std; 16 int n,N,k; 17 int Rank[maxn]; 18 int A[maxn]; 19 int B[maxn]; 20 int SS[maxn]; 21 int temp[maxn]; 22 int SA[maxn]; 23 //比较(rank[k,i],rank[k,i+k])与(rank[k,j],rank[k,j+k]) 24 bool cmp_sa(int i,int j){ 25 if(Rank[i]!=Rank[j])return Rank[i]<Rank[j]; 26 else{ 27 int ri = i+k<=n?Rank[i+k]:-1; 28 int rj = j+k<=n?Rank[j+k]:-1; 29 return ri<rj; 30 } 31 } 32 void construct_sa(int *S,int len, int *sa){ 33 n = len ; 34 //长度为1的字符串,Rank直接取字符的编码 35 vector<int>v; 36 for(int i=0;i<=n;++i)v.push_back(S[i]); 37 sort(v.begin(),v.end()); 38 v.erase(unique(v.begin(),v.end()),v.end()); 39 for(int i=0;i<=n;++i){ 40 sa[i]=i; 41 int pos = lower_bound(v.begin(),v.end(),S[i])-v.begin(); 42 Rank[i]=i<n?pos:-1; 43 } 44 //利用对长度为k的字符串排序计算长度位2k的顺序 45 for(k=1;k<=n;k*=2){//注意这里不是 int k 46 sort(sa,sa+n+1,cmp_sa); 47 //计算新的rank，暂存到temp中 48 temp[sa[0]]=0; 49 //调整rank,相同的字符串的rank时一样的 50 for(int i=1;i<=n;++i){ 51 temp[sa[i]] = temp[sa[i-1]]+ (cmp_sa(sa[i-1],sa[i])?1:0); 52 } 53 //存回rank 54 for(int i=0;i<=n;++i){ 55 Rank[i]=temp[i]; 56 } 57 } 58 } 59 void solve(){ 60 reverse_copy(A,A+N,B); 61 construct_sa(B,N,SA); 62 int p1; 63 for(int i=0;i<=N;++i){ 64 p1 = N-SA[i]; 65 if(N-p1>=2&&p1>=1)break; 66 } 67 int m = N-p1; 68 reverse_copy(A+p1,A+N,B); 69 reverse_copy(A+p1,A+N,B+m); 70 construct_sa(B,m*2,SA); 71 int p2; 72 for(int i=0;i<=2*m;++i){ 73 p2 = p1+m-SA[i]; 74 if(p2-p1>=1&&N-p2>=1)break; 75 } 76 reverse(A,A+p1); 77 reverse(A+p1,A+p2); 78 reverse(A+p2,A+N); 79 for(int i=0;i<N;++i)printf("%d\n",A[i]); 80 } 81 82 int main(){ 83 //freopen("in.txt","r",stdin); 84 //freopen("out.txt","w",stdout); 85 ios::sync_with_stdio(false); 86 scanf("%d",&N); 87 for(int i=0;i<N;++i)scanf("%d",&A[i]); 88 solve(); 89 }
