「2017 山东三轮集训 Day7」Easy
一棵带边权的树,多次询问 $x$ 到编号为 $[l,r]$ 的点最短距离是多少
$n \leq 100000$
sol:
动态点分治,每层重心维护到所有点的距离
查询的时候在管辖这个点的 log 层线段树里查就可以了
因为这样每一层的答案只会漏而不会错,所以正确性有保障
不会写点分治了...orz
#include <bits/stdc++.h>
#define LL long long
#define rep(i, s, t) for (register int i = (s), i##end = (t); i <= i##end; ++i)
#define dwn(i, s, t) for (register int i = (s), i##end = (t); i >= i##end; --i)
using namespace std;
inline int read() {
int x = , f = ;
char ch;
for (ch = getchar(); !isdigit(ch); ch = getchar())
if (ch == '-')
f = -f;
for (; isdigit(ch); ch = getchar()) x = * x + ch - '';
return x * f;
}
const int maxn = , inf = 1e9;
int n, m;
int first[maxn], to[maxn << ], nx[maxn << ], val[maxn], cnt;
inline void add(int u, int v, int w) {
to[++cnt] = v;
val[cnt] = w;
nx[cnt] = first[u];
first[u] = cnt;
}
int sig, root, size[maxn], pa[maxn], f[maxn], vis[maxn];
inline void findroot(int x, int fa) {
size[x] = ;
f[x] = ;
for (int i = first[x]; i; i = nx[i]) {
if ((to[i] == fa) || (vis[to[i]]))
continue;
findroot(to[i], x);
size[x] += size[to[i]];
f[x] = max(f[x], size[to[i]]);
}
f[x] = max(f[x], sig - size[x]);
if (f[x] < f[root])
root = x;
}
inline void build(int x) {
vis[x] = ;
for (int i = first[x]; i; i = nx[i]) {
if (vis[to[i]])
continue;
root = ;
sig = size[to[i]];
findroot(to[i], x);
pa[root] = x;
build(root);
}
}
namespace sp_lca {
int size[maxn], dep[maxn], fa[maxn], bl[maxn], dp[maxn];
inline void dfs1(int x) {
size[x] = ;
for (int i = first[x]; i; i = nx[i]) {
if (to[i] == fa[x])
continue;
fa[to[i]] = x;
dep[to[i]] = dep[x] + ;
dp[to[i]] = dp[x] + val[i];
dfs1(to[i]);
size[x] += size[to[i]];
}
}
inline void dfs2(int x, int col) {
int k = ;
bl[x] = col;
for (int i = first[x]; i; i = nx[i])
if (to[i] != fa[x] && size[to[i]] > size[k])
k = to[i];
if (!k)
return;
dfs2(k, col);
for (int i = first[x]; i; i = nx[i])
if (to[i] != fa[x] && to[i] != k)
dfs2(to[i], to[i]);
}
inline void init() {
dfs1();
dfs2(, );
}
inline int lca(int x, int y) {
while (bl[x] != bl[y]) {
if (dep[bl[x]] < dep[bl[y]])
swap(x, y);
x = fa[bl[x]];
}
return dep[x] < dep[y] ? x : y;
}
} // namespace sp_lca
inline int dis(int x, int y) { return sp_lca::dp[x] + sp_lca::dp[y] - (sp_lca::dp[sp_lca::lca(x, y)] << ); }
int rt[maxn], ls[maxn << ], rs[maxn << ], va[maxn << ], ToT;
inline void Insert(int &x, int l, int r, int pos, int v) {
if (!x)
x = ++ToT, va[x] = inf;
va[x] = min(va[x], v);
if (l == r)
return;
int mid = (l + r) >> ;
if (pos <= mid)
Insert(ls[x], l, mid, pos, v);
else
Insert(rs[x], mid + , r, pos, v);
}
inline int query(int x, int l, int r, int L, int R) {
if (!x)
return inf;
if (L <= l && r <= R)
return va[x];
int mid = (l + r) >> , res = inf;
if (L <= mid)
res = min(res, query(ls[x], l, mid, L, R));
if (R > mid)
res = min(res, query(rs[x], mid + , r, L, R));
return res;
}
int main() {
// freopen("a.in","r",stdin);
// freopen("a.out","w",stdout);
n = read();
rep(i, , n) {
int u = read(), v = read(), w = read();
add(u, v, w);
add(v, u, w);
}
f[] = inf;
root = ;
sig = n;
findroot(, );
build(root);
sp_lca ::init();
rep(i, , n) for (int x = i; x; x = pa[x]) Insert(rt[x], , n, i, dis(x, i));
// rep(i, 1, n) cout << pa[i] << endl;
int q = read();
while (q--) {
int u = read(), v = read(), x = read();
int cur = x, ans = inf;
while (cur) {
ans = min(ans, query(rt[cur], , n, u, v) + dis(cur, x));
cur = pa[cur];
}
cout << ans << endl;
}
}
「2017 山东三轮集训 Day7」Easy的更多相关文章
- 【loj6145】「2017 山东三轮集训 Day7」Easy 动态点分治+线段树
题目描述 给你一棵 $n$ 个点的树,边有边权.$m$ 次询问,每次给出 $l$ .$r$ .$x$ ,求 $\text{Min}_{i=l}^r\text{dis}(i,x)$ . $n,m\le ...
- #6145. 「2017 山东三轮集训 Day7」Easy 动态点分治
\(\color{#0066ff}{题目描述}\) JOHNKRAM 最近在参加 C_SUNSHINE 举办的聚会. C 国一共有 n 座城市,这些城市由 n−1 条无向道路连接.任意两座城市之间有且 ...
- LOJ #6145. 「2017 山东三轮集训 Day7」Easy 点分树+线段树
这个就比较简单了~ Code: #include <cstdio> #include <algorithm> #define N 100004 #define inf 1000 ...
- 「2017 山东三轮集训 Day7 解题报告
「2017 山东三轮集训 Day7」Easy 练习一下动态点分 每个点开一个线段树维护子树到它的距离 然后随便查询一下就可以了 注意线段树开大点... Code: #include <cstdi ...
- 「2017 山东三轮集训 Day1」Flair
模拟赛的题 好神仙啊 题面在这里 之前的Solution很蠢 现在已经update.... 题意 有$ n$个商品价格均为$ 1$,您有$ m$种面值的货币,面值为$ C_1..C_m$ 每种物品你有 ...
- 【loj6142】「2017 山东三轮集训 Day6」A 结论题+Lucas定理
题解: 当奇数 发现答案就是C(n,1)^2+C(n,3)^2+...C(n,n)^2 倒序相加,发现就是C(2n,n) 所以答案就是C(2n,n)/2 当偶数 好像并不会证 打表出来可以得到 2.当 ...
- loj #6138. 「2017 山东三轮集训 Day4」Right
题目: 题解: 暴力一波 \(SG\) 函数可以发现这么一个规律: \(p\) 为奇数的时候 : \(SG(n) = n \% 2\) \(p\) 为偶数的时候 : \(SG(n) = n \% (p ...
- loj #6136. 「2017 山东三轮集训 Day4」Left
题目: 题解: 我们可以发现所有的交换器都是一个位置连接着下一层左侧的排序网络,另一个位置连着另一侧的排序网络. 而下一层是由两个更低阶的排序网络构成的. 两个网络互不干扰.所以我们可以通过第一行和最 ...
- Loj #6142. 「2017 山东三轮集训 Day6」A
link: https://loj.ac/problem/6142 推完一波式子之后发现求的是:ΣC(N,i)^2, 其中i是偶数. 然后就可以卢卡斯乱搞了,分奇偶和之前的答案合并就好了233. #i ...
随机推荐
- 内置函数(Day16)
现在python一共为我们提供了68个内置函数.它们就是python提供给你直接可以拿来使用的所有函数 内置函数 abs() divmod() input() open() stati ...
- VMware虚拟机NAT模式的具体配置
NAT模式的具体配置 NAT方式:虚拟机可以上外网,可以访问宿主计算机所在网络的其他计算机(反之不行). 1.1.1. 查看虚拟机的网络参数 1) 打开虚拟机,选择菜单&q ...
- Video标签的使用
现在如果要在页面中使用video标签,需要考虑三种情况,支持Ogg Theora或者VP8(如果这玩意儿没出事的话)的(Opera.Mozilla.Chrome),支持H.264的(Safari.IE ...
- Python3 optparse模块
Python 有两个内建的模块用于处理命令行参数: 一个是 getopt,<Deep in python>一书中也有提到,只能简单处理 命令行参数: 另一个是 optparse,它功能强大 ...
- ASP.NET 4.0 ListView等容器控件中获取ClientID值与HTML中自动生成ID字符串不一样问题。
ASP.NET 4.0 中 ClientIDMode的属性 可以设置获取不同ID格式的值. 项目中遇到的问题: 1.ListView1 ItemDataBound事件中,获取ClientID结果与自动 ...
- Springboot WebSocket例子
Springboot整合WebSocket 1.application.properties #设置服务端口号 server.port=8080 #thymeleaf配置 #是否启用模板缓存. spr ...
- 【leetcode刷提笔记】Search Insert Position
Given a sorted array and a target value, return the index if the target is found. If not, return the ...
- Django用户注册、邮箱验证实践
算法流程如下:1)处理用户注册数据,存入数据库,is_activity字段设置为False,用户未认证之前不允许登陆2)产生token,生成验证连接URL3)发送验证邮件4)用户通过认证邮箱点击验证连 ...
- INSPIRED启示录 读书笔记 - 第22章 原型测试
物色测试者 1.如果你已经拥有一批特约用户,可以邀请他们参加测试 2.如果是企业级产品,同类产品的展销会是寻找目标用户的好去处 3.可以在分类信息网站上发布广告,征集测试者.征集要求可以写得笼统些,不 ...
- Go Log模块生成日志文件
使用log模块示例代码: package main import ( "fmt" "time" "log" "os" ) ...