[COCI2016-2017#1] Mag
[COCI2016-2017#1] Mag 题解
题目描述
你将获得一棵由无向边连接的树。树上每个节点都有一个魔力值。
我们定义，一条路径的魔力值为路径上所有节点魔力值的乘积除以路径上的节点数。
例如，若一条路径包含两个魔力值分别为\(3,5\)的节点，则这条路径的魔力值为 \(3\times 5/2=7.5\)。
请你计算，这棵树上魔力值最小的路径的魔力值。
输入格式
第一行一个整数\(n\)，表示树共有\(n\)个节点，编号为\(1\dots n\)。
接下来\(n-1\)行，每行两个整数\(a_i,b_i\)，表示编号为\(a_i,b_i\)的两个节点由一条无向边连接。
接下来\(n\)行，每行一个整数\(x_i\)，表示编号为\(i\)的节点的魔力值。
输出格式
一行，一个既约分数\(p/q\)。
输入输出样例
输入 #1
2
1 2
3
4
输出 #1
3/1
输入 #2
5
1 2
2 4
1 3
5 2
2
1
1
1
3
输出 #2
1/2
说明/提示
【样例解释】
样例 1 解释
注意，路径可以只包含一个节点。
这棵树上魔力值最小的路径的包含节点\(1\)，其魔力值为\(3/1\)。
样例 2 解释
这棵树上魔力值最小的路径的包含节点\(2,4\)，其魔力值为\(1\times 1/2=1/2\)。
数据规模与约定
对于\(100\%\)的数据，\(1\le n\le 10^6\)，\(1\le a_i,b_i\le n\)，\(1\le x_i\le 10^9\)。
数据保证，\(p,q\)不会超过\(10^{18}\)。
说明
题目译自 COCI2016-2017 CONTEST #1 T4 Mag。
思路
首先来证明一点。满足题意的最优解一定为由全部是1，或仅含有1个2的链来组成的（当且仅当2两边的）。
a b
——————2——————
1 1 1 1 1 1 1 2 1 1 1 1 1 1 1
证明：
- 当\(a=b\)时，含有2的链的价值为：\(\frac{2}{a+b+1}=\frac{2}{2a +1}\)，而不含有2的链价值为\(\frac{1}{a}\)，很明显，含有二的链比含有1的链优。
- 当\(a≠b\)时，不妨设\(a＞b\)，设\(a-b=k\)，含有2的链价值为\(\frac{2}{a+b+1}=\frac{2}{2a-k +1}\)，而含1的链为\(\frac{1}{a-k}=\frac{2}{2a-2k}\)，明显含有1的链更优。在1链和1链中添加任意一个大于2的数，则都会上述情况一样，没有含1的链更优，可以以上述情况来推广。
有了上述的证明，就可以进行树形DP了。就称只含有1的链为1链，其中含了一的1链为2链。
设\(dp[i][0|1][0|1]\)的第一维表示是哪一个节点的状态。对应的价值最小。第二维表示该节点参与的链中是否含有2这个点，若含有2，则该维度对应值为0，否则为1。第三维表示分子分母，1为分子，0为分母。
\(sec[0|1][0|1]\)对应的是次小值。其中的第一维与\(dp[i]\)的第二维意义相同，第二维与其第三维相同。
\(id[0|1][0|1]\)的第一维表示该节点参与的链中是否含有2这个点，若含有2，则该维度对应值为0，否则为1。第二维表示最小一次对应次小与最小。
先来说说如何状态转移。
如下面的代码所示（主体部分）。
for(int i = 0; i < SIZ; i++) {
int next = v[now][i];
if(next == fa)
continue;
DP(next, now);
if(val[now] == 1) {
if(val[next] > 2)
continue;
double son1 = (dp[next][1][1] * 1.0) / ((dp[next][1][0] + 1) * 1.0);
double son2 = (dp[next][0][1] * 1.0) / ((dp[next][0][0] + 1) * 1.0);
double self1 = (dp[now][1][1] * 1.0) / (dp[now][1][0] * 1.0);
double self2 = (dp[now][0][1] * 1.0) / (dp[now][0][0] * 1.0);
double num1 = (sec[1][1] * 1.0) / (sec[1][0] * 1.0);
double num2 = (sec[0][1] * 1.0) / (sec[0][0] * 1.0);
if(son1 < self1) {
sec[1][1] = dp[now][1][1];
sec[1][0] = dp[now][1][0];
id[1][0] = id[1][1];
dp[now][1][1] = dp[next][1][1];
dp[now][1][0] = dp[next][1][0] + 1;
id[1][1] = next;
}
else if(son1 < num1) {
id[1][0] = next;
sec[1][1] = dp[next][1][1];
sec[1][0] = dp[next][1][0] + 1;
}
if(son2 < self2) {
sec[0][1] = dp[now][0][1];
sec[0][0] = dp[now][0][0];
id[0][0] = id[0][1];
dp[now][0][1] = dp[next][0][1];
dp[now][0][0] = dp[next][0][0] + 1;
id[0][1] = next;
}
else if(son2 < num2) {
sec[0][1] = dp[next][0][1];
sec[0][0] = dp[next][0][0] + 1;
id[0][0] = next;
}
}
else if(val[now] == 2) {
if(val[next] != 1)
continue;
double son = (dp[next][1][1] * 2.0) / ((dp[next][1][0] + 1) * 1.0);
double self = (dp[now][0][1] * 1.0) / (dp[now][0][0] * 1.0);
double num = (sec[0][1] * 1.0) / (sec[0][0] * 1.0);
if(son < self) {
sec[0][1] = dp[now][0][1];
sec[0][0] = dp[now][0][0];
id[0][0] = id[0][1];
dp[now][0][1] = dp[next][1][1] * 2;
dp[now][0][0] = dp[next][1][0] + 1;
id[0][1] = next;
}
else if(son < num) {
sec[0][1] = dp[next][1][1] * 2;
sec[0][0] = dp[next][1][0] + 1;
id[0][0] = next;
}
}
}
情况一
设当前正在遍历的节点为now，对于now的每一个子节点next，若有now的权值为1，则有：
- son1：当前子节点带有的最小1链加上父节点的价值。
- son2：当前子节点带有的最小2链加上父节点的价值。
- self1：最小1链的价值。
- self2：最小2链的价值。
- num1：次小1链的价值。
- num2：次小2链的价值。
若next的权值大于2，则没有资格更新自己的父节点。需要先判断是否有资格更新父节点的值。
if(val[next] > 2)
continue;
又有几种情况：
- 当son1 < self1时，即是当前子节点的1链可以更新最小链的1链。先使用最小1链来更新次小1链的值。在使用子节点来更新最小1链的值。代码如下。
if(son1 < self1) {
sec[1][1] = dp[now][1][1];
sec[1][0] = dp[now][1][0];
id[1][0] = id[1][1];
dp[now][1][1] = dp[next][1][1];
dp[now][1][0] = dp[next][1][0] + 1;
id[1][1] = next;
}
- 在不满足上述情况时，当前子节点仅仅只能更新次小链的值，那么就用次小链更新最小链的值。代码如下。
else if(son1 < num1) {
id[1][0] = next;
sec[1][1] = dp[next][1][1];
sec[1][0] = dp[next][1][0] + 1;
}
- 更新最小二链与最次小二链，与更新1链的方法相似。更新最小2链的代码如下。
if(son2 < self2) {
sec[0][1] = dp[now][0][1];
sec[0][0] = dp[now][0][0];
id[0][0] = id[0][1];
dp[now][0][1] = dp[next][0][1];
dp[now][0][0] = dp[next][0][0] + 1;
id[0][1] = next;
}
- 更新次小2链的价值
else if(son2 < num2) {
sec[0][1] = dp[next][0][1];
sec[0][0] = dp[next][0][0] + 1;
id[0][0] = next;
}
情况二
设当前正在遍历的节点为now，对于now的每一个子节点next，若有now的权值为2。因为这条链中必会有2，所以就不需要考虑更新1链的价值。则有：
- son：当前子节点带有的最小1链加上父节点的价值。
- self2：最小2链的价值。
- num：次小2链的价值。
若当前子节点的子节点的权值不为1，则也没有资格更新父节点。
if(val[next] != 1)
continue;
更新方式与上述相同。
但最小2链只能由儿子的1链来更新，因为一条2链中只能含有1个2，而当前节点就是2。
- 更新最小
if(son < self) {
sec[0][1] = dp[now][0][1];
sec[0][0] = dp[now][0][0];
id[0][0] = id[0][1];
dp[now][0][1] = dp[next][1][1] * 2;
dp[now][0][0] = dp[next][1][0] + 1;
id[0][1] = next;
}
- 更新次小
else if(son < num) {
sec[0][1] = dp[next][1][1] * 2;
sec[0][0] = dp[next][1][0] + 1;
id[0][0] = next;
}
现在来更新答案
又分了几种情况
- 只用最小2链。若当前父节点只含有1个子节点，则只有1条链来更新答案。
double Num = (dp[now][0][1] * 1.0) / (dp[now][0][0] * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][0][1];
ans2 = dp[now][0][0];
}
- 只用最小1链，理由与2链相同。
Num = (dp[now][1][1] * 1.0) / (dp[now][1][0] * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][1][1];
ans2 = dp[now][1][0];
}
- 最长1链加最长2链。当且仅当最小2链与最小1链来自于不同的子节点。
if(id[1][1] != id[0][1]) {
Num = ((dp[now][0][1] * dp[now][1][1]) * 1.0) / ((dp[now][0][0] + dp[now][1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][0][1] * dp[now][1][1];
ans2 = dp[now][0][0] + dp[now][1][0] - 1;
}
}
- 若不满足最小2链与最小1链来自于不同的子节点，则使用次小1，2链与最小2，1链相结合来更新答案（注意顺序）。
else {
if(id[1][0] != id[0][1]) {
Num = ((dp[now][0][1] * sec[1][1]) * 1.0) / ((dp[now][0][0] + sec[1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][0][1] * sec[1][1];
ans2 = dp[now][0][0] + sec[1][0] - 1;
}
}
if(id[0][0] != id[1][1]) {
Num = ((sec[0][1] * dp[now][1][1]) * 1.0) / ((sec[0][0] + dp[now][1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = sec[0][1] * dp[now][1][1];
ans2 = sec[0][0] + dp[now][1][0] - 1;
}
}
}
- 用最小1链与次小1链来更新答案。
Num = ((dp[now][1][1] * sec[1][1]) * 1.0) / ((dp[now][1][0] + sec[1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][1][1] * sec[1][1];
ans2 = dp[now][1][0] + sec[1][0] - 1;
}
输出
记得分子分母需要互质。
C++代码
#include <cstdio>
#include <vector>
using namespace std;
#define INF 1e8
#define Min(a, b) ((a) < (b) ? (a) : (b))
void Quick_Read(int &N) {
N = 0;
char c = getchar();
int op = 1;
while(c < '0' || c > '9') {
if(c == '-')
op = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
N = (N << 1) + (N << 3) + c - 48;
c = getchar();
}
N *= op;
}
const int MAXN = 1e6 + 5;
bool flag;
vector<int> v[MAXN];
int down[MAXN], dp[MAXN][2][2];
int val[MAXN];
int n, minn;
double ans;
int ans1, ans2;
int GCD(int a, int b) {
return b == 0 ? a : GCD(b, a % b);
}
void DP(int now, int fa) {
int sec[2][2];
int id[2][2];
for(int i = 0; i <= 1; i++)
for(int j = 0; j <= 1; j++)
sec[i][j] = id[i][j] = INF;
int SIZ = v[now].size();
for(int i = 0; i < SIZ; i++) {
int next = v[now][i];
if(next == fa)
continue;
DP(next, now);
if(val[now] == 1) {
if(val[next] > 2)
continue;
double son1 = (dp[next][1][1] * 1.0) / ((dp[next][1][0] + 1) * 1.0);
double son2 = (dp[next][0][1] * 1.0) / ((dp[next][0][0] + 1) * 1.0);
double self1 = (dp[now][1][1] * 1.0) / (dp[now][1][0] * 1.0);
double self2 = (dp[now][0][1] * 1.0) / (dp[now][0][0] * 1.0);
double num1 = (sec[1][1] * 1.0) / (sec[1][0] * 1.0);
double num2 = (sec[0][1] * 1.0) / (sec[0][0] * 1.0);
if(son1 < self1) {
sec[1][1] = dp[now][1][1];
sec[1][0] = dp[now][1][0];
id[1][0] = id[1][1];
dp[now][1][1] = dp[next][1][1];
dp[now][1][0] = dp[next][1][0] + 1;
id[1][1] = next;
}
else if(son1 < num1) {
id[1][0] = next;
sec[1][1] = dp[next][1][1];
sec[1][0] = dp[next][1][0] + 1;
}
if(son2 < self2) {
sec[0][1] = dp[now][0][1];
sec[0][0] = dp[now][0][0];
id[0][0] = id[0][1];
dp[now][0][1] = dp[next][0][1];
dp[now][0][0] = dp[next][0][0] + 1;
id[0][1] = next;
}
else if(son2 < num2) {
sec[0][1] = dp[next][0][1];
sec[0][0] = dp[next][0][0] + 1;
id[0][0] = next;
}
}
else if(val[now] == 2) {
if(val[next] != 1)
continue;
double son = (dp[next][1][1] * 2.0) / ((dp[next][1][0] + 1) * 1.0);
double self = (dp[now][0][1] * 1.0) / (dp[now][0][0] * 1.0);
double num = (sec[0][1] * 1.0) / (sec[0][0] * 1.0);
if(son < self) {
sec[0][1] = dp[now][0][1];
sec[0][0] = dp[now][0][0];
id[0][0] = id[0][1];
dp[now][0][1] = dp[next][1][1] * 2;
dp[now][0][0] = dp[next][1][0] + 1;
id[0][1] = next;
}
else if(son < num) {
sec[0][1] = dp[next][1][1] * 2;
sec[0][0] = dp[next][1][0] + 1;
id[0][0] = next;
}
}
}
double Num = (dp[now][0][1] * 1.0) / (dp[now][0][0] * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][0][1];
ans2 = dp[now][0][0];
}
Num = (dp[now][1][1] * 1.0) / (dp[now][1][0] * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][1][1];
ans2 = dp[now][1][0];
}
if(id[1][1] != id[0][1]) {
Num = ((dp[now][0][1] * dp[now][1][1]) * 1.0) / ((dp[now][0][0] + dp[now][1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][0][1] * dp[now][1][1];
ans2 = dp[now][0][0] + dp[now][1][0] - 1;
}
}
else {
if(id[1][0] != id[0][1]) {
Num = ((dp[now][0][1] * sec[1][1]) * 1.0) / ((dp[now][0][0] + sec[1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][0][1] * sec[1][1];
ans2 = dp[now][0][0] + sec[1][0] - 1;
}
}
if(id[0][0] != id[1][1]) {
Num = ((sec[0][1] * dp[now][1][1]) * 1.0) / ((sec[0][0] + dp[now][1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = sec[0][1] * dp[now][1][1];
ans2 = sec[0][0] + dp[now][1][0] - 1;
}
}
}
Num = ((dp[now][1][1] * sec[1][1]) * 1.0) / ((dp[now][1][0] + sec[1][0] - 1) * 1.0);
if(Num < ans) {
ans = Num;
ans1 = dp[now][1][1] * sec[1][1];
ans2 = dp[now][1][0] + sec[1][0] - 1;
}
}
void Read() {
ans = INF;
minn = INF;
int A, B;
Quick_Read(n);
for(int i = 1; i < n; i++) {
Quick_Read(A);
Quick_Read(B);
v[A].push_back(B);
v[B].push_back(A);
}
for(int i = 1; i <= n; i++) {
Quick_Read(val[i]);
if(val[i] == 1) {
flag = true;
}
minn = Min(minn, val[i]);
}
for(int i = 1; i <= n; i++) {
dp[i][0][1] = dp[i][1][1] = INF;
if(val[i] == 1) {
dp[i][1][1] = 1;
dp[i][0][1] = 1;
}
else if(val[i] == 2) {
dp[i][1][1] = 2;
dp[i][0][1] = 2;
}
dp[i][0][0] = dp[i][1][0] = 1;
}
}
int main() {
Read();
if(!flag) {
printf("%d/1", minn);
return 0;
}
DP(1, -1);
int gcd = GCD(ans1, ans2);
ans1 /= gcd;
ans2 /= gcd;
printf("%d/%d", ans1, ans2);
return 0;
}
[COCI2016-2017#1] Mag的更多相关文章
- 「COCI2016/2017 Contest #2」Bruza
「COCI2016/2017 Contest #2」Bruza 解题思路 : 首先对于任意时刻 \(i\) ,硬币一定移动到了深度为 \(i\) 的节点,所以第 \(i\) 时刻 Danel 一定染掉 ...
- bjwc Day1 暴力大战
今天终于有题了... 题目是COCI2016/2017 Round #4 T1一看就是NP问题,k<=50,开始想暴力,想了个n^4的,大概能过,就没去管它 T2想得太naive,丢了100分给 ...
- CI Weekly #10 | 2017 DevOps 趋势预测
2016 年的最后几个工作日,我们对 flow.ci Android & iOS 项目做了一些优化与修复: iOS 镜像 cocoapods 版本更新: fir iOS上传插件时间问题修复: ...
- 猖獗的假新闻：2017年1月1日起iOS的APP必须使用HTTPS
一.假新闻如此猖獗 刚才一位老同事 打电话问:我们公司还是用的HTTP,马上就到2017年了,提交AppStore会被拒绝,怎么办? 公司里已经有很多人问过这个问题,回答一下: HTTP还是可以正常提 ...
- iOS的ATS配置 - 2017年前ATS规定的适配
苹果规定 从2017年1月1日起,新提交的 app 不允许使用NSAllowsArbitraryLoads来绕过ATS(全称:App Transport Security)的限制. 以前为了能兼容ht ...
- 深入研究Visual studio 2017 RC新特性
在[Xamarin+Prism开发详解三:Visual studio 2017 RC初体验]中分享了Visual studio 2017RC的大致情况,同时也发现大家对新的Visual Studio很 ...
- Xamarin+Prism开发详解三：Visual studio 2017 RC初体验
Visual studio 2017 RC出来一段时间了,最近有时间就想安装试试,随带分享一下安装使用体验. 1,卸载visual studio 2015 虽然可以同时安装visual studio ...
- Microsoft Visual Studio 2017 for Mac Preview 下载+安装+案例Demo
目录: 0. 前言 1. 在线安装器 2. 安装VS 3. HelloWorld 4. ASP.NET MVC 5. 软件下载 6. 结尾 0. 前言: 工作原因,上下班背着我的雷神,一个月瘦了10斤 ...
- Create an offline installation of Visual Studio 2017 RC
Create an offline installation of Visual Studio 2017 RC 2016年12月7日 ...
- .NET Core 2.0版本预计于2017年春季发布
英文原文: NET Core 2.0 Planned for Spring 2017 微软项目经理 Immo Landwerth 公布了即将推出的 .NET Core 2.0 版本的细节,该版本预计于 ...
随机推荐
- Spring In Action 5th中的一些错误
引言 最近开始学习Spring,了解到<Spring实战>已经出到第五版了,遂打算跟着<Spring实战(第五版)>来入门Spring,没想到这书一点也不严谨,才看到第三章就发 ...
- 简单记录几个wpf学习上的问题[ObservableQueue]
我想给我的程序加一个下载队列,当我点击一个下载按钮的时候,他应该把这个插件信息(对象)加到一个队列中,然后队列里去实现下载和删除任务,下载完成则删除对象 首先我想到了Queue类型,然后我在我的vie ...
- js工厂函数
经常会遇到工厂函数这个词,不过javascript不是严格的面向对象语言,不像java,C#拥有类,所以工厂函数会变得比较模糊. 简单来讲,就是创建一个可以用来创建实例的函数,这样每一个实例都是独立的 ...
- List集合，对象根据某个相同的属性，合并另外属性
描述及实现: 1.List中有 Modular字段值有重复一样的2.Modular字段一样的话去重并且OrderAmount相加 HashMap<String,DataReport> te ...
- 使用 k8s 搭建 confluence 6.10.x 版本
将公司中已有的 confluence 服务迁移到 k8s 集群中,需要保留当前已有的数据.整体需要分为如下几个步骤: 备份 mysql 数据 备份 confluence 安装目录 备份 conflue ...
- Java学习的第四十九天
1.构造函数 import java.util.Scanner; public class Cjava { public static void main(String[]args) { Time t ...
- ModelSerializer 高级使用
前言 ModelSerializer中还具有一些高级用法,如批量更新.批量删除.批量创建等. 但是批量过来的数据格式都需要与前端做好协商,什么样的数据格式是单条操作,什么样的数据格式是批量操作. 如下 ...
- Java线程池原理及分析
线程池是很常用的并发框架,几乎所有需要异步和并发处理任务的程序都可用到线程池. 使用线程池的好处如下: 降低资源消耗:可重复利用已创建的线程池,降低创建和销毁带来的消耗: 提高响应速度:任务到达时,可 ...
- 力扣 - 19. 删除链表的倒数第N个节点
目录 题目 思路1 代码实现 思路2 代码实现 题目 给定一个链表,删除链表的倒数第 n 个节点,并且返回链表的头结点. 示例: 给定一个链表: 1->2->3->4->5, ...
- 重拾python所要知道的一些主干知识点
前言:因为有一段时间没有用python了,最近需要用到,只能回头过去看B站视频补一补,因为语言都是相通的,而且一些细节都可以去查表解决,所以呢,我们只需要知道一些python与其他语言的不同和常见的优 ...