首页 > 代码库 > [NOIP2013]花匠
[NOIP2013]花匠
Description
花匠栋栋种了一排花,每株花都有自己的高度。花儿越长越大,也越来越挤。栋栋决定把这排中的一部分花移走,将剩下的留在原地,使得剩下的花能有空间长大,同时,栋栋希望剩下的花排列得比较别致。 具体而言,栋栋的花的高度可以看成一列整数h1, h2, … , hn。设当一部分花被移走后,剩下的花的高度依次为g1, g2, … , gm,则栋栋希望下面两个条件中至少有一个满足: 条件 A:对于所有的1 ≤ i < m / 2,g_2i > g_2i?1,且g_2i > g_2i+1; 条件 B:对于所有的1 ≤ i < m / 2,g_2i < g_2i?1,且g_2i < g_2i+1。 此处2i及2i-1,2i+1都为下标。 注意上面两个条件在m = 1时同时满足,当m > 1时最多有一个能满足。 请问,栋栋最多能将多少株花留在原地。
Input
输入的第一行包含一个整数 n,表示开始时花的株数。 第二行包含 n 个整数,依次为h1, h2,… , hn,表示每株花的高度。
Output
输出一行,包含一个整数 m,表示最多能留在原地的花的株数。
Sample Input
5 5 3 2 1 2
Sample Output
3
Hint
对于 20%的数据,n ≤ 10; 对于 30%的数据,n ≤ 25; 对于 70%的数据,n ≤ 1000,0 ≤ hi ≤ 1000; 对于 100%的数据,1 ≤ n ≤ 100,000,0 ≤ hi ≤ 1,000,000,所有的h_i随机生成,所有随机数服从某区间内的均匀分布。
思路:{按a,b两种情况暴力dp,O(n^2),然后用权值线段树优化即可,最终复杂度O(nlogn)。}
1 #include<iostream> 2 #include<cstdio> 3 #include<cstdlib> 4 #include<cstring> 5 #include<algorithm> 6 #define RG register 7 #define lowbit ((k)&(-k)) 8 #define rs ((o<<1)|1) 9 #define ls (o<<1) 10 #define mid ((l+r)>>1) 11 #define maxx 1000003 12 #include<queue> 13 using namespace std; 14 int a[maxx],dp[maxx],sum[maxx],n,tree1[maxx*4],tree2[maxx*4]; 15 void Insert1(int o,int l,int r,int p,int num){ 16 if(l==r){tree1[o]=num;return;} 17 if(mid<p)Insert1(rs,mid+1,r,p,num); 18 else Insert1(ls,l,mid,p,num); 19 tree1[o]=max(tree1[rs],tree1[ls]); 20 } 21 void Insert2(int o,int l,int r,int p,int num){ 22 if(l==r){tree2[o]=num;return;} 23 if(mid<p)Insert2(rs,mid+1,r,p,num); 24 else Insert2(ls,l,mid,p,num); 25 tree2[o]=max(tree2[rs],tree2[ls]); 26 } 27 int ask1(int o,int l,int r,int L,int R){ 28 if(l>=L&&r<=R)return tree1[o]; 29 if(mid<L)return ask1(rs,mid+1,r,L,R); 30 else if(mid>=R)return ask1(ls,l,mid,L,R); 31 else return max(ask1(rs,mid+1,r,L,R),ask1(ls,l,mid,L,R)); 32 } 33 int ask2(int o,int l,int r,int L,int R){ 34 if(l>=L&&r<=R)return tree2[o]; 35 if(mid<L)return ask2(rs,mid+1,r,L,R); 36 else if(mid>=R)return ask2(ls,l,mid,L,R); 37 else return max(ask2(rs,mid+1,r,L,R),ask2(ls,l,mid,L,R)); 38 } 39 int main(){ 40 scanf("%d",&n); 41 for(RG int i=1;i<=n;++i)scanf("%d",&a[i]),a[i]++; 42 dp[1]=sum[1]=1;Insert1(1,0,maxx,a[1],1); 43 for(RG int i=2;i<=n;++i){int Max1=ask1(1,0,maxx,0,a[i]-1),Max2=ask2(1,0,maxx,a[i]+1,maxx); 44 if(Max1>Max2){ 45 dp[i]=Max1+1; 46 Insert2(1,0,maxx,a[i],dp[i]); 47 } 48 else dp[i]=Max2+1,Insert1(1,0,maxx,a[i],dp[i]); 49 } 50 int ans=dp[n];memset(tree1,0,sizeof(tree1));memset(tree2,0,sizeof(tree2)); 51 memset(dp,0,sizeof(dp)); 52 dp[1]=sum[1]=1;Insert1(1,0,maxx,a[1],1); 53 for(RG int i=2;i<=n;++i){int Max1=ask1(1,0,maxx,a[i]+1,maxx),Max2=ask2(1,0,maxx,0,a[i]-1); 54 if(Max1>Max2){ 55 dp[i]=Max1+1; 56 Insert2(1,0,maxx,a[i],dp[i]); 57 } 58 else dp[i]=Max2+1,Insert1(1,0,maxx,a[i],dp[i]); 59 } 60 ans=max(ans,dp[n]);cout<<ans; 61 return 0; 62 }
[NOIP2013]花匠
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。