【题解】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 ...
随机推荐
- netty之粘包分包的处理
1.netty在进行字节数组传输的时候,会出现粘包和分包的情况.当个数据还好,如果数据量很大.并且不间断的发送给服务器,这个时候就会出现粘包和分包的情况. 2.简单来说:channelBuffer在接 ...
- Unity Container中的几种注册方式与示例
1.实例注册 最简单的注册方式就是实例注册,Unity 容器负责维护对一个类型的单例引用,比如: 有如下的实际类型: namespace ConsoleSample { public class Sa ...
- 05-JVM对象探秘
一.对象的内存布局 以Hotspot虚拟机为例,对象在内存中的结构可以分为三部分:对象头(header).实例数据(instance data).对齐填充(padding). 1.1. ...
- MySQL☞数值处理函数
1.round():四舍五入函数 round(数值,参数):如果参数的值为正数,表示保留几位小数,如果参数的值为0,则只保留正数部分们如果参数的值为负数,表示对小数点前第几位进行四舍五入. Eg:(1 ...
- 通过 Python_Faker 生成测试数据
通过 Python_Faker 生成测试数据 一.介绍 在软件需求.开发.测试过程中,有时候需要使用一些测试数据,针对这种情况,我们一般要么使用已有的系统数据,你不可能通过手工来生成(最傻的方法)可能 ...
- Anyproxy抓包工具
1.安装Nodejs和AnyProxy以及安装模拟器(移动端抓包工具) 第一步:nodejs下载地址:http://nodejs.cn/download/ 下载Windows版本,直接运行安装即可, ...
- lintcode735. Replace With Greatest From Right
Given an array of integers, replace every element with the next greatest element (greatest element o ...
- Spring Boot - Filter实现简单的Http Basic认证
Copy自http://blog.csdn.net/sun_t89/article/details/51916834 @SpringBootApplicationpublic class Spring ...
- 如何让thinkpad X1C 用U盘 安装上专业版win10
1 BIOS内置了文件 会导致win10 iso默认装家庭版 2 给iso 的resouse 目录中增加文件ei.cfg 3 内容如下 [EditionID]Professional[Channel] ...
- POJ 3525/UVA 1396 Most Distant Point from the Sea(二分+半平面交)
Description The main land of Japan called Honshu is an island surrounded by the sea. In such an isla ...