P5958 【[POI2017]Sabotaż】

题意描述

在一棵以1号节点为根节点的树上,有很多纯洁的白点,

BUT,突然有一个黑点出现(可能在任意位置)

它要染黑尽可能多的节点,而在一棵子树中,

只有当黑点的比例\(>x\)才可以染黑根节点(即整棵子树)

求x的最小值,使得整棵树中被染黑的节点数不超过\(k\)个

如果你看不懂请走传送门

算法分析

一道很裸的树形DP,但思路很巧

显然本题有以下性质:

  1. 最坏情况下,最开始的叛徒是叶子结点
  2. 因为一个节点被染黑了,一起为根节点的子树将全黑,所以最终被染黑的一定是一颗子树

先设计状态:\(f(i)\)表示使得\(i\)不变黑的最小\(x\)

易得:\(f(i)\)也是使得\(i\)变黑的最大\(x\)

可知\(f(i)\)仅与\(f(son{i})\)以及\(soni\)的大小有关(这里的\(soni\)表示\(i\)的子节点)

那么我们用\(sum(i)\)表示以\(i\)为根节点的子树的大小,\(sum(i)\)是需要提前用dfs预处理的

显然\(i\)被染黑仅必须满足以下两种情况:

  1. \(f(soni)<=x\),即\(i\)的某棵子树被染黑
  2. \(sum(soni)/(sum(i)-1)<=x\),即\(soni\)的被染黑足以导致\(i\)的被染黑

根据以上规律可以推出如下方程:

\(f(i)=max(f(i),min(f(soni),sum(soni)/(sum(i)-1)))\)(不会用LaTeX写公式的蒟蒻瑟瑟发抖)

取\(min\)是因为需要同时满足条件1和条件2,

取\(max\)是因为需要答案最优

(貌似貌似到这里就结束了呢)

其实还有以下三点细节需要注意:

  1. 对于叶子结点\(a\),显然有\(f(a)=1\)
  2. 对于\(ans\)当\(sum(a)<k\)时是不需要考虑的,因为它合法,所以即使它被染黑也无所谓
  3. 针对上一条结论,易得\(ans=max(ans,f(a))\),其中\(sum(a)>=k\)

然后就去快乐\(AC\)吧......

代码实现

警告:不要直接复制代码,如果你\(RE\)了你知道是为什么的

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<vector>
#define maxn 500050
using namespace std; int n,k,sum[maxn];
double f[maxn],ans;
vector<int>v[maxn];//用vector存图笑哈哈 void dfs(int now,int fa){
sum[now]=1;
for(int i=0;i<v[now].size();i++){
int to=v[now][i];
dfs(to,now);
sum[now]+=sum[to];
}//dfs预处理sum[]
if(sum[now]==1){f[now]=1.0;return;}//叶子结点的处理
for(int i=0;i<v[now].size();i++){
int to=v[now][i];
f[now]=max(f[now],min(f[to],(double)sum[to]/(sum[now]-1)));
}//dp
if(sum[now]>k) ans=max(ans,f[now]);//统计答案,注意if
return;
} int main(){
scanf("%d %d",&n,&k);
int x;
for(int i=2;i<=n;i++){
scanf("%d",&x);
v[x].push_back(i);
}
dfs(1,0);
printf("%.8lf",ans);//注意精度问题
while(1){}//防作弊
return 0;
}

结语

安利dalao博客

(光速逃...

P5958 【[POI2017]Sabotaż】的更多相关文章

  1. 【POI2017||bzoj4726】Sabota?

    上学期putsnan过了一次,这学期认真写了一遍…… #include<bits/stdc++.h> #define N 500010 using namespace std; ]; ,n ...

  2. 【POI2017||bzoj4726】Flappy Birds

    外国人很良心的啊,这题比NOIP那题还简单…… 不用管他最后的位置,因为移动的次数肯定是恒定的,所以维护在每一个柱子的位置能飞到的范围,递推下去即可. #include<bits/stdc++. ...

  3. 【BZOJ4726】[POI2017]Sabota? 树形DP

    [BZOJ4726][POI2017]Sabota? Description 某个公司有n个人, 上下级关系构成了一个有根树.其中有个人是叛徒(这个人不知道是谁).对于一个人, 如果他 下属(直接或者 ...

  4. 【AR实验室】mulberryAR : ORBSLAM2+VVSION

    本文转载请注明出处 —— polobymulberry-博客园 0x00 - 前言 mulberryAR是我业余时间弄的一个AR引擎,目前主要支持单目视觉SLAM+3D渲染,并且支持iOS端,但是该引 ...

  5. 【.net 深呼吸】细说CodeDom(1):结构大观

    CodeDom 是啥东东?Html Dom听过吧,XML Dom听过吧.DOM一般可翻译为 文档对象模型,那 Code + DOM呢,自然是指代码文档模型了.如果你从来没接触过 CodeDom,你大概 ...

  6. 【Machine Learning】KNN算法虹膜图片识别

    K-近邻算法虹膜图片识别实战 作者:白宁超 2017年1月3日18:26:33 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现的深入理解.本系列文章是作者结 ...

  7. 【前端性能】高性能滚动 scroll 及页面渲染优化

    最近在研究页面渲染及web动画的性能问题,以及拜读<CSS SECRET>(CSS揭秘)这本大作. 本文主要想谈谈页面优化之滚动优化. 主要内容包括了为何需要优化滚动事件,滚动与页面渲染的 ...

  8. 【深入浅出jQuery】源码浅析--整体架构

    最近一直在研读 jQuery 源码,初看源码一头雾水毫无头绪,真正静下心来细看写的真是精妙,让你感叹代码之美. 其结构明晰,高内聚.低耦合,兼具优秀的性能与便利的扩展性,在浏览器的兼容性(功能缺陷.渐 ...

  9. 【.net 深呼吸】程序集的热更新

    当一个程序集被加载使用的时候,出于数据的完整性和安全性考虑,程序集文件(在99.9998%的情况下是.dll文件)会被锁定,如果此时你想更新程序集(实际上是替换dll文件),是不可以操作的,这时你得把 ...

随机推荐

  1. Java知识系统回顾整理01基础01第一个程序07Eclipse使用----找不到类如何解决?

    一.现象 有时候会碰到如图所示的问题,分明有Hello这个类,并且也有主方法,可是运行就会出现找不到或者无法加载类Hello,或者Class Not Found 异常. 出现这个状况,有多种原因造成, ...

  2. Java知识系统回顾整理01基础01第一个程序03Eclipse下载安装

    Eclipse是最流行的java 集成开发环境IDE(Integrated Development Environment) 下载安装Eclipse两种方式 一.方式1:Eclipse官网下载安装 链 ...

  3. 【奇淫巧技】sqlmap绕过过滤的tamper脚本分类汇总

    sqlmap绕过过滤的tamper脚本分类汇总

  4. Arduino Uno微控制器采用的是Atmel的ATmega328

    参考:https://www.yiboard.com/thread-831-1-1.html 在本篇文章中,我们将详细介绍Arduino开发板的硬件电路部分,具体来说,就是介绍Arduino Uno开 ...

  5. matlab中ceil朝正无穷大四舍五入

    来源:https://ww2.mathworks.cn/help/matlab/ref/ceil.html?searchHighlight=ceil&s_tid=doc_srchtitle 本 ...

  6. 《C++ primer plus》第3章练习题

    注:有的题设条件自己改动了一下,比如英寸英尺改成米和厘米,部分题目自己加了点额外要求. 1.要求用户输入身高(m),用下划线提示用户输入,将身高转化成"米"加"厘米&qu ...

  7. Python装饰器实现带参数和不带参数

    1 def log(text=None): 2 3 if isinstance(text, str): 4 def decorator(func): 5 @functools.wraps(func) ...

  8. 基于raft共识搭建的Fabric1.4.4多机网络环境

    1准备工作介绍 1各个主机ip以及节点分配情况 各个主机的节点分配情况 ip地址 orderer0.example.com,peer0.org1.example.com 172.17.3.60 ord ...

  9. spring boot:redis+lua实现生产环境中可用的秒杀功能(spring boot 2.2.0)

    一,秒杀需要具备的功能: 秒杀通常是电商中用到的吸引流量的促销活动方式 搭建秒杀系统,需要具备以下几点: 1,限制每个用户购买的商品数量,(秒杀价格为吸引流量一般会订的很低,不能让一个用户全部抢购到手 ...

  10. Mac下面 matplotlib 中文无法显示解决

    一.环境描述 python 3.7 mac 10.14.5 二.问题描述 如下图所示,当使用matplotlib绘制图片的时候,所有的中文字符无法正常显示. 三.解决方法 1.下载字体ttf文件 链接 ...