2020-9-29 T3
题意:给定一颗大小为 \(n(n \le 5 \times 10 ^ 4)\) 的树,保证树的生成方式随机,你需要选定两个点 \(x, y\),最小化:
\]
你会发现这个问题很类似 [APIO2015]八邻旁之桥,那么同样的我们可以先思考只有单点的情况。
即考虑如何选定 \(x\) 最小化:
\]
那么是不是只要能快速算出每个点 \(x\) 的 \(w_x = \sum\limits_{i = 1} ^ n a_i \times dis_{x, i}\) 即可呢?
来尝试一下,因为树高是 \(\log n\) 的,因此我们要尽可能地使用和树高有关的做法。
那么根据两点之间的距离公式 \(dis_{x, i} = dep_{x} + dep_i - 2 \times dep_{lca}\),显然 \(dep_x\) 是固定的。
其中最难处理的就是 \(dep_{lca}\) 了,因此我们要想办法固定住 \(dep_{lca}\)。
结合之前树高是 \(\log n\) 的性质,不难发现我们可以一步步向上跳,将每次跳到的点记为 \(y\)。
如果此时 \(y\) 是 \(lca_{x, i}\),显然 \(i\) 就应该是除了 \(x\) 所在子树以外的所有点。
那么那么此时这些点 \(i\) 到 \(x\) 的距离就是其到 \(y\) 的距离在加上 \(y\) 到 \(x\) 的距离。
这部分多与的距离很好维护,直接令 \(S_x = \sum\limits_{i \in x} a_i\),然后减一下乘距离即可。
那么这些点到 \(y\) 的距离怎么维护呢,不难发现实际上这是 \(y\) 的整颗子树到其距离减去 \(x\) 所在子树到 \(y\) 的距离即可。
那么就只需要维护 \(f_x = \sum\limits_{i \in x} a_i \times dis_{x, i}\) 以及 \(g_x = \sum\limits_{i \in x} a_i \times dis_{i, fa_x}\),即可。
上面这两个值通过子树合并即可计算。
通过最开始的流程,不难发现 \(w_x\) 的计算方法。
值得一提的是,上面这个做法本质上使用的是点分治的思想。
那么现在一维的情况已经解决,可以来考虑二维的情况。
和最开始提到的那道题很类似的是,对于选定的两个点 \(x, y\) 一定是从这条链的中点分开一部分全部走去 \(x\),另一部分全部走去 \(y\)。
那么我们可以类似与那道题的做法,断去树中的一条边,两边的树就只需要跑上面哪个一维的子问题了。
但是你会发现的是,如果我们继续重新对两边重新计算 \(w_x\),复杂度还是 \(O(n ^ 2)\) 的,同样的我们还是希望能找到一个和树高有关的做法。
可以继续沿用上面递推的思想,看一看能否先计算出一个点的 \(w_x\) 然后沿着树边往下递推。
可以考虑一个简单的情况,如果 \(x\) 往其一个儿子 \(y\) 走 \(w_y\) 的变化情况。
不难发现会有:\(w_y = w_x + S - 2 \times S_y, S\) 为整颗树的权值和。
那么如果 \(w_y\) 会变小当且仅当 \(S - 2 \times S_y < 0\),移项可得:\(S_y > \frac{S}{2}\)。
可以惊奇地发现对于显然这样的 \(y\) 只会存在一个。
那么如果存在这样的点 \(y\) 就一定可以从当前的根已知往下走直到 \(S - 2 \times S_y > 0\) 为止。
可以发现上述这个过程在数随机生成的情况下是可以保证为 \(O(\log n)\) 的。
那么现在的任务就在于如何计算出一个断边 \((u, v)(fa_v = u)\) 后的 \(w_u, w_v\) 即可。
不难发现有 \(w_v = f_v, w_u = w_u - g_v\)。
并且可以发现断边以后将 \(u\) 做为新树的根只会影响 \(u\) 到根路径上的点的 \(S\) 值,直接暴力修改即可。
于是我们就在 \(O(n \log n)\) 的复杂度下完成了本题,事实上从大佬哪里听说貌似使用动态点分治就可以不依赖于树随机生成的性质,以后学习了动态点分治再说吧。
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define rep(i, l, r) for (int i = l; i <= r; ++i)
#define Next(i, u) for (int i = h[u]; i; i = e[i].next)
const int N = 50000 + 5;
const int inf = 1e18;
struct edge {
int v, next;
}e[N << 1];
int n, u, v, tot, top, ans = inf, cnt1, cnt2, h[N], a[N], f[N], g[N], w[N], S[N], fa[N], st[N], dep[N];
int read() {
char c; int x = 0, f = 1;
c = getchar();
while (c > '9' || c < '0') { if(c == '-') f = -1; c = getchar();}
while (c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
void add(int u, int v) {
e[++tot].v = v, e[tot].next = h[u], h[u] = tot;
e[++tot].v = u, e[tot].next = h[v], h[v] = tot;
}
void dfs(int u, int fa, int topf) {
f[topf] += (dep[u] - dep[topf]) * a[u];
Next(i, u) if(e[i].v != fa) dfs(e[i].v, u, topf);
}
void Prefix(int u, int Fa) {
S[u] = a[u], dep[u] = dep[Fa] + 1, fa[u] = Fa;
Next(i, u) if(e[i].v != Fa) Prefix(e[i].v, u), S[u] += S[e[i].v];
dfs(u, Fa, u), g[u] = f[u] + S[u], w[u] = f[u];
}
void access(int u, int fa, int &cnt, int tmp, int topf) {
cnt = tmp;
Next(i, u) {
int v = e[i].v; if(v == fa) continue;
if(S[topf] - 2 * S[v] < 0) access(v, u, cnt, tmp + S[topf] - 2 * S[v], topf);
}
}
void solve(int u, int Fa) {
Next(i, u) {
int v = e[i].v; if(v == Fa) continue;
st[u] = S[u], S[u] = S[1] - S[v];
for (int j = fa[u], k = u; j; k = j, j = fa[j]) st[j] = S[j], S[j] = S[1] - st[k];
access(v, u, cnt1, f[v], v), access(u, v, cnt2, w[u] - g[v], u);
for (int j = u; j; j = fa[j]) S[j] = st[j];
ans = min(ans, cnt1 + cnt2), solve(v, u);
}
}
signed main() {
n = read();
rep(i, 1, n - 1) u = read(), v = read(), add(u, v);
rep(i, 1, n) a[i] = read();
Prefix(1, 0);
rep(i, 1, n) for (int j = i; fa[j]; j = fa[j]) {
w[i] += f[fa[j]] - g[j] + (dep[i] - dep[fa[j]]) * (S[fa[j]] - S[j]);
}
solve(1, 0);
printf("%lld", ans);
return 0;
}
值得一提的是,当树高期望为 \(\log n\) 时,一个很好运用性质的方式就是从叶子开始往上递推或者从上往下递推。
2020-9-29 T3的更多相关文章
- OpenWrt(LEDE)2020.4.29更新 UPnP+NAS+多拨+网盘+DNS优化+帕斯沃 无缝集成+软件包
交流群:QQ 1030484865 电报: t_homelede 固件说明 基于Lede OpenWrt R2020.4.8版本(源码截止2020.4.29)Lienol Feed及若干自行维护 ...
- 学习 Gin 问题总结 2020.12.29
学习 Gin 问题总结 2020.12.29 数据绑定与解析 BindXXX,ShouldBindXXX和ShouldBindWith区别 BindXXX 会自动返回信息,输入无效时,在header写 ...
- 2020/1/29 PHP代码审计之进一步学习XSS【持续更新】
0x00 上午学习了XSS漏洞,中午吃饭想了想,还是思考的太浅层了,这种老生常谈的东西对于现在的我意义不大.现在我需要的是思考.于是就有了这个随笔.在本文中,我会持续更新一些XSS的深入思考,payl ...
- Test 6.29 T3 小学生
问题描述 "不错,不错!那么,准备好迎接下一道题了么?"一道白光闪过,CJK 眼前出现了 1e100 个小学生."他们中,有一些人轨了我的机子.现在,我需要你在 1S 之 ...
- 2020/1/29 PHP代码审计之XSS漏洞
0x00 XSS漏洞简介 人们经常将跨站脚本攻击(Cross Site Scripting)缩写为CSS,但这会与层叠样式表(Cascading Style Sheets,CSS)的缩写混淆.因此,有 ...
- 2020年AI、CV、NLP顶会最全时间表
2020年AI.CV.NLP顶会最全时间表 2019-09-01 14:04:19 weixin_38753768 阅读数 40 2020 AI.CV.NLP主流会议时间表,包含会议举办的时间.地 ...
- 软路由OpenWrt(LEDE)2020.5.10更新 UPnP+NAS+多拨+网盘+DNS优化
近期更新:2020.05.10更新-基于OpenWrt R2020.5.9版本,源码截止2020.05.10. 交流群:QQ 1030484865 电报 t.me/t_homelede 版本说 ...
- OpenWrt(LEDE)2020.4.12编译 UnPnP+NAS+多拨+网盘+DNS优化+帕斯沃 无缝集成
固件说明 基于Lede OpenWrt R2020.4.8版本(源码截止2020.4.12)Lienol Feed及若干自行维护的软件包 结合家庭x86软路由场景需要定制 按照家庭应用场景对固件及软件 ...
- HomeLede 2020.5.27更新 UPnP+NAS+多拨+网盘+DNS优化+帕斯沃/Clash 无缝集成+软件包
交流群:QQ 1030484865 电报 t.me/t_homelede 固件说明 基于Lede OpenWrt R2020.5.20版本(源码截止2020.5.27)及若干自行维护的软件包 结合 ...
- OpenGL学习日志(2020.4之前)
咳咳,原本这个日志是本机上随便写的一些记录,也没怎么注意可读性和格式,有用信息密度很小,所以实用价值并不大.暂时由于不可抗因素得先鸽一段落了... 后续的日志会升格为模块化的学习记录,(应该)将会有很 ...
随机推荐
- Dubbo 的设计思想
在java远程调用多年的沉淀 <1>首先是socket调用.在orderService中开放socket服务,在userService中进行远程调用. 优点:解决了单机调用的问题. 缺点: ...
- 洛谷——P1980 [NOIP2013 普及组] 计数问题
题目描述 试计算在区间 11 到 nn的所有整数中,数字x(0 ≤ x ≤ 9)x(0≤x≤9)共出现了多少次?例如,在 11到1111中,即在 1,2,3,4,5,6,7,8,9,10,111,2, ...
- Java语言程序设计复习提纲
这是我在准备Java考试时整理的提纲,如果是通过搜索引擎搜索到这篇博客的师弟师妹,建议还是先参照PPT和课本,这个大纲也不是很准确,自己总结会更有收获,多去理解含义,不要死记硬背,否则遇到概念辨析题 ...
- 使用所学Spring知识,实现简易的图书查询系统功能。实现查询全部图书。 根据书籍编号查询信息。 根据书名查询书籍信息。 根据状态查询书籍信息。
相关 知识 >>> 相关 练习 >>> 实现要求: 使用所学Spring知识,实现简易的图书查询系统功能. 要求如下: 查询全部图书. 根据书籍编号查询信息. 根据 ...
- iNeuOS工业互联网操作系统,增加DTU与平台实时交互的应用场景
目 录 1. 概述... 2 2. 平台演示... 2 3. 硬件设置... 2 4. 应用过程... 3 1. 概述 DTU向下连接硬件传感器 ...
- ProtoBuf3语法指南(Protocol Buffers)_上
0.说明 ProtoBuf3语法指南, 又称为proto3, 是谷歌的Protocol Buffers第3个版本. 本文基于官方英文版本翻译, 加上了自己的理解少量修改, 一共分为上下两部分. 1.序 ...
- openmesh - impl - Remove Duplicated Vertices
openmesh - impl - Remove Duplicated Vertices 关于openmesh元素删除实现的介绍参见:openmesh - src - trimesh delete a ...
- python 自动化测试框架unittest与pytest的区别
前言: 有使用过unittest单元测试框架,再使用pytest单元测试框架,就可以明显感觉到pytest比unittest真的简洁.方便很多. unittest与pytest的区别: 主要从用例编写 ...
- python 数据库查询,将查询的列名也保存下来,返回字典类型的数据
前言: python + pymysql:通过pymysql,查询到符合要求的数据,不做类型转化的,返回的时元组类型的数据,可以根据下标,取想要的值. self.cursor.execute(sql) ...
- Linux中常用的监控性能的命令(sar、mpstat,vmstat, iostat,)详解
Linux中常用的监控性能的命令有: sar:能查看CPU的平均信息,还能查看指定CPU的信息.与mpstat相比,sar能查看CPU历史信息 mpstat:能查看所有CPU的平均信息,还能查看指定C ...