虚树详解

虚树是一种优化树形DP的方法,通过构建包含询问点和不超过询问点数减一个LCA的树来减少点的数量。算法主要应用于处理组询问,如每次询问w个数时,总询问次数为300000。建树过程按DFN排序后依次插入,保证深度关系并避免无用节点。时间复杂度为O(∑i=1qwi×logn)。清除图可以使用DFS清空头指针或迭代实现。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

虚树大意:通过建一棵只包含询问点和不超过询问点数减一个lca的树来减少点的个数,降低时间复杂度

这东西不难,常用于辅助树形dp,难点在于dp…
使用该算法的标志为不确定组的询问次数与给出的询问点数和 例如:
每次询问w个数, ∑ i = 1 q w i = 300000 \sum_{i=1}^qw_i=300000\quad i=1qwi=300000


算法:

构建过程十分简单,放心食用
举个例子:P2495 [SDOI2011]消耗战
就是找一些边割掉使得根与所有选定点不连通,使割掉的边权和最小
首先给出一棵树:(原谅我丑陋的图)

在这里插入图片描述
红色的点是我们要询问的,明显我们不可能直接在原树上搞O(n)dp,我要是全部点都选不就gg了。。
所以我们考虑省掉一些点来优化总点数。

很容易观察到,3和7号节点并没有对答案造成贡献,如果只是找最小边权只需要用一个前缀记录就可以,那么如何优化呢?

明显我们需要保证虚树仍保持节点深度关系的同时保证不加入无用点。对选定点按照dfn排序再依次插入,用栈记录当前的一条链,那么可见我们新加入一个点x会有以下两种情况:

  1. 栈顶是x的祖先,直接加进去
  2. 栈顶不是x的祖先,求栈顶和x的lca,不停出栈,并连一条stk[top-1]->stk[top]的边(构建当前链),直到找到栈内最后一个点的dfn>=lca的dfn,即栈顶上一个点的dfn<lca的dfn。如果栈顶的点不是lca,则连一条lca到栈顶的边,再把栈顶改成lca,加入x,保持深度递增。

其实我们的虚树的深度就是当前根到x路径上的选定点个数,加入lca是因为要使得同一深度的点保持这种关系。举个栗子,x和y的深度都是4,你连x->y或者y->x都会改变原树形态,此时就要lca->x,lca->y。
我们又按照了dfn排序,此时选定点有序,只需要处理相邻两点之间的位置问题即可。

注意一下,第一个点永远置为1,方便搜索;最后将栈内的全部出栈加边。

可以通过上面的图理解:

首先点排序后长这样:4 10 11 9 5

  1. stk:1 (初始)
  2. stk:1 4 此时加入了4,栈顶1是4的祖先
  3. stk:1 4 10 同上
  4. stk:1 4 7 11 此时加入11,11和栈顶10的祖ca是7,而4的深度小于七,那么把栈顶改为7,加入11
  5. stk:1 4 9 此时加入9,9和栈顶11的lca是4,将11和7出栈,此时栈顶就是lca,无需改动,加入9
  6. stk:1 5 此时加入5,5和栈顶9的lca是1,将4和9出栈,此时栈顶就是lca,无需改动,加入5
  7. stk:1 此时将栈内元素出完

现在看看建出的树:
1 无
2 无
3 无
4 7->10
5 7->11 4->7
6 4->9 1->4
7 1->5

画出来长这样 :

在这里插入图片描述
少了一堆点,效率up

时间复杂度分析:

刚才讲了只需要处理相邻两点的lca,即加入不超过w-1(w是询问点数)个lca。
那么每次询问只需要加入不超过2w-1个点。

总时间复杂度为 O ( ∑ i = 1 q w i × l o g n ) O(\sum_{i=1}^qw_i\times logn) O(i=1qwi×logn),于是就可以快乐地解决问题了。


技巧:

建完虚树清除图是个问题,毕竟也不可能memset,瞬间上天。。。
我们可以dfs一边,每到一个点就清除其head,就不需要耗费过多时间。

还有一种方法,是类似迭代实现,就不需要建边。我们在建树过程中其实就是在遍历虚树,那么我们只要把加边操作替换为更新父节点信息,出栈时更新ans或者直接输出一号节点信息(看题目)即可。


代码:

我是用树剖求lca的。

#include<bits/stdc++.h>
using namespace std;
int dfn[250010],dep[250010],siz[250010],n,q;
int top[250010],fa[250010],stk[250010],son[250010],tot,tp;
int e[500010],w[500010],nxt[500010],head[250010],cnt;
int a[250010];
void add(int x,int y,int z){
   
	cnt++;
	e[cnt]=y;
	w[cnt]=z;
	nxt[cnt]=head[x];
	head[x]=cnt;
}</
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值