洛谷 P3233 [HNOI2014]世界树(虚树+dp)
题面
题解
数据范围已经告诉我们是虚树了,考虑如何在虚树上面\(dp\)
以下摘自hzwer博客:
构建虚树以后两遍dp处理出虚树上每个点最近的议事处
然后枚举虚树上每一条边,考虑其对两端点的答案贡献
可以用倍增二分出分界点
如果a,b的分界点为mid,a,b路径上a的第一个儿子为x
则对a的贡献是size[x]-size[mid]
对b的贡献是size[mid]-size[b]
还要算上没被考虑的点
Code
// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define LL long long
#define RG register
using namespace std;
template<class T> inline void read(T &x) {
x = 0; RG char c = getchar(); bool f = 0;
while (c != '-' && (c < '0' || c > '9')) c = getchar(); if (c == '-') c = getchar(), f = 1;
while (c >= '0' && c <= '9') x = x*10+c-48, c = getchar();
x = f ? -x : x;
return ;
}
template<class T> inline void write(T x) {
if (!x) {putchar(48);return ;}
if (x < 0) x = -x, putchar('-');
int len = -1, z[20]; while (x > 0) z[++len] = x%10, x /= 10;
for (RG int i = len; i >= 0; i--) putchar(z[i]+48);return ;
}
int n;
const int N = 300010;
struct node {
int to, next;
}g[N<<1];
int last[N], gl;
inline void add(int x, int y) {
g[++gl] = (node) {y, last[x]};
last[x] = gl;
return ;
}
int dfn[N], cnt, siz[N], dep[N], anc[N][21], rem[N], bel[N];
void init(int u, int fa) {
dfn[u] = ++cnt; siz[u] = 1;
anc[u][0] = fa;
for (int i = 1; i <= 20; i++)
anc[u][i] = anc[anc[u][i-1]][i-1];
for (int i = last[u]; i; i = g[i].next) {
int v = g[i].to;
if (v == fa) continue;
dep[v] = dep[u]+1;
init(v, u);
siz[u] += siz[v];
}
return ;
}
int lca(int x, int y) {
if (dep[x] < dep[y]) swap(x, y);
for (int i = 20; i >= 0; i--)
if (dep[x]-(1<<i) >= dep[y])
x = anc[x][i];
if (x == y) return x;
for (int i = 20; i >= 0; i--)
if (anc[x][i] != anc[y][i])
x = anc[x][i], y = anc[y][i];
return anc[x][0];
}
int dis(int x, int y) {
return dep[x]+dep[y]-2*dep[lca(x, y)];
}
int top, len, m, a[N], b[N], s[N], c[N], f[N];
bool cmp(int a, int b) {
return dfn[a] < dfn[b];
}
inline void insert(int x) {
if (top == 1) {s[++top] = x; return ;}
int o = lca(x, s[top]);
while (top > 1 && dfn[s[top-1]] >= dfn[o]) add(s[top-1], s[top]), top--;
if (o != s[top]) add(o, s[top]), s[top] = o;
s[++top] = x;
return ;
}
void dfs1(int x) {
c[++len] = x; rem[x] = siz[x];
for (int i = last[x]; i; i = g[i].next) {
dfs1(g[i].to);
if (!bel[g[i].to]) continue;
int t1 = dis(bel[g[i].to], x), t2 = dis(bel[x], x);
if ((t1 == t2 && bel[g[i].to] < bel[x]) || t1 < t2 || !bel[x])
bel[x] = bel[g[i].to];
}
return ;
}
void dfs2(int x) {
for (int i = last[x]; i; i = g[i].next) {
int t1 = dis(bel[x], g[i].to), t2 = dis(bel[g[i].to], g[i].to);
if ((t1 == t2 && bel[g[i].to] > bel[x]) || t1 < t2 || !bel[g[i].to])
bel[g[i].to] = bel[x];
dfs2(g[i].to);
}
return ;
}
void solve(int a, int b) {
int x = b, mid = b;
for (int i = 20; i >= 0; i--)
if (dep[anc[x][i]] > dep[a])
x = anc[x][i];
rem[a] -= siz[x];
if (bel[a] == bel[b]) {
f[bel[a]] += siz[x]-siz[b];
return ;
}
for (int i = 20; i >= 0; i--) {
int nxt = anc[mid][i];
if (dep[nxt] <= dep[a]) continue;
int t1 = dis(bel[a], nxt), t2 = dis(bel[b], nxt);
if (t1 > t2 || (t1 == t2 && bel[b] < bel[a])) mid = nxt;
}
f[bel[a]] += siz[x]-siz[mid];
f[bel[b]] += siz[mid]-siz[b];
return ;
}
void query() {
top = len = gl = 0;
read(m);
for (int i = 1; i <= m; i++) read(a[i]), b[i] = a[i];
for (int i = 1; i <= m; i++) bel[a[i]] = a[i];
sort(a+1, a+1+m, cmp);
if (bel[1] != 1) s[++top] = 1;
for (int i = 1; i <= m; i++) insert(a[i]);
for (int i = 1; i < top; i++) add(s[i], s[i+1]);
dfs1(1); dfs2(1);
for (int i = 1; i <= len; i++)
for (int j = last[c[i]]; j; j = g[j].next)
solve(c[i], g[j].to);
for (int i = 1; i <= len; i++) f[bel[c[i]]] += rem[c[i]];
for (int i = 1; i <= m; i++) write(f[b[i]]), putchar(' ');
putchar('\n');
for (int i = 1; i <= len; i++) f[c[i]] = bel[c[i]] = last[c[i]] = 0;
return ;
}
int main() {
read(n);
for (int i = 1; i < n; i++) {
int x, y;
read(x); read(y);
add(x, y); add(y, x);
}
init(1, 0);
memset(last, 0, sizeof(last));
int q; read(q);
while (q--) query();
return 0;
}
洛谷 P3233 [HNOI2014]世界树(虚树+dp)的更多相关文章
- 洛谷3233 HNOI2014(虚树+dp)
膜拜一发\(mts\_246,forever\_shi\) 这两位爷是真的无敌! 首先来看这个题,一看题目的数据范围和"关键点"字眼,我们就能得知这是一道虚树题 那就先一如既往的建 ...
- BZOJ 3572 [HNOI2014]世界树 (虚树+DP)
题面:BZOJ传送门 洛谷传送门 题目大意:略 细节贼多的虚树$DP$ 先考虑只有一次询问的情况 一个节点$x$可能被它子树内的一个到x距离最小的特殊点管辖,还可能被管辖fa[x]的特殊点管辖 跑两次 ...
- ●洛谷P3233 [HNOI2014]世界树
题链: https://www.luogu.org/problemnew/show/P3233题解: 虚树,dp,倍增. 首先对于每个询问,要把虚树建出来,这一步就从略了.这里着重分享一下如何求答案. ...
- bzoj3572[Hnoi2014] 世界树 虚树+dp+倍增
[Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1921 Solved: 1019[Submit][Status][Dis ...
- BZOJ 3572: [Hnoi2014]世界树 [虚树 DP 倍增]
传送门 题意: 一棵树,多次询问,给出$m$个点,求有几个点到给定点最近 写了一晚上... 当然要建虚树了,但是怎么$DP$啊 大爷题解传送门 我们先求出到虚树上某个点最近的关键点 然后枚举所有的边$ ...
- 洛谷P3233 [HNOI2014]世界树
虚树= = #include<cstdio> #include<cstdlib> #include<algorithm> #include<cstring&g ...
- bzoj 3572世界树 虚树+dp
题目大意: 给一棵树,每次给出一些关键点,对于树上每个点,被离它最近的关键点(距离相同被标号最小的)控制 求每个关键点控制多少个点 分析: 虚树+dp dp过程如下: 第一次dp,递归求出每个点子树中 ...
- 【BZOJ】3572: [Hnoi2014]世界树 虚树+倍增
[题意]给定n个点的树,m次询问,每次给定ki个特殊点,一个点会被最近的特殊点控制,询问每个特殊点控制多少点.n,m,Σki<=300000. [算法]虚树+倍增 [题解]★参考:thy_asd ...
- BZOJ 3572: [Hnoi2014]世界树 虚树 树形dp
https://www.lydsy.com/JudgeOnline/problem.php?id=3572 http://hzwer.com/6804.html 写的时候参考了hzwer的代码,不会写 ...
随机推荐
- CSS 内容生成
原文地址:http://www.zhangxinxu.com/wordpress/?p=739 一.哗啦哗啦的简介 zxx://这里“哗啦哗啦”的作用是为了渲染一种氛围.content属性早在CSS2 ...
- spring aop自动代理注解配置之二
<?xml version="1.0" encoding="UTF-8"?> <beans xmlns="http://www.sp ...
- 【转】Provisional headers are shown
在chrome开发者工具的 Network 面板中,某些请求头后面会跟着下面这行文字: Provisional headers are shown 这种请求实际上根本没有产生,对应的请求头当然也不应该 ...
- Makefile 编写规则 - 1
Makefilen内容 1. 显示规则:显示规则说明了,如何生成一个或多个目标.这是由Makefile指出要生成的文件和文件依赖的文件.2. 隐晦规则:基于Makefile的自动推导功能3. 变量的定 ...
- UOJ#22. 【UR #1】外星人
传送门 分析 我们发现一个很神的性质,就是对于一个数如果放在它之前的数小于它那它一定对答案没有贡献 于是我们用dp[i][j]表示从大往小考虑了前i个数,当前答案是j的方案数 我们知道它由两种情况转移 ...
- ServletContext作用功能详解.RP
ServletContext ServletContext,是一个全局的储存信息的空间,服务器开始,其就存在,服务器关闭,其才释放.request,一个用户可有多个:session,一个用户一个:而s ...
- Maven——依赖
一.依赖的范围 依赖的范围有几个可选值,我们用得到的是:compile.test.provided 三个. [1]从项目结构角度理解 compile 和 test 的区别 [2]从开发和运行这两个不同 ...
- Terminologies in MVC: Part 2 (Razor Engine Syntax vs Web Form)
By Abhishek Jaiswal :) on Mar 21, 2015 In this article we learn about Razor Engine Syntax vs Web For ...
- using JSTL
http://docs.oracle.com/javaee/5/tutorial/doc/bnake.html JSTL(JSP Standard Tag Library)
- Java IO输入输出流File 字节流
public static void main(String[] args) throws IOException { //创建源 File f =new File("file/stream ...