原题目地址(牛客)

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[题解]的更多相关文章

  1. UVALive 6124 Hexagon Perplexagon 题解

    http://vjudge.net/problem/viewProblem.action?id=37480 East Central Regional Contest Problem C: Hexag ...

  2. leetcode & lintcode 题解

    刷题备忘录,for bug-free 招行面试题--求无序数组最长连续序列的长度,这里连续指的是值连续--间隔为1,并不是数值的位置连续 问题: 给出一个未排序的整数数组,找出最长的连续元素序列的长度 ...

  3. Codeforces Round #469 Div. 2题解

    A. Left-handers, Right-handers and Ambidexters time limit per test 1 second memory limit per test 25 ...

  4. AtCoder Beginner Contest 089完整题解

    A - Grouping 2 Time limit : 2sec / Memory limit : 256MB Score : 100 points Problem Statement There a ...

  5. ICPC — International Collegiate Programming Contest Asia Regional Contest, Yokohama, 2018–12–09 题解

    目录 注意!!此题解存在大量假算法,请各位巨佬明辨! Problem A Digits Are Not Just Characters 题面 题意 思路 代码 Problem B Arithmetic ...

  6. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  7. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  8. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  9. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  10. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

随机推荐

  1. Yolo:实时目标检测实战(下)

    Yolo:实时目标检测实战(下) YOLO:Real-Time Object Detection After a few minutes, this script will generate all ...

  2. Amazon SageMaker和NVIDIA NGC加速AI和ML工作流

    Amazon SageMaker和NVIDIA NGC加速AI和ML工作流 从自动驾驶汽车到药物发现,人工智能正成为主流,并迅速渗透到每个行业.但是,开发和部署AI应用程序是一项具有挑战性的工作.该过 ...

  3. ISP算法高水平分析(下)

    ISP算法高水平分析(下)  十.LSC(Lens Shade Correction)------镜头阴影矫正 Lens Shading指画面四角由于入射光线不足形成的暗角,同时,由于不同频率的光折射 ...

  4. 自监督学习(Self-Supervised Learning)多篇论文解读(下)

    自监督学习(Self-Supervised Learning)多篇论文解读(下) 之前的研究思路主要是设计各种各样的pretext任务,比如patch相对位置预测.旋转预测.灰度图片上色.视频帧排序等 ...

  5. 不是都需要ARM吗?

    不是都需要ARM吗? ARM系统架构简介 什么是ARM处理器,为什么没有听说? ARM-缩写:Advanced RISC Machines 该处理器起源于1984年的英格兰.在成立之初,ARM代表Ac ...

  6. Windows下Qt VS 打包程序 到他人电脑安装运行出现的问题

    1.可能缺程序依赖的Qt动态库   ------>    使用Qt自带的windeployqt进入安装程序所在的文件夹内进行自动配置 将程序安装在C盘之外的盘,这样可以方便windeployqt ...

  7. 【NX二次开发】Block UI 表达式

    属性说明 常规         类型 描述     BlockID     String 控件ID     Enable     Logical 是否可操作     Group     Logical ...

  8. 并发王者课-铂金2:豁然开朗-“晦涩难懂”的ReadWriteLock竟如此妙不可言

    欢迎来到<并发王者课>,本文是该系列文章中的第15篇. 在上篇文章中,我们介绍了Java中锁的基础Lock接口.在本文中,我们将介绍Java中锁的另外一个重要的基本型接口,即ReadWri ...

  9. Linux中测试网络命令

    ping IP -t 是持续性查看网络状态

  10. OO unit1 summary

    Unit 1 summary 一.前言 三周左右的学习,OO第一单元顺利结束了,个人认为有必要写个blog来反思总结一下自己第一单元的学习情况,以便更好地进行后面的学习. 之前从来没有写blog的习惯 ...