CH1812 生日礼物
题意
描述
ftiasch 18岁生日的时候,lqp18_31给她看了一个神奇的序列 A1, A2, ..., AN. 她被允许选择不超过 M 个连续的部分作为自己的生日礼物。ftiasch想要知道选择元素之和的最大值。你能帮助她吗?
输入格式
第一行两个整数N,M。
第二行N个整数A1~AN。
输出格式
一个整数,表示答案。
样例输入
5 2
2 -3 2 -1 2
样例输出
5
数据范围与约定
对于100%的数据: N,M<=10^5, |Ai|<=10^4
来源
北京大学
分析
参照solider98的题解。
设输入序列为A[1...n], 将A的首尾的连续的非正数去除的序列B[1...m], 将序列B所有相邻的正数合并, 所有相邻的非正数合并.得到序列C[1...p], C[1...p]为正, 非正交替的序列且C[1] > 0, C[p] > 0, C[1...p]中正数的个数大于M, 设C[i]为C[1...p]中绝对值最小的数, 那么有如下结论成立:
结论1: 如果i = 1或i = p那么一定存在不包含C[i]和与C[i]相邻的元素的最优选择方案(所选则的至多M个连续的部分的和最大的选择方案, 根据C[i]的绝对值最小, 且与C[i]相邻的元素必为绝对值不小于C[i]的负数较易证明此结论)
结论2: 如果1 < i < p, 那么存在最优的选择方案同时包含C[i - 1], C[i], C[i + 1]
证明:
(1)当C[i] > 0时, 假设存在最优方案R只包含三者中的一个, 显然不可能只包含Ci - 1或Ci + 1, 即只包含C[i], 由于C中正数的个数大于M, 故可使用任意未被选择的整数C[j]替换C[i], 替换之后所有选择的数的和不会减小.假设存在最优方案R只包含三者中的两个, 显然不能为C[i - 1]和C[i + 1], 不妨设包含C[i]和C[i + 1], 此时去掉C[i]和C[i + 1]剩余选择的元素之和不会减小. 故C[i] > 0时, 存在最优方案同时包含C[i - 1], C[i]和C[i + 1]或同时不包含.
(2)当C[i] <= 0时, 假设存在最优方案R只包含三者中的一个, 显然不可能只包含C[i], 不妨设包含C[i - 1], 此时将C[i]和C[i - 1]包含进来, 同时选择的元素之和不会减小, 假设最优方案R只包含三者中的两个, 显然不可能包含C[i - 1]和C[i]或C[i + 1]和C[i], 即只能包含C[i - 1]和C[i + 1], 此时可包含进C[i]使得选择的连续的部分个数减少1, 因此可在C中选择任意一个之前R中未被选择的正数, 且所选择的元素之和不会减小, 故C[i] <= 0时, 存在最优方案同时包含C[i - 1], C[i]和C[i + 1]或同时不包含.
基于结论2, 可不断减少序列中正数的个数, 直至正数的个数不超过M(此时最优方案选择所有的正数)
时间复杂度\(O(N \log N)\)
代码
#include<bits/stdc++.h>
#define co const
#define il inline
#define rg register
#define fi first
#define se second
#define mp make_pair
template<class T>il T read(){
rg T data=0,w=1;
rg char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') w=-1;
ch=getchar();
}
while(isdigit(ch))
data=data*10+ch-'0',ch=getchar();
return data*w;
}
template<class T>il T read(rg T&x){
return x=read<T>();
}
typedef long long LL;
using namespace std;
co int MAX=1e5+5;
int A[MAX],N,M;
pair<int,int> he[MAX];int size;
pair<int,int> li[MAX<<1];int head,tot,nex[MAX<<1],pre[MAX<<1];
bool cmp(int a,int b){
return abs(he[a].fi)<abs(he[b].fi);
}
void hswap(int a,int b){
swap(li[he[a].se].se,li[he[b].se].se),swap(he[a],he[b]);
}
void up(int a){
for(int q=a,p=q>>1;p&&cmp(q,p);hswap(q,p),q=p,p=q>>1);
}
void down(int a){
for(int p=a,l=p<<1,r=l+1,t;
l<=size&&(cmp(l,p)||r<=size&&cmp(r,p));
t=r>size||min(l,r,cmp)==l?l:r,hswap(p,t),p=t,l=p<<1,r=l+1);
}
void hins(pair<int,int> pa){
he[++size]=pa,up(size);
}
void hdel(int a){
hswap(a,size--);if(a!=size+1) up(a),down(a);
}
void lins(int a,pair<int,int> b){
li[++tot]=b;
if(!a) nex[tot]=head,pre[head]=tot,head=tot;
else nex[tot]=nex[a],pre[tot]=a,pre[nex[a]]=tot,nex[a]=tot;
}
void ldel(int a){
if(a==head) head=nex[head],pre[head]=0;
else nex[pre[a]]=nex[a],pre[nex[a]]=pre[a];
}
int main(){
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
read(N),read(M);
for(int i=1;i<=N;++i) read(A[i]);
int be=1,end=N;
while(be<=end&&A[be]<=0) ++be; while(end>=be&&A[end]<=0) --end;
if(be>end){
puts("0"); return 0;
}
A[1]=A[be];int len=1;
for(int i=be+1;i<=end;++i)
if(A[i-1]>0&&A[i]>0||A[i-1]<=0&&A[i]<=0) A[len]+=A[i];
else A[++len]=A[i];
for(int i=1;i<=len;++i)
li[i].fi=A[i],nex[i]=i+1,pre[i]=i-1;
nex[len]=0,pre[1]=0,head=1,tot=len;
for(int i=1;i<=len;++i) li[i].se=size+1,hins(mp(li[i].fi,i));
while(size/2+1>M){
int u=he[1].fi,v=he[1].se;
if(!pre[v]) hdel(li[nex[v]].se),hdel(1),ldel(nex[v]),ldel(v);
else if(!nex[v]) hdel(li[pre[v]].se),hdel(1),ldel(pre[v]),ldel(v);
else{
int t=u+li[pre[v]].fi+li[nex[v]].fi,pp=pre[pre[v]];
hdel(li[pre[v]].se),hdel(li[nex[v]].se),hdel(1),
ldel(pre[v]),ldel(nex[v]),ldel(v),
lins(pp,mp(t,size+1)),hins(mp(t,tot));
}
}
int ans=0;
for(int i=head;i;i=nex[i]) if(li[i].fi>0) ans+=li[i].fi;
printf("%d\n",ans);
return 0;
}
CH1812 生日礼物的更多相关文章
- 【BZOJ】【1293】【SCOI2009】生日礼物
二分/堆 求一个最小的区间使得包含所有的颜色(并不一定只出现一次)$n\leq 10^6$ 我想的做法是:二分这个最小的长度(满足单调性……好久才想到QAQ),然后O(n)判断是否有可行的区间,这一步 ...
- 华为上机:Tom的生日礼物
Tom的生日礼物 描述: 四月一日快到了,Tom想了个愚人的好办法——送礼物.嘿嘿,不要想的太好,这礼物可没那么简单,Tom为了愚人,准备了一堆盒子,其中有一个盒子里面装了礼物.盒子里面可以再放零个或 ...
- 【链表】BZOJ 2288: 【POJ Challenge】生日礼物
2288: [POJ Challenge]生日礼物 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 382 Solved: 111[Submit][S ...
- BZOJ2288: 【POJ Challenge】生日礼物
2288: [POJ Challenge]生日礼物 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 284 Solved: 82[Submit][St ...
- 1293: [SCOI2009]生日礼物
1293: [SCOI2009]生日礼物 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1096 Solved: 584[Submit][Statu ...
- CJOJ 2485 UVa 11991 生日礼物 / UVa 11991 Easy Problem from Rujia Liu?
CJOJ 2485 UVa 11991 生日礼物 / UVa 11991 Easy Problem from Rujia Liu? Description (原题来自刘汝佳<训练指南>Pa ...
- TurnipBit:和孩子一起动手DIY“滚动”的生日礼物
当孩子的生日来临之时,做父母的总是会为该为孩子准备什么礼物而烦恼.下面就教家长朋友们利用TurnipBit开发板DIY一份"特殊"的生日礼物,不仅能增加与孩子的互动性还能提升孩子在 ...
- BZOJ 1293: [SCOI2009]生日礼物【单调队列】
1293: [SCOI2009]生日礼物 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2534 Solved: 1383[Submit][Stat ...
- bzoj 2288 【POJ Challenge】生日礼物 双向链表+堆优化
2288: [POJ Challenge]生日礼物 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1003 Solved: 317[Submit][ ...
随机推荐
- Django RF:学习笔记(8)——快速开始
Django RF:学习笔记(8)——快速开始 安装配置 1.使用Pip安装Django REST Framework: pip install djangorestframework 2.在Sett ...
- 对java沙箱机制的一点了解
1. 引入 我们都知道,程序员编写一个Java程序,默认的情况下可以访问该机器的任意资源,比如读取,删除一些文件或者网络操作等.当你把程序部署到正式的服务器上,系统管理员要为服务器的安全承担责任, ...
- superset可视化不同算法的点击率
1. 首先我们通过superset的SQL Editor来编辑语句,语句没有写完整 2. 得到的结果为: 3. 然后点击Visualize,如图所示: 4. 因为要在图中显示不同算法的点击率,需要把d ...
- WdatePicker显示乱码
1.修改zh-cn.js内容: var $lang={ errAlertMsg: "不合法的日期格式或者日期超出限定范围,需要撤销吗?", aWeekStr: ["周&q ...
- 分布式session的管理
在分布式架构或微服务架构下,必须保证一个应用服务器上保存Session后,其它应用服务器可以同步或共享这个Session,可能会出现在A1系统登录后创建并保存Session,再次发起请求,请求被转发到 ...
- c++第二十九天
p143~p151:其他隐式类型转换1.数组转换成指针,大多数表达式自动转换成指向数组首元素的指针. 2.指针的转换. 3.转换成布尔类型,例如在if (condition) 中. 4.转换成常量. ...
- JAVA面试题整理(5)-数据库
数据库 1.Oracle/mysql分页有什么优化 2.悲观锁.乐观锁 悲观锁(Pessimistic Concurrency Control,PCC):假定会发生并发冲突,屏蔽一切可能违反数据完整性 ...
- 初探动态规划(DP)
学习qzz的命名,来写一篇关于动态规划(dp)的入门博客. 动态规划应该算是一个入门oier的坑,动态规划的抽象即神奇之处,让很多萌新 萌比. 写这篇博客的目标,就是想要用一些容易理解的方式,讲解入门 ...
- 使用 v-cloak 防止页面加载时出现 vue.js 的变量名
知识点:使用 v-cloak 防止页面加载时出现 vue.js 的变量名 场景:在使用vue语法,实现下拉框功能时,展示数据列表之前,出现对应的 vuejs 变量名 代码: var vm = new ...
- Learning Ansible with Vagrant
介绍 Vagrant介绍 Vagrant是一款构建虚拟开发环境的工具,Vagrant的使用依赖于VirtualBox,VMware等一些虚拟机,通过Vagrant能更加方便的管理虚拟机,同时还能够通过 ...