题解 P3233 [HNOI2014]世界树
解题思路
正解当然是虚树了。
首先对于原树以及虚树各开一个结构体存边,这个不用多说。
然后我们先 DFS 一遍,求出各个节点的时间戳,子树大小,深度以及父亲节点,并初始化倍增 LCA 。
对于每一次的操作,我们都建一棵虚树(注意数组的清空),为了方便,我们此后操作的 DFS 都从 1 节点开始,如果 1 节点不是临时议事处,我们也把它给加入到虚树中并做一下标记。
显然,在建虚树之前,要先把所有的临时议事处按照时间戳排序,然后建虚树。
此后没,我们进行两遍 DFS :
第一遍:用于更新虚树中所有节点的 dp 以及其所属的临时议事处。
因为虚树的深度较小的点不一定是临时议事处,而所有的叶子节点一定是临时议事处。
所以,先遍历到叶子节点,再用叶子节点来更新父亲节点。
第二遍:利用所有的虚树中的值进而求出原树中的值。
如果子树中没有临时议事处,那么这个节点管辖的点数就要加上子树的siz。
在原树上是一条链,并且两端都是临时议事处,那么直接一半一半分就好了。
- 注意:此过程也需要更新虚树上点的 dp 值以及所属临时议事处,但是不同与第一次 DFS 的是由父亲节点更新叶子节点。
最后输出时注意判断 1 节点是否是临时议事处就好了。
code
#include<bits/stdc++.h>
using namespace std;
const int N=3e5+10,M=N<<1,INF=1e9;
int n,m,Q,tim,s[N],dp[N],g[N],dfn[N],q[N],ans[N],f[N][25],siz[N],dep[N];
int top,sta[N];
bool flag,vis[N];
struct Edge
{
int tot,head[N],nxt[M],ver[M];
void add(int x,int y)
{
ver[++tot]=y;
nxt[tot]=head[x];
head[x]=tot;
}
}e1,e2;
void dfs(int x,int fat)
{
f[x][0]=fat;
siz[x]=1;
dfn[x]=++tim;
for(int i=e1.head[x];i;i=e1.nxt[i])
{
int to=e1.ver[i];
if(to==fat)
continue;
dep[to]=dep[x]+1;
dfs(to,x);
siz[x]+=siz[to];
}
}
void LCA_init()
{
dep[1]=1;
dfs(1,0);
for(int i=1;i<=20;i++)
for(int j=1;j<=n;j++)
f[j][i]=f[f[j][i-1]][i-1];
}
int LCA_ask(int x,int y)
{
if(x==y)
return x;
if(dep[x]>dep[y])
swap(x,y);
for(int i=20;i>=0;i--)
if(dep[f[y][i]]>=dep[x])
y=f[y][i];
if(x==y)
return x;
for(int i=20;i>=0;i--)
if(f[x][i]!=f[y][i])
{
x=f[x][i];
y=f[y][i];
}
return f[x][0];
}
bool comp(int x,int y)
{
return dfn[x]<dfn[y];
}
void build(int x)
{
if(!top)
{
sta[++top]=x;
return ;
}
int lca=LCA_ask(x,sta[top]);
while(top>1&&dep[lca]<dep[sta[top-1]])
{
e2.add(sta[top-1],sta[top]);
e2.add(sta[top],sta[top-1]);
top--;
}
if(dep[lca]<dep[sta[top]])
{
e2.add(lca,sta[top]);
e2.add(sta[top],lca);
top--;
}
if(!top||sta[top]!=lca)
sta[++top]=lca;
sta[++top]=x;
}
void dfs1(int x,int fa)
{
dp[x]=INF;
for(int i=e2.head[x];i;i=e2.nxt[i])
{
int to=e2.ver[i];
if(to==fa)
continue;
dfs1(to,x);
int dis=dep[to]-dep[x];
if(dp[x]>dp[to]+dis)
{
g[x]=g[to];
dp[x]=dp[to]+dis;
}
else if(dp[x]==dp[to]+dis)
g[x]=min(g[x],g[to]);
}
if(vis[x])
{
dp[x]=0;
g[x]=x;
}
}
void work(int x,int y)
{
int b=y;
for(int i=20;i>=0;i--)
{
int l=dep[y]-dep[f[b][i]]+dp[y],r=dep[f[b][i]]-dep[x]+dp[x];
if(dep[f[b][i]]>dep[x]&&(l<r||(l==r&&g[y]<g[x])))
b=f[b][i];
}
ans[g[y]]+=siz[b]-siz[y];
ans[g[x]]-=siz[b];
}
void dfs2(int x,int fa)
{
for(int i=e2.head[x];i;i=e2.nxt[i])
{
int to=e2.ver[i];
if(to==fa)
continue;
int dis=dep[to]-dep[x];
if(dp[to]>dp[x]+dis)
{
g[to]=g[x];
dp[to]=dp[x]+dis;
}
else if(dp[to]==dp[x]+dis)
g[to]=min(g[to],g[x]);
work(x,to);
dfs2(to,x);
}
ans[g[x]]+=siz[x];
}
int main()
{
scanf("%d",&n);
for(int i=1,x,y;i<n;i++)
{
scanf("%d%d",&x,&y);
e1.add(x,y);
e1.add(y,x);
}
LCA_init();
scanf("%d",&Q);
while(Q--)
{
flag=top=e2.tot=0;
memset(vis,false,sizeof(vis));
memset(e2.head,0,sizeof(e2.head));
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
scanf("%d",&q[i]);
vis[q[i]]=true;
ans[q[i]]=0;
}
if(!vis[1])
{
flag=true;
q[++m]=1;
}
for(int i=1;i<=m;i++)
s[i]=q[i];
sort(s+1,s+m+1,comp);
for(int i=1;i<=m;i++)
build(s[i]);
if(top)
while(--top)
{
e2.add(sta[top],sta[top+1]);
e2.add(sta[top+1],sta[top]);
}
dfs1(1,0);
dfs2(1,0);
for(int i=1;i<=m;i++)
if(q[i]!=1||!flag)
printf("%d ",ans[q[i]]);
printf("\n");
}
return 0;
}
题解 P3233 [HNOI2014]世界树的更多相关文章
- ●洛谷P3233 [HNOI2014]世界树
题链: https://www.luogu.org/problemnew/show/P3233题解: 虚树,dp,倍增. 首先对于每个询问,要把虚树建出来,这一步就从略了.这里着重分享一下如何求答案. ...
- 洛谷 P3233 [HNOI2014]世界树(虚树+dp)
题面 luogu 题解 数据范围已经告诉我们是虚树了,考虑如何在虚树上面\(dp\) 以下摘自hzwer博客: 构建虚树以后两遍dp处理出虚树上每个点最近的议事处 然后枚举虚树上每一条边,考虑其对两端 ...
- 洛谷P3233 [HNOI2014]世界树
虚树= = #include<cstdio> #include<cstdlib> #include<algorithm> #include<cstring&g ...
- luogu P3233 [HNOI2014]世界树
传送门 我是什么时候写的这题的qwq 首先,发现关键点的总数被限制了,很自然想到虚树,并且,对于一个关键点,他管理的点显然是一个联通块 然后把虚树先建出来,然后两次dfs,第一次是向祖先更新离每个点最 ...
- 【题解】HNOI2014世界树
脑子不清醒的时候千万别写题.写题写不下去了千万别死扛,重构才是你唯一的出路QAQ 昨天很想快点写道题,思路没有很清晰的时候就写了,结果……今天一怒之下决定重整思路重构代码,其实不过是半个小时的事情…… ...
- P3233 [HNOI2014]世界树
传送门 看到指定的总节点数小于等于 300000 就知道要搞虚树了 考虑如何在虚树确定每个议事处控制的节点数量 可以两遍dfs 第一遍求儿子对父亲的影响,第二遍求父亲对儿子影响 注意搜索顺序,这样就可 ...
- [题解] [HNOI2014] 世界树
题面 [HNOI2014]世界树 题解 从数据范围很容易看出是个虚树DP(可惜看出来了也还是不会做) 虚树大家应该都会, 不会的话自己去搜吧, 我懒得讲了, 我们在这里只需要考虑如何DP即可 首先我们 ...
- [BZOJ3572][Hnoi2014]世界树
[BZOJ3572][Hnoi2014]世界树 试题描述 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森,在他们的信条 ...
- 【BZOJ3572】[Hnoi2014]世界树 虚树
[BZOJ3572][Hnoi2014]世界树 Description 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森 ...
随机推荐
- Kafka源码分析系列-目录(收藏不迷路)
持续更新中,敬请关注! 目录 <Kafka源码分析>系列文章计划按"数据传递"的顺序写作,即:先分析生产者,其次分析Server端的数据处理,然后分析消费者,最后再补充 ...
- webpack 快速入门 系列 —— 实战一
实战一 准备本篇的环境 虽然可以仅展示核心代码,但笔者认为在一个完整的环境中边看边做,举一反三,效果更佳. 这里的环境其实就是初步认识 webpack一文完整的示例,包含 webpack.devSer ...
- [论文阅读笔记] Fast Network Embedding Enhancement via High Order Proximity Approximati
[论文阅读笔记] Fast Network Embedding Enhancement via High Order Proximity Approximation 本文结构 解决问题 主要贡献 主要 ...
- 转: inline关键字使用
1.inline用在函数声明时,还是函数定义时?还是两边都加? 首先,内联函数声明和定义最好在同一个文件中,其它的情况没有实用上的意义. 只要在同一个文件中,声明和定义至少其一加"inlin ...
- [Python] 条件 & 循环
条件语句 不加 () 结尾加 : elif else 和 if 成对使用 省略判断条件 String:空字符串为False,其余为True int:0为False,其余为True Bool:True为 ...
- UltraISO制作启动盘安装CentOS7
UltraISO制作启动盘安装CentOS7 发表于 2020-03-10 | 分类于 DevOps | 没有评论 简单几个步骤即可完成启动盘的制作,非常便利 准备工具 准备8G优盘(启动盘制 ...
- sed -i '14s/yes/no/' tftp
修改tftp 内容 # cd /etc/xinetd.d/[root@localhost xinetd.d]# cp tftp tftp.bak[root@localhost xinetd.d]# c ...
- 010.Ansible_palybook 循环语句
Ansible循环语句 1 简介 我们在编写playbook的时候,不可避免的要执行一些重复性操作,比如指安装软件包,批量创建用户,操作某个目录下的所有文件等.正如我们所说,ansible一门简单的自 ...
- 关于Ajax 的 cache 属性 (Day_34)
最近做项目,在某些页面显示,ajax刷新总是拿不到新内容,时常需要清除缓存,才能到达想要的效果. 经过再次查看文档,最后加了一行属性:cache:false 即可解决问题 我们先看下文档的说明: 可以 ...
- protege 构建本体
这里我们使用的是Protégé-OWL规范. 推理机后的内容主要是实操内容,根据推理机来对protege本体模型的一个操作过程,以加深本体模型的一个规范认识. 一.什么是本体(Ontologie) 本 ...