Description

题库链接

给出一棵 \(n\) 个节点的树, \(q\) 次询问,每次给出 \(k\) 个关键点。树上所有的点会被最靠近的关键点管辖,若距离相等则选编号最小的那个。求每个关键点管辖多少个节点。

\(1\leq n,q,\sum k\leq 300000\)

Solution

构出虚树后,我们能用简单的树形 \(dp\) 求出每个点离他最近的关键点。大体是做两遍 \(dfs\) 。第一遍用儿子更新父亲,第二遍用父亲更新儿子。

处理好这个之后,对于虚树上每个点。他的子树有两种:一个是虚树里的,一个是不在虚树里的。不在虚树里的后代肯定和他共用同一个关键点;

对于虚树上的一条边 \((u,v)\) ,我们需要找到 \((u,v)\) 边上的所有点以及他们连出去的块的最近点,更新答案。

如果 \(u, v\) 的最近点相同,那么这条边所代表的所有点的最近点肯定就是 \(u,v\) 的最近点;否则,可以用倍增找到临界点,计算贡献。

Code

//It is made by Awson on 2018.2.21
#include <bits/stdc++.h>
#define LL long long
#define dob complex<double>
#define Abs(a) ((a) < 0 ? (-(a)) : (a))
#define Max(a, b) ((a) > (b) ? (a) : (b))
#define Min(a, b) ((a) < (b) ? (a) : (b))
#define Swap(a, b) ((a) ^= (b), (b) ^= (a), (a) ^= (b))
#define writeln(x) (write(x), putchar('\n'))
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int N = 300000;
const int INF = ~0u>>1;
void read(int &x) {
char ch; bool flag = 0;
for (ch = getchar(); !isdigit(ch) && ((flag |= (ch == '-')) || 1); ch = getchar());
for (x = 0; isdigit(ch); x = (x<<1)+(x<<3)+ch-48, ch = getchar());
x *= 1-2*flag;
}
void print(int x) {if (x > 9) print(x/10); putchar(x%10+48); }
void write(int x) {if (x < 0) putchar('-'); print(Abs(x)); } int n, lim, u, v, fa[N+5][20], dep[N+5], size[N+5], dfn[N+5], times;
int flag[N+5], lst[N+5], kp[N+5], belong[N+5], dist[N+5], k, q, ans[N+5], S[N+5], top;
struct graph {
struct tt {int to, next; }edge[(N<<1)+5];
int path[N+5], top;
void add(int u, int v) {edge[++top].to = v, edge[top].next = path[u], path[u] = top; }
void dfs1(int o, int depth) {
dep[o] = depth, size[o] = 1, dfn[o] = ++times; for (int i = 1; i <= lim; i++) fa[o][i] = fa[fa[o][i-1]][i-1];
for (int i = path[o]; i; i = edge[i].next)
if (dfn[edge[i].to] == 0) fa[edge[i].to][0] = o, dfs1(edge[i].to, depth+1), size[o] += size[edge[i].to];
}
void dfs2(int o) {
belong[o] = 0, dist[o] = INF>>1;
if (flag[o]) belong[o] = o, dist[o] = 0;
for (int i = path[o]; i; i = edge[i].next) {
dfs2(edge[i].to);
if ((dist[edge[i].to]+dep[edge[i].to]-dep[o] < dist[o]) || (dist[edge[i].to]+dep[edge[i].to]-dep[o] == dist[o] && belong[o] > belong[edge[i].to])) dist[o] = dist[edge[i].to]+dep[edge[i].to]-dep[o], belong[o] = belong[edge[i].to];
}
}
void dfs3(int o) {
for (int i = path[o]; i; i = edge[i].next) {
if ((dist[o]+dep[edge[i].to]-dep[o] < dist[edge[i].to]) || (dist[o]+dep[edge[i].to]-dep[o] == dist[edge[i].to] && belong[o] < belong[edge[i].to])) dist[edge[i].to] = dist[o]+dep[edge[i].to]-dep[o], belong[edge[i].to] = belong[o];
dfs3(edge[i].to);
}
}
void dfs4(int o) {
int rem = size[o];
for (int &i = path[o]; i; i = edge[i].next) {
int x = edge[i].to; for (int j = lim; j >= 0; j--) if (dep[fa[x][j]] > dep[o]) x = fa[x][j];
rem -= size[x];
if (belong[edge[i].to] == belong[o]) ans[belong[o]] += size[x]-size[edge[i].to];
else {
int v = edge[i].to;
for (int j = lim; j >= 0; j--)
if (dep[fa[v][j]] >= dep[o])
if ((dist[edge[i].to]+dep[edge[i].to]-dep[fa[v][j]] < dist[o]+dep[fa[v][j]]-dep[o]) || (dist[edge[i].to]+dep[edge[i].to]-dep[fa[v][j]] == dist[o]+dep[fa[v][j]]-dep[o] && belong[edge[i].to] < belong[o])) v = fa[v][j];
ans[belong[o]] += size[x]-size[v];
ans[belong[edge[i].to]] += size[v]-size[edge[i].to];
}
dfs4(edge[i].to);
}
ans[belong[o]] += rem;
}
}g1, g2;
bool comp(const int &a, const int &b) {return dfn[a] < dfn[b]; }
int get_lca(int u, int v) {
if (dep[u] < dep[v]) Swap(u, v);
for (int i = lim; i >= 0; i--) if (dep[fa[u][i]] >= dep[v]) u = fa[u][i];
if (u == v) return u;
for (int i = lim; i >= 0; i--) if (fa[u][i] != fa[v][i]) u = fa[u][i], v = fa[v][i];
return fa[u][0];
} void work() {
read(n); lim = log(n)/log(2);
for (int i = 1; i < n; i++) read(u), read(v), g1.add(u, v), g1.add(v, u);
g1.dfs1(1, 1); read(q);
while (q--) {
read(k); for (int i = 1; i <= k; i++) read(lst[i]), flag[kp[i] = lst[i]] = 1;
top = g2.top = 0; sort(lst+1, lst+1+k, comp);
S[++top] = 1;
for (int i = 1; i <= k; i++) {
int lca = get_lca(lst[i], S[top]);
while (dfn[lca] < dfn[S[top]]) {
if (dfn[lca] >= dfn[S[top-1]]) {
g2.add(lca, S[top]), --top;
if (S[top] != lca) S[++top] = lca;
break;
}
g2.add(S[top-1], S[top]), --top;
}
if (S[top] != lst[i]) S[++top] = lst[i];
}
while (top > 1) g2.add(S[top-1], S[top]), --top;
g2.dfs2(1), g2.dfs3(1), g2.dfs4(1);
for (int i = 1; i < k; i++) write(ans[kp[i]]), putchar(' ');
writeln(ans[kp[k]]);
for (int i = 1; i <= k; i++) flag[kp[i]] = ans[kp[i]] = 0;
}
}
int main() {
work(); return 0;
}

[HNOI 2014]世界树的更多相关文章

  1. 数据结构(虚树,动态规划):HNOI 2014 世界树

    Hnoi2014 世界树 Description 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森,在他们的信条里,公平 ...

  2. 解题:HNOI 2014 世界树

    题面 首先建虚树 DFS求虚树上每个点所属的点和到它所属点的距离,然后在=考虑虚树所有的边(对应原树一条链).如果两个端点所属节点不同就倍增出分界点统计答案,否则不用管(之后会统计到的):注意根节点特 ...

  3. HNOI 2014 米特运输(图论)

    HNOI 2014 米特运输 题目大意 给一棵树,每个点有自己的权值,要求更改一些点的权值,使得整棵树满足两个条件: 同一个父亲的所有子节点权值相同 父节点的取值为所有子节点的和 答案输出最少要更改的 ...

  4. HNOI 2014

    D1T1:画框 frame 题意:给你两个n阶正整数方阵,请你求最大的\( \sum_{i = 1}^{n} A_{i, p_i}\times \sum_{i = 1}^{n} B_{i, p_i}  ...

  5. 图论(KM算法,脑洞题):HNOI 2014 画框(frame)

    aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAABPoAAANFCAIAAABtIwXVAAAgAElEQVR4nOydeVxTV/r/n9ertaJEC4

  6. [HNOI 2014]画框

    Description 题库链接 \(T\) 组询问,每组询问给你个 \(2\times N\) 的带权二分图,两个权值 \(a,b\) ,让你做匹配使得 \[\sum a\times \sum b\ ...

  7. [HNOI 2014]江南乐

    Description 题库链接 给你指定一个数 \(f\) ,并给你 \(T\) 组游戏,每组有 \(n\) 堆石子, \(A,B\) 两人轮流对石子进行操作,每次你可以选择其中任意一堆数量不小于 ...

  8. [HNOI 2014]道路堵塞

    Description A国有N座城市,依次标为1到N.同时,在这N座城市间有M条单向道路,每条道路的长度是一个正整数.现在,A国 交通部指定了一条从城市1到城市N的路径,并且保证这条路径的长度是所有 ...

  9. [HNOI 2014]米特运输

    Description 米特是D星球上一种非常神秘的物质,蕴含着巨大的能量.在以米特为主要能源的D星上,这种米特能源的运输和储 存一直是一个大问题.D星上有N个城市,我们将其顺序编号为1到N,1号城市 ...

随机推荐

  1. 云+社区分享——腾讯云OCR文字识别

    欢迎大家前往腾讯云+社区,获取更多腾讯海量技术实践干货哦~ 本文由云+社区运营团队发布在腾讯云+社区 前言 2018年3月27日腾讯云云+社区联合腾讯云智能图像团队共同在客户群举办了腾讯云OCR文字识 ...

  2. JavaEE Servlet 核心方法及生命周期

    做JavaWeb开发,免不了要和Servlet打交道.Servlet是Sun(Oracle)官方定义的一个Web开发规范,所有Servlet开发都必须遵守.自己以前也没有从头做过Web开发,所以这方面 ...

  3. 通过cmd窗口导入导出mysql数据库

    1.导入数据库 使用source命令 首先要在cmd窗口中连接数据库,然后再用source命令进行导入操作 mysql>use 数据库名 mysql>source d:/dbname.sq ...

  4. Beta冲刺合集

    Beta冲刺序列: Beta凡事预则立 :Beta No.0 Beta冲刺Day1:Beta No.1 Beta冲刺Day2:Beta No.2 Beta冲刺Day3:Beta No.3 Beta冲刺 ...

  5. 实验四Java Android简易开发

    实验准备 Android Studio下载 Android Studio安装 实验内容 Android Stuidio的安装测试 Android Stuidio的安装测试: 参考<Java和An ...

  6. 201621123057 《Java程序设计》第4周学习总结

    1. 本周学习总结 1.1 写出你认为本周学习中比较重要的知识点关键词 答: (普通方法 / 构造函数)重载. static . final.继承与多态.extends.object类.abstrac ...

  7. 2017-2018-1 我爱学Java 第八周 作业

    团队六七周作业 团队分工 UML图 工具选择 小编(金立清)有话说 参考资料 团队分工 返回目录 UML图 用例图 类图 活动图 状态图 返回目录 工具选择 ProcessOn - 免费在线作图,实时 ...

  8. 位图(BitMap)

    位图 引入:2.给40亿个不重复的无符号整数,没排过序.给一个无符号整数,如何快速判断一个数是否在这40亿个数中. [腾讯] 思路:如果内存够的话,40亿个整型使用位图存储需要500M左右的空间. 分 ...

  9. 构建微服务开发环境3————Java应用的优秀管理工具Maven的下载安装及配置

    [内容指引] 下载安装包: MacOS下Maven的安装及配置: Windows下Maven的安装及配置. 一.下载安装包 进入Maven的官方下载地址:http://maven.apache.org ...

  10. c# 运算符:? ,??

    参考微软帮助 1  ?  空值条件运算符,用于在执行成员访问 (?.) 或索引 (?[) 操作之前,测试是否存在 NULL. // ? 空值条件运算符 string str = null; Conso ...