P1642 规划

P1642 规划

题意:

某地方有N个工厂,有N-1条路连接它们,且它们两两都可达。每个工厂都有一个产量值和一个污染值。现在工厂要进行规划,拆除其中的M个工厂,使得剩下的工厂依然连成一片且 总产量/总污染 的值最大。

题解:

01分数规划+树形dp
01分数规划这里就不细讲了,详细看这里
对我来说难的是树形dp(捂脸)
本题要求所选的工厂依然可以连成一片,我一开始想的是直接dfs判断,但是可以判断是否连通,但是无法求出最大值来
有什么东西既可以保证连通又可以求出最大值
额,树形dp
我们设dp[u][j]表示以u为根的子树中,有一个点数位j的连通块时的最大值。这个连通块一定是包含u的
转移方程:
f[x][j] = max(f[x][j], f[x][j - k] + f[v][k]);
v是x的儿子

代码:

#include<bits/stdc++.h>
#define N 1000005
#define INF 23333333
using namespace std;
const double eps = 0.000001;
struct edge{	int v, nxt;
}e[N];
int p[N], eid;
void init(){memset(p, -1, sizeof p);eid = 0;
}
void insert(int u, int v){e[eid].v = v;e[eid].nxt = p[u];p[u] = eid ++;
}
double f[105][105], d[105];
int vis[105], size[105], n, m, a[105], b[105];
void dfs(int x){//简单的树dp, f[i][j]表示第 i 个点, 选 j 个的方案数, 其中 i 必须选 vis[x] = 1;size[x] = 1;f[x][1] = d[x];f[x][0] = 0;for(int i = p[x]; i + 1; i = e[i].nxt){int v = e[i].v;if(vis[v]) continue;dfs(v); size[x] += size[v];for(int j = min(m, size[x]); j >= 1; j --){for(int k = 0; k < j; k ++){f[x][j] = max(f[x][j], f[x][j - k] + f[v][k]);}}}
}
int check(double x){for(int i = 1; i <= n; i ++) {d[i] = a[i] * 1.0 - b[i] * 1.0 * x;vis[i] = 0;size[i] = 0; //求 d[i] 	}for(int i = 1; i <= n; i ++)for(int j = 0; j <= m; j ++) f[i][j] = -INF;dfs(1);for(int i = 1; i <= n; i ++){if(f[i][m] > -eps) return 1; //如果有满足要求的就退出 }return 0;
}
int main(){init();scanf("%d%d", &n, &m);m = n - m;for(int i = 1; i <= n; i ++) scanf("%d", &a[i]);for(int i = 1; i <= n; i ++) scanf("%d", &b[i]);for(int i = 1; i < n; i ++){int u, v;scanf("%d%d", &u, &v);insert(u, v);insert(v, u);}	double l = 0, r = 100000;while(l + eps < r){ //二分答案 double mid = (l + r) * 1.0 / 2.0;if(check(mid)) l = mid;else r = mid;}printf("%.1f", l);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/317878.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.netcore下的微服务、容器、运维、自动化发布

微服务1.1 基本概念1.1.1 什么是微服务&#xff1f;微服务架构是SOA思想某一种具体实现。是一种将单应用程序作为一套小型服务开发的方法&#xff0c;每种应用程序都在其自己的进程中运行&#xff0c;并采用轻量级的通讯机制&#xff08;TCP&#xff09;进行通信。这…

51nod1227-平均最小公倍数【杜教筛,欧拉函数】

正题 题目链接:http://www.51nod.com/Challenge/Problem.html#problemId1227 题目大意 定义 F(a)∑i1alcm(a,i)aF(a)\frac{\sum_{i1}^a lcm(a,i)}{a}F(a)a∑i1a​lcm(a,i)​ 给出l,rl,rl,r求∑ilrF(i)\sum_{il}^rF(i)∑ilr​F(i) 解题思路 好久没做数论题了 直接拆成两个前缀…

CF1017G The Tree(统计+树剖)

CF1017G The Tree 给定一棵树&#xff0c;维护以下 \(3\) 个操作&#xff1a; 1 x 如果节点 \(x\) 为白色&#xff0c;则将其染黑。否则对这个节点的所有儿子递归进行相同操作2 x 将以节点 \(x\) 为根的子树染白。3 x 查询节点 \(x\) 的颜色 \(n,q\le 10^5\)。 \(\bigstar\text…

P1768 天路

P1768 天路 题意&#xff1a; 小X为所有的路径定义了两个值&#xff0c;Vi和Pi&#xff0c;分别表示火车线路的风景趣味度和乘坐一次的价格。现在小X想知道&#xff0c;乘客从任意一个景点开始坐火车走过的一条回路上所有的V之和与P之和的比值的最大值。以便为顾客们推荐一条…

LOJ:出纳员问题(差分约束)

解析 毒瘤题 思路的关键是利用前缀和建图&#xff0c;枚举sum[24]点值 &#xff08;其实可以二分&#xff09; 主要是细节的处理不够清晰 使下标从1开始会一下子好做起来 然后把0当做源点 差分约束一定要有源点&#xff01;&#xff01; 代码 #include<bits/stdc.h> u…

[选拔赛1]花园(矩阵快速幂),JM的月亮神树(最短路),保护出题人(斜率优化)

多年不考试&#xff0c;一夜回到解放前T1&#xff1a;花园titlesolutioncodeT2&#xff1a;月亮神树titlesolutioncodeT3&#xff1a;保护出题人titlesolutioncodeT1&#xff1a;花园 title 小 L 有一座环形花园&#xff0c;沿花园的顺时针方向&#xff0c;他把各个花圃编号为…

.NET的关于人脸识别引擎分享(C#)

最近在Github上找合适的人脸识别引擎&#xff0c;想要本地化用&#xff0c;不用开放的一些API&#xff08;比如腾讯AI、百度AI&#xff09;,有些场景联不了网&#xff0c;一开始搜索的是时候&#xff08;关键字&#xff1a;face recognition&#xff09;就找到了最出名的face_r…

UOJ#84-[UR #7]水题走四方【dp】

正题 题目链接:https://uoj.ac/problem/84 题目大意 有nnn个点的一棵树&#xff0c;111为根&#xff0c;两个人从根节点往下走&#xff08;只能从深度小的点走到深度大的点&#xff09;。 两个人每一秒都可以一条边&#xff08;也可以不移动&#xff09;&#xff0c;或者不消…

361. 观光奶牛

361. 观光奶牛 题意&#xff1a; 给一个带点权和边权的有向图&#xff0c;求图中一个 Σ点权/Σ边权 最大的环。 题解&#xff1a; 01分数规划裸题 可以拿来练练手 代码&#xff1a; #include<bits/stdc.h> #define MAXN 2005 #define inf int(1e9) #define eps (…

LOJ洛谷P3225:矿场搭建(割点、点双)

文章目录解析代码解析 我的实现不够简洁&#xff0c;导致bug越调越谜… 关键就是**找被割点分开的分量的方法&#xff1a;dfs&#xff01;&#xff01;&#xff01; 一下子就豁然开朗了 代码 #include<bits/stdc.h> using namespace std; #define ll long long const …

Docker最全教程——数据库容器化(十一)

终于按时完成第二篇。本来准备着手讲一些实践&#xff0c;但是数据库部分没有讲到&#xff0c;部分实践会存在一些问题&#xff0c;于是就有了此篇以及后续——数据库容器化。本篇将从SQL Server容器化实践开始&#xff0c;并逐步讲解其他数据库的容器化实践&#xff0c;中间再…

CF1416E-Split【dp,set】

正题 题目链接:https://www.luogu.com.cn/problem/CF1416E 题目大意 给出nnn个正整数的一个序列aia_iai​&#xff0c;你要把aia_iai​拆成两个正整数的和b2i,b2i1b_{2i},b_{2i1}b2i​,b2i1​&#xff0c;要求使得bbb的相同连续段最少。 1≤n≤5105,1≤ai≤1091\leq n\leq 5\…

[选拔赛2 NOIP2018雅礼集训 Day3 u,v,w]玩个三角形(二维差分),玩个球(状压DP+map),玩个树(树上DP)

文章目录T1&#xff1a;玩个三角形titlesolutioncodeT2&#xff1a;玩个球titlesolutioncodeT3&#xff1a;玩个树titlesolutioncodeT1&#xff1a;玩个三角形 title 题目描述 考虑一个 n ∗ n 的矩阵 A&#xff0c;初始所有元素均为 0。 执行 q 次如下形式的操作&#xff…

LOJ:蚂蚁之旅(欧拉回路)

题意描述 给你无向图的NNN个点和MMM条边&#xff0c;保证这 条边都不同且不会存在同一点的自环边&#xff0c;现在问你至少要几笔才能所有边都画一遍。&#xff08;一笔画的时候笔不离开纸&#xff09; 解析 无向图中&#xff0c;若一个图是连通的&#xff0c;且奇数度个数的…

比赛记录

CF1039 2022.10.31 B 随机化概率分析&#xff0c;发现 \(\left(\dfrac{1}{40}\right)^{4000}\) 等于 \(0\)&#xff0c;所以会蒙对一次。C 比较容易分析&#xff0c;发现若 \(c_u\oplus x\notc_v\)&#xff0c;则对这条边没有任何影响&#xff1b;若 \(c_u\oplus xc_v\)&#x…

2019 年起如何开始学习 ABP 框架系列文章-开篇有益

阅读文本大概需要 3.3 分钟。本系列文章推荐阅读地址为&#xff1a;52ABP 开发文档https://www.52abp.com/Wiki/52abp/latest/Welcome-to-52abp本文的目的是为了让刚刚接触 ABP 框架的同学或者准备接触 ABP 框架的同学&#xff0c;能够理解和搞明白 ABP 框架到底是怎么回事&…

【模板】差分约束算法

【模板】差分约束算法 题意&#xff1a; 题解&#xff1a; 模板题 算法讲解 给出一组包含 m 个不等式&#xff0c;有 n 个未知数。求任意一组满足这个不等式组的解&#xff0c;或判定无解。 连边之后跑最短路&#xff0c;保证每个连通块都没有负环即可。 也可以建源点s 0&a…

数论一之定理证明——裴蜀/威尔逊/费马/扩展欧几里得/[扩展]欧拉/[扩展]中国剩余定理,欧拉函数,逆元,剩余系,筛法

打死没想到会在H老师处学懂数论同余&#xff0c;整除模运算埃式筛法欧拉筛法最大公约数和最小公倍数辗转相除法更相减损术裴蜀定理威尔逊定理费马定理同余等价类、剩余系、缩系欧拉函数欧拉定理扩展欧拉定理区间逆元扩展欧几里得中国剩余定理扩展中国剩余定理利用以上所有知识进…

LOJ:相框(欧拉回路、分类讨论)

解析 本题是一道if的板子题 抓住关键&#xff1a;使所有点的度数全变为2 首先对于度数大于2的点&#xff0c;把它分为若干2度点&#xff08;和可能的一个单点&#xff09; 现在我们只剩下单点和二度点了 接下来分来讨论一下 若有多个连通块&#xff0c;我们要把它们变成链再拼…

CF1342F Make It Ascending(状压+求过程-求结果)

CF1342F Make It Ascending 给予一个包含 \(n\) 个元素的数组 \(a\)&#xff0c;你可以进行以下操作&#xff1a; 选择两个不同的元素 \(a_i,a_j\)(\(1 \le i,j \le n\)&#xff0c;\(i \ne j\))将 \(a_j\) 的值加上 \(a_i\)&#xff0c;并移除 \(a\) 中的第 \(i\) 个元素。 求…