[题解] [HNOI2014] 世界树
题面
题解
从数据范围很容易看出是个虚树DP(可惜看出来了也还是不会做)
虚树大家应该都会, 不会的话自己去搜吧, 我懒得讲了, 我们在这里只需要考虑如何DP即可
首先我们需要求出每个点被哪个点所控制, 设\(u\)点被\(bl[u]\)所控制, 两遍DFS即可, 考虑儿子对父亲的影响和父亲对儿子的影响
代码细节相信不要我说, 能做这道题的总不可能不会DFS吧
还是贴一下自己看吧
void dfs1(int u, int fa)
{
for(int i = head[u]; i; i = e[i].next)
{
int v = e[i].to; if(v == fa) continue;
dfs1(v, u); int dv = dep[bl[v]] - dep[u], du = bl[u] ? dep[bl[u]] - dep[u] : 0x3f3f3f3f;
if(dv < du || (dv == du && bl[v] < bl[u])) bl[u] = bl[v];
}
}
void dfs2(int u, int fa)
{
for(int i = head[u]; i; i = e[i].next)
{
int v = e[i].to; if(v == fa) continue;
int dv = dis(bl[v], v), du = dis(bl[u], v);
if(du < dv || (du == dv && bl[u] < bl[v])) bl[v] = bl[u];
dfs2(v, u);
}
}
//分别在递归前和递归后处理一下就完事了
然后考虑如何计算答案, 对虚树上每一条边讨论
①: 边的两端被同一个节点所控制, 加上这两个点不在虚树中的儿子的sz即可
②: 边的两端被不同点控制, 我们需要找出一个分界点, 满足此分界点归下面那个点的\(bl[]\)所控制, 此分界点的父亲被上面那个点的\(bl[]\)所控制, 用个数据结构维护一下或者倍增跳一下就可以了
至于每个点不在虚树中的儿子的sz, 拿当前点的sz减去所有他在虚树中的儿子的sz即可
具体实现细节参见代码(参考了一下题解的思路嘿嘿嘿)
代码实现
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#define N 300005
using namespace std;
int n, m, Q, head[N], sz[N], son[N], dep[N], dfn[N], f[N][21], top[N], cnt, tp, a[N], b[N], stk[N], bl[N], con[N], ans[N], l[N];
struct edge { int from, to, next; } e[N << 1];
inline int read()
{
int x = 0, w = 1;
char c = getchar();
while(c < '0' || c > '9') { if (c == '-') w = -1; c = getchar(); }
while(c >= '0' && c <= '9') { x = x * 10 + c - '0'; c = getchar(); }
return x * w;
}
inline void add(int u, int v) { e[++cnt] = (edge) { u, v, head[u] }; head[u] = cnt; }
void dfs_sz(int u, int fa)
{
sz[u] = 1; dep[u] = dep[fa] + 1; f[u][0] = fa;
for(int i = head[u]; i; i = e[i].next)
{
int v = e[i].to; if(v == fa) continue;
dfs_sz(v, u); sz[u] += sz[v]; if(sz[son[u]] < sz[v]) son[u] = v;
}
}
void dfs_top(int x, int y)
{
dfn[x] = ++cnt; top[x] = y;
if(!son[x]) return; dfs_top(son[x], y);
for(int i = head[x]; i; i = e[i].next) if(e[i].to != f[x][0] && e[i].to != son[x]) dfs_top(e[i].to, e[i].to);
}
int LCA(int x, int y)
{
while(top[x] != top[y])
{
if(dep[top[x]] < dep[top[y]]) swap(x, y);
x = f[top[x]][0];
}
return dep[x] < dep[y] ? x : y;
}
bool cmp(int x, int y) { return dfn[x] < dfn[y]; }
void dfs1(int u, int fa)
{
for(int i = head[u]; i; i = e[i].next)
{
int v = e[i].to; if(v == fa) continue;
dfs1(v, u); int dv = dep[bl[v]] - dep[u], du = bl[u] ? dep[bl[u]] - dep[u] : 0x3f3f3f3f;
if(dv < du || (dv == du && bl[v] < bl[u])) bl[u] = bl[v];
}
}
int dis(int x, int y) { return dep[x] + dep[y] - 2 * dep[LCA(x, y)]; }
void dfs2(int u, int fa)
{
for(int i = head[u]; i; i = e[i].next)
{
int v = e[i].to; if(v == fa) continue;
int dv = dis(bl[v], v), du = dis(bl[u], v);
if(du < dv || (du == dv && bl[u] < bl[v])) bl[v] = bl[u];
dfs2(v, u);
}
}
void dp(int u)
{
for(int s, mid, nt, dv, du, i = head[u]; i; i = e[i].next)
{
int v = e[i].to; dp(v); s = mid = v;
for(int j = l[dep[v]]; j >= 0; j--) if(dep[f[s][j]] > dep[u]) s = f[s][j];
con[u] -= sz[s];
if(bl[u] == bl[v]) { ans[bl[u]] += sz[s] - sz[v]; continue; }
for(int j = l[dep[v]]; j >= 0; j--)
{
nt = f[mid][j]; if(dep[nt] <= dep[u]) continue;
dv = dis(bl[v], nt), du = dis(bl[u], nt);
if(dv < du || (dv == du && bl[v] < bl[u])) mid = nt;
}
ans[bl[u]] += sz[s] - sz[mid];
ans[bl[v]] += sz[mid] - sz[v];
}
ans[bl[u]] += con[u];
}
void query()
{
m = read(); cnt = tp = 0;
for(int i = 1; i <= m; i++) bl[a[++cnt] = b[i] = read()] = b[i];
sort(a + 1, a + cnt + 1, cmp);
for(int i = 1; i < m; i++) a[++cnt] = LCA(a[i], a[i + 1]);
a[++cnt] = 1; sort(a + 1, a + cnt + 1, cmp);
int len = unique(a + 1, a + cnt + 1) - a - 1;
cnt = 0; for(int i = 1; i <= len; i++) head[a[i]] = 0, con[a[i]] = sz[a[i]];
for(int i = 1; i <= len; i++)
{
while(tp && dfn[a[i]] >= dfn[stk[tp]] + sz[stk[tp]]) tp--;
if(tp) add(stk[tp], a[i]); stk[++tp] = a[i];
}
dfs1(1, 0); dfs2(1, 0); dp(1);
for(int i = 1; i <= m; i++) printf("%d%c", ans[b[i]], i == m ? '\n' : ' ');
for(int i = 1; i <= len; i++) bl[a[i]] = ans[a[i]] = con[a[i]] = 0;
}
int main()
{
n = read(); for(int i = 2; i <= n; i++) l[i] = l[i >> 1] + 1;
for(int i = 1; i < n; i++) { int u = read(), v = read(); add(u, v); add(v, u); }
cnt = 0; dfs_sz(1, 0); dfs_top(1, 0);
for(int i = 1; i <= n; i++)
for(int j = 1; j <= 20; j++)
f[i][j] = f[f[i][j - 1]][j - 1];
Q = read(); while(Q--) query();
return 0;
}
[题解] [HNOI2014] 世界树的更多相关文章
- [BZOJ3572][Hnoi2014]世界树
[BZOJ3572][Hnoi2014]世界树 试题描述 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森,在他们的信条 ...
- 【BZOJ3572】[Hnoi2014]世界树 虚树
[BZOJ3572][Hnoi2014]世界树 Description 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森 ...
- bzoj3572[Hnoi2014] 世界树 虚树+dp+倍增
[Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1921 Solved: 1019[Submit][Status][Dis ...
- BZOJ 3572: [Hnoi2014]世界树
BZOJ 3572: [Hnoi2014]世界树 标签(空格分隔): OI-BZOJ OI-虚数 OI-树形dp OI-倍增 Time Limit: 20 Sec Memory Limit: 512 ...
- bzoj 3572: [Hnoi2014]世界树 虚树 && AC500
3572: [Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 520 Solved: 300[Submit][Status] ...
- bzoj 3572 [Hnoi2014]世界树(虚树+DP)
3572: [Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 645 Solved: 362[Submit][Status] ...
- <虚树+树型DP> HNOI2014世界树
<虚树+树型DP> HNOI2014世界树 #include <iostream> #include <cstdio> #include <cstring&g ...
- BZOJ3572:[HNOI2014]世界树——题解
+++++++++++++++++++++++++++++++++++++++++++ +本文作者:luyouqi233. + +欢迎访问我的博客:http://www.cnblogs.com/luy ...
- 【题解】HNOI2014世界树
脑子不清醒的时候千万别写题.写题写不下去了千万别死扛,重构才是你唯一的出路QAQ 昨天很想快点写道题,思路没有很清晰的时候就写了,结果……今天一怒之下决定重整思路重构代码,其实不过是半个小时的事情…… ...
随机推荐
- c# mvc权限系统
登陆界面可以用index1,现在这个是仿照这个写的 登陆后的界面.模板都是套的,数据那个easyui自己写的 后台直接三层架构,dapper连处理数据后面加了autofac demo地址:https: ...
- 你不知道的css各类布局(三)之自适应布局
自适应布局 概念 自适应布局(Adaptive Layout)是对凡是有自适应特性的一类布局的统称 自适应布局使用media query来检测当前浏览器的宽度进而通过CSS样式调整页面大小.自适应布局 ...
- JavaScript函数式编程——柯里化
柯里化原理 如何实现柯里化 柯里化的应用 一.柯里化原理 柯里化:在数学和计算机科学中,柯里化是一种使用多个参数的一个函数转换成一系列使用一个参数的函数的技术. 前端使用柯里化的用途主要就应该是简化代 ...
- linux环境下,springboot jar启动方式
linux环境下,springboot jar启动方式 一.前台启动(ctrl+c会关掉进程) java -jar application.jar 二.后台启动(ctrl+c不会关闭) java -j ...
- TVM使用问题记录
1.numpy提示repeat错误 错误信息为 One method of fixing this is to repeatedly uninstall numpy until none is fou ...
- ONNX源码安装
ONNX是facebook提出的一个 Open Neural Network Exchange协议,能够让训练好的模型在不同的框架间进行交互. ONNX的安装相对来说不是特别麻烦,麻烦的是其依赖库的安 ...
- 【异常】Phoenix异常:java.lang.ArithmeticException: Rounding necessary
1 异常sql upsert into WMBIGDATA.ODS_ES_CHARGING_STATION(id,evcosType,address,serviceTel,supportOrder,o ...
- 利用协程和socket实现并发
服务端代码 from gevent import monkey monkey.patch_all() from gevent import spawn import socket def commun ...
- Centos下安装compoer, 解决vendor目录问题
laravel部署时使用git拉代码的话,vendor目录部署是没有提交到git上的,如果没有该目录,服务会报错如下 Warning: require(/opt/soft/pro/bootstrap/ ...
- 复杂sql优化步骤与技巧
数据管理型系统,由于用户的要求或者系统设计要求,会出现大量表进行join,还要进行大量统计性数据查询展示,甚至数据权限控制等操作.最后会导致sql异常复杂,随着数据量增加,或者只是应用到生产环境(正式 ...