I-Identical Day[题解]
Identical Day
题目大意
给定一个长度为 \(n\) 的 \(01\) 串,对于每段长度为 \(l\) 的连续的 \(1\) ,其权值为 \(\frac{l\times(l+1)}{2}\) ,对于给出的 \(01\) 串显然可得其权值和。现在你能够将该串中任意的 \(1\) 改为 \(0\) ,问最少改多少次你能使得新串的权值和不超过 \(k\) 。
分析
首先,对于此题,需要知道,对于一个连续 \(1\) 串,在其中更改的时候如何使得效率更高?要到达最好的效果,我们必须要让该串在更改后被分成的几个连续 \(1\) 串含有 \(1\) 的数量尽量平均。
例如:
11111 \\只能更改两个1,如何更改使得效果最好?
10101 \\更改前,权值为15,现在权值为3
11111111111111111 \\只能更改两个1,如何使得更改效果最好?
11101110111101111 \\更改前,权值为105,现在权值为32
注意要尽量平均,如第二个例子,就不能是 \(11101110111011111\) ,这样权值是 \(33\) 。
知道了这一点,接下来我们想想这道题怎么做?
一开始,我们或许会想到使用一个大根堆,将每一段入堆,每次取出最长的一段来平均分,分出来的小段又分别入队。
这么贪心显然是错误的,错在哪里?有可能我们的方案并非最优。
例如:
1111111 \\设其为原串中的一段连续1串
1110111 \\第一次更改
1010111 \\按照我们上述的贪心方案第二次更改会变成这样
1011011 \\实际上更改两个最优应该是这样
那应该怎么办?
保留最开始的大根堆思路,我们可以将一个 \(01\) 串内的每个连续 \(1\) 串单独包装,计算其砍一刀的贡献,砍两刀的贡献(注意砍两刀的贡献是相对于已经砍了一刀的权值)。我们可以重载大根堆运算为贡献大的优先,这样,我们一刀一刀砍,肯定是最优的。
这么做的好处是:
1111111 \\设其为原串中的一段连续1串
1110111 \\第一次更改
1011011 \\若我们发现该段还需要二次更改我们不会管第一次更改,直接在最开始的基础上更改两次,得到最优方案
CODE
#include <bits/stdc++.h>
using namespace std;
const int N=1e5+10;
inline int read()
{
int s=0,w=1;
char ch=getchar();
while(ch<'0'||ch>'9') { if(ch=='-') w*=-1; ch=getchar(); }
while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
return s*w;
}
struct node{
int l; //该段长度
int v; //当前阶段改变数
int nex; //至下一阶段变化量
int val;
bool operator < (const node &x)const{ //重载运算符
if(nex<x.nex) return 1;
return 0;
}
};
int n,k,ans;
int unhappy; //储存目前的不开心值
int a[N];
priority_queue<node> q;
inline void solve()
{
while(unhappy>k){
node now=q.top(),add;
q.pop(); //出堆
ans+=1; //更新答案
unhappy-=now.nex; //更新不快乐值
int ad=(now.l-now.v-1)/(now.v+2),lv=(now.l-now.v-1)%(now.v+2),temp=0;
for(int i=1;i<=lv;i++) temp=temp+(ad+1)*(ad+2)/2;
for(int i=1;i<=now.v+2-lv;i++) temp=temp+ad*(ad+1)/2;
add.val=temp;
add.nex=now.val-temp;
add.l=now.l;
add.v=now.v+1;
q.push(add);
}
}
int main()
{
n=read(),k=read();
for(int i=1;i<=n;i++){
char x=getchar();
a[i]=x-'0';
}
int now=0;
for(int i=1;i<=n+1;i++){
if(a[i]) now++;
else if(now){
int temp=now*(now+1)/2; //记录该段在不改变情况下对不开心值的贡献
int l1=(now-1)/2,l2=(now-1)-l1; //最佳的两段
unhappy+=temp;
node add; //处理入堆元素
add.l=now,now=0;
add.v=1;
add.nex=temp-(l1*(l1+1)+l2*(l2+1))/2; //改变一次对答案的最大贡献
add.val=(l1*(l1+1)+l2*(l2+1))/2;
q.push(add); //入大根堆
}
}
solve();
printf("%d\n",ans);
return 0;
}
I-Identical Day[题解]的更多相关文章
- UVALive 6124 Hexagon Perplexagon 题解
http://vjudge.net/problem/viewProblem.action?id=37480 East Central Regional Contest Problem C: Hexag ...
- leetcode & lintcode 题解
刷题备忘录,for bug-free 招行面试题--求无序数组最长连续序列的长度,这里连续指的是值连续--间隔为1,并不是数值的位置连续 问题: 给出一个未排序的整数数组,找出最长的连续元素序列的长度 ...
- Codeforces Round #469 Div. 2题解
A. Left-handers, Right-handers and Ambidexters time limit per test 1 second memory limit per test 25 ...
- AtCoder Beginner Contest 089完整题解
A - Grouping 2 Time limit : 2sec / Memory limit : 256MB Score : 100 points Problem Statement There a ...
- ICPC — International Collegiate Programming Contest Asia Regional Contest, Yokohama, 2018–12–09 题解
目录 注意!!此题解存在大量假算法,请各位巨佬明辨! Problem A Digits Are Not Just Characters 题面 题意 思路 代码 Problem B Arithmetic ...
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- noip2016十连测题解
以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- Codeforces Round #353 (Div. 2) ABCDE 题解 python
Problems # Name A Infinite Sequence standard input/output 1 s, 256 MB x3509 B Restoring P ...
- 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解
题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...
随机推荐
- Yolo:实时目标检测实战(下)
Yolo:实时目标检测实战(下) YOLO:Real-Time Object Detection After a few minutes, this script will generate all ...
- Amazon SageMaker和NVIDIA NGC加速AI和ML工作流
Amazon SageMaker和NVIDIA NGC加速AI和ML工作流 从自动驾驶汽车到药物发现,人工智能正成为主流,并迅速渗透到每个行业.但是,开发和部署AI应用程序是一项具有挑战性的工作.该过 ...
- ISP算法高水平分析(下)
ISP算法高水平分析(下) 十.LSC(Lens Shade Correction)------镜头阴影矫正 Lens Shading指画面四角由于入射光线不足形成的暗角,同时,由于不同频率的光折射 ...
- 自监督学习(Self-Supervised Learning)多篇论文解读(下)
自监督学习(Self-Supervised Learning)多篇论文解读(下) 之前的研究思路主要是设计各种各样的pretext任务,比如patch相对位置预测.旋转预测.灰度图片上色.视频帧排序等 ...
- 不是都需要ARM吗?
不是都需要ARM吗? ARM系统架构简介 什么是ARM处理器,为什么没有听说? ARM-缩写:Advanced RISC Machines 该处理器起源于1984年的英格兰.在成立之初,ARM代表Ac ...
- Windows下Qt VS 打包程序 到他人电脑安装运行出现的问题
1.可能缺程序依赖的Qt动态库 ------> 使用Qt自带的windeployqt进入安装程序所在的文件夹内进行自动配置 将程序安装在C盘之外的盘,这样可以方便windeployqt ...
- 【NX二次开发】Block UI 表达式
属性说明 常规 类型 描述 BlockID String 控件ID Enable Logical 是否可操作 Group Logical ...
- 并发王者课-铂金2:豁然开朗-“晦涩难懂”的ReadWriteLock竟如此妙不可言
欢迎来到<并发王者课>,本文是该系列文章中的第15篇. 在上篇文章中,我们介绍了Java中锁的基础Lock接口.在本文中,我们将介绍Java中锁的另外一个重要的基本型接口,即ReadWri ...
- Linux中测试网络命令
ping IP -t 是持续性查看网络状态
- OO unit1 summary
Unit 1 summary 一.前言 三周左右的学习,OO第一单元顺利结束了,个人认为有必要写个blog来反思总结一下自己第一单元的学习情况,以便更好地进行后面的学习. 之前从来没有写blog的习惯 ...