码迷,mamicode.com
首页 > 其他好文 > 详细

poj 3581

时间:2014-10-19 14:17:20      阅读:291      评论:0      收藏:0      [点我收藏+]

标签:style   blog   http   color   io   os   ar   使用   for   

题意:给你一个长度为n的数组A,第一个最大,要求你把它切成三段然后分别翻转,问你翻转完后的字典序最小的数组

分析:切成三段需要确定两个分割点,对于第一个分割点由于第一个数是最大的,那么只要求一下第一段翻转后字典序最小的就是答案,

求这个字典序最小第一段的方法就是对翻转后的A建立后缀数组,取第一个符合要求的就是第一个分割点(第一段长度至少1,后两段至少为2)

对于第二个分割点,不能简单的直接去最小的字典序,第二段会影响到第三段。将一个序列分两段然后分别翻转的结果可以看做,

将原序列拼接在一起后在翻转的字串。

bubuko.com,布布扣

这样就可以继续像前面一样用后缀数组搞了。

但是这题有两个坑点:

1、有多余数据,不能使用while(scanf("%d",&N))多case读入,而是直接scanf("%d",&N) 单case

2、A的数据范围没有给,是由负数的,所以在计算后缀数组时,要进行一定的修改,在找rank值的时候。

我的做法是在初始化rank的时候对A做一次离散化的处理,然后初始化rank。计算后缀数组的算法是倍增。具体的可以看代码。

bubuko.com,布布扣
 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 }
View Code

 

poj 3581

标签:style   blog   http   color   io   os   ar   使用   for   

原文地址:http://www.cnblogs.com/shuzy/p/4034818.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!