【题解】HNOI2014世界树
脑子不清醒的时候千万别写题。写题写不下去了千万别死扛,重构才是你唯一的出路QAQ
昨天很想快点写道题,思路没有很清晰的时候就写了,结果……今天一怒之下决定重整思路重构代码,其实不过是半个小时的事情……
提示很明显,总点数限制了范围。建立出虚树,在虚树上面 dp。 在虚树上面我们dp两遍,两遍一起处理出每一个(虚树上的)节点被谁管辖,在连接两者的路径上找到管理的分界点即可。处理两点之间的距离一定要上 ST 表啊……
#include <bits/stdc++.h>
using namespace std;
#define maxn 350000
#define CNST 20
int n, dep[maxn], gra[maxn][CNST];
int size[maxn], hson[maxn], ans[maxn];
int top, S[maxn], a[maxn], b[maxn];
int timer, dfn[maxn], sum[maxn];
bool mark[maxn], vis[maxn];
int T, Log[maxn * 2], bit[21], ST[maxn * 2][CNST], pos[maxn]; int read()
{
int x = 0, k = 1;
char c;
c = getchar();
while(c < '0' || c > '9') { if(c == '-') k = -1; c = getchar(); }
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * k;
} struct edge
{
int cnp, to[maxn * 2], last[maxn * 2], head[maxn];
edge() { cnp = 1; }
void add(int u, int v)
{
if(u == v) return;
to[cnp] = v, last[cnp] = head[u], head[u] = cnp ++;
to[cnp] = u, last[cnp] = head[v], head[v] = cnp ++;
}
}E1, E2; bool cmp(int x, int y) { return dfn[x] < dfn[y]; } void dfs(int u, int fa)
{
size[u] = 1; dfn[u] = ++ timer;
gra[u][0] = fa; dep[u] = dep[fa] + 1; ST[++ T][0] = dep[u]; pos[u] = T;
for(int i = 1; i < CNST; i ++) gra[u][i] = gra[gra[u][i - 1]][i - 1];
for(int i = E1.head[u]; i; i = E1.last[i])
{
int v = E1.to[i];
if(v == fa) return;
dfs(v, u); size[u] += size[v];
ST[++ T][0] = dep[u];
}
} int LCA(int x, int y)
{
if(dep[x] < dep[y]) swap(x, y);
for(int i = CNST - 1; ~i; i --)
if(dep[gra[x][i]] >= dep[y]) x = gra[x][i];
for(int i = CNST - 1; ~i; i --)
if(gra[x][i] != gra[y][i]) x = gra[x][i], y = gra[y][i];
return (x == y) ? x : gra[x][0];
} int Get_Dis(int x, int y)
{
int ret = dep[x] + dep[y];
if(!x) return dep[y]; else if(!y) return dep[x];
x = pos[x], y = pos[y];
if(x > y) swap(x, y);
int k = Log[y - x + 1];
ret = ret - 2 * min(ST[x][k], ST[y - bit[k] + 1][k]);
return ret;
} int Get_Gra(int x, int y)
{
for(int i = CNST - 1; ~i; i --)
if(dep[gra[x][i]] > dep[y]) x = gra[x][i];
return x;
} int Work(int u, int v, int k)
{
bool flag = 0;
int x = u;
for(int i = CNST - 1; i >= 0; i --)
{
int tem = gra[x][i]; if(dep[tem] < dep[k]) continue;
int d1 = Get_Dis(u, tem), d2 = Get_Dis(v, tem);
if(d1 < d2) x = gra[x][i];
else if(d1 == d2 && u < v) x = gra[x][i];
}
return x;
} void DP(int u, int fa)
{
hson[u] = 0; sum[u] = 0;
for(int i = E2.head[u]; i; i = E2.last[i])
{
int v = E2.to[i]; if(!v) continue;
if(v == fa) continue;
DP(v, u);
int t = hson[v];
int dis = Get_Dis(u, t), dis2 = Get_Dis(u, hson[u]);
if(dis < dis2) hson[u] = t;
else if(dis == dis2 && t < hson[u]) hson[u] = t;
if(!hson[u]) hson[u] = t;
}
if(mark[u]) hson[u] = u;
} void DP2(int u, int fa)
{
int tot = 0;
if(hson[fa])
{
int t = hson[fa];
int dis = Get_Dis(u, t), dis2 = Get_Dis(u, hson[u]);
if(dis < dis2) hson[u] = t;
else if(dis == dis2 && t < hson[u]) hson[u] = t;
if(!hson[u]) hson[u] = t;
}
for(int i = E2.head[u]; i; i = E2.last[i])
{
int v = E2.to[i];
if(v == fa) continue; DP2(v, u);
int y = Get_Gra(v, u);
int x = Work(hson[v], hson[u], u);
int tem = size[x] - size[v];
ans[hson[v]] += tem, ans[hson[u]] += size[y] - size[x];
tot += size[y];
}
ans[hson[u]] += size[u] - tot;
E2.head[u] = 0;
} void Solve()
{
int K = read();
for(int i = 1; i <= K; i ++) b[i] = a[i] = read(), mark[a[i]] = 1;
sort(a + 1, a + 1 + K, cmp); S[top = 1] = 1; E2.cnp = 1; ans[0] = 0;
for(int i = 1; i <= K; i ++)
{
int lca = LCA(S[top], a[i]);
while(2333)
{
if(dep[S[top - 1]] <= dep[lca])
{
E2.add(lca, S[top --]);
if(lca != S[top]) S[++ top] = lca;
break;
}
else E2.add(S[top], S[top - 1]), top --;
}
S[++ top] = a[i];
}
while(top > 1) E2.add(S[top], S[top - 1]), top --;
sum[1] = 0; DP(1, 0); DP2(1, 0);
for(int i = 1; i <= K; i ++) printf("%d ", ans[b[i]]), mark[b[i]] = 0, ans[b[i]] = 0;
puts("");
} void init()
{
bit[0] = 1; for(int i = 1; i <= 20; i ++) bit[i] = bit[i - 1] << 1;
Log[0] = -1; for(int i = 1; i < maxn * 2; i ++) Log[i] = Log[i >> 1] + 1;
} int main()
{
init(); n = read();
for(int i = 1; i < n; i ++)
{
int x = read(), y = read();
E1.add(x, y);
}
dfs(1, 0);
for(int i = 1; i <= Log[T]; i ++)
for(int j = 1; j + bit[i - 1] <= T; j ++)
ST[j][i] = min(ST[j][i - 1], ST[j + bit[i - 1]][i - 1]);
int q = read();
for(int i = 1; i <= q; i ++) Solve();
return 0;
}
【题解】HNOI2014世界树的更多相关文章
- [题解] [HNOI2014] 世界树
题面 [HNOI2014]世界树 题解 从数据范围很容易看出是个虚树DP(可惜看出来了也还是不会做) 虚树大家应该都会, 不会的话自己去搜吧, 我懒得讲了, 我们在这里只需要考虑如何DP即可 首先我们 ...
- [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 ...
随机推荐
- C# 组装XML传给webserver+XML 返回获取多个xml,根据多个XML 返回dataset类型
大致流程介绍: 传值给 webserver+XML ,得到webserver+XML多个返回值,组装成dataset形式返回 首先创建所需要的类型 DataSet ds = new DataSet() ...
- 浅析JVM内存区域及垃圾回收
一.JVM简介 JVM,全称Java Virtual Machine,即Java虚拟机.以Java作为编程语言所编写的应用程序都是运行在JVM上的.JVM是一种用于计算设备的规范,它是一个虚构出来的计 ...
- 原生js常用方法
原生JavaScript设置cookie值 function setCookie(name, value, Hours) { var d = new Date(); var offset = 8; v ...
- leetcode-零钱兑换—int溢出
零钱兑换 给定不同面额的硬币 coins 和一个总金额 amount.编写一个函数来计算可以凑成总金额所需的最少的硬币个数.如果没有任何一种硬币组合能组成总金额,返回 -1. 示例 1: 输入: c ...
- python 文件编译成exe可执行文件。
pyinstaller打包方法: pyinstaller安装参考地址:http://www.pyinstaller.org/ pywin32的下载地址:https://sourceforge.net/ ...
- python leveldb 文档
标签(空格分隔): python leveldb import leveldb db = leveldb.LevelDB('./db') db.Put('hello', 'world') print ...
- 【转载】inotify+rsync实时同步 解决同步慢问题 (转载备记)
原文地址:http://www.ttlsa.com/web/let-infotify-rsync-fast/ 背景 我们公司在用inotify+rsync做实时同步,来解决分布式集群文件一致性的问题. ...
- js如何使浏览器允许脚本异步加载
js如何使浏览器允许脚本异步加载 如果脚本体积很大,下载和执行的时间就会很长,因此造成浏览器堵塞,用户会感觉到浏览器“卡死”了,没有任何响应.这显然是很不好的体验,所以浏览器允许脚本异步加载,下面就是 ...
- 数论的欧拉定理证明 & 欧拉函数公式(转载)
欧拉函数 :欧拉函数是数论中很重要的一个函数,欧拉函数是指:对于一个正整数 n ,小于 n 且和 n 互质的正整数(包括 1)的个数,记作 φ(n) . 完全余数集合:定义小于 n 且和 n 互质的数 ...
- Dom的样式操作和属性操作
如果说web的研究对象是html和css,那么整个dom结构,包含html树和dom树的dom结构才是研究对象,而在整个页面呈现上面,js起到的作用则是异步的用户行为. 按照上面整个思路,获取dom元 ...