【二次扫描与换根】总结

N o . 0 No.0 No.0 前言

前置知识:树形DP。

N o . 1 No.1 No.1 引入&思想

换根 DP 通常与树形 DP 相联系,在进行树形 DP 时,根节点通常是确定的,解决方案大多是从根节点开始递归。但一类题目不会指定根结点,且根结点的变化会对答案产生影响,暴力枚举每一个节点为根节点又会在时间复杂度上出问题,这时可以考虑使用换根 DP 解决。

换根 DP 的思路是使用两次 DFS,第一次确定以某一节点为根节点时的答案,第二次以刚才确定的根节点出发,在每次递归前都自顶向下进行状态转移,用父节点的状态更新子节点的状态,计算出“换根”后的结果。

接下来以例题为背景分析换根 DP 的思路。

N o . 2 No.2 No.2 例题

[POI2008] STA-Station

题意清晰明了,不再重复。

考虑转换根节点会发生什么。

如图,当树的根节点从5变成4后,原以4为根的子树的节点深度全部减一,除开这些节点其他节点深度全部加一。由此推导得出状态转移方程: d p v = d p u − s i v + ( n − s i v ) dp_v = dp_u - si_v + (n - si_v) dpv=dpusiv+(nsiv),其中 u u u v v v 的父亲节点, d p dp dp 代表总深度, s i u si_u siu 代表以 u u u 为根的子树节点个数。

分析第二个 DFS ,也就是换根的实现过程:

void Tr_DP(ll u, ll fa) {
	for (int i = 0; i < Gra[u].size(); i ++){
		ll v = Gra[u][i];
		if (v == fa){
			continue;
		}
		dp[v] = dp[u] - si[v] + (n - si[v]);//根据父节点u求出以v为根时的答案
		Tr_DP(v, u);//已得到以v为根时的答案,根据此继续解决v的子节点
	}
}

[APIO2014] 连珠线

题目大意: 开始只有一个结点。两种操作:可以用一条红线将一个新点连向旧点。或者是选择一条红线,用一个新点和两条蓝线连接原先红线连接的两个点。现在给出最终的连线情况,问蓝线可能的最大长度是多少。

分析:因为始终出现的新节点会只连上一个旧节点,或在两个已经有关系的旧点中间出现,可以分析得出最后的图为一棵树(严谨看待,如果题目没有说明,应该会出现将树认成图的情况)。

所以蓝线的连线情况只可能为两种: f a − u − s o n fa - u - son fauson s o n 1 − u − s o n 2 son_1 - u - son_2 son1uson2

但如图,如果为后者,也就是 e d g e 2 , 3   , e d g e 2 , 6 edge_2,_{3}\ ,edge_2,_{6} edge2,3 ,edge2,6 为蓝边,就可以 3 或 6 号节点为根,将原图转换为前者,所以我们只考虑第一种情况。

考虑设计状态:

d p i , 0 dp_i,_{0} dpi,0 i i i 节点为蓝线中的 u u u 节点时最大值

d p i , 1 dp_i,_{1} dpi,1 i i i 节点只与子节点用红线连接最大值

d p i , 2 dp_i,_{2} dpi,2 i i i 节点为蓝线中的父节点时最大值

不难发现, d p i , 1 , d p i , 2 dp_i,_{1},dp_i,_{2} dpi,1,dpi,2 都与接下来的转移结果没有影响,换句话说, i i i 节点为蓝线中的父节点时可以与多条蓝线连接,所以可以合并为: d p i , 1 dp_i,_{1} dpi,1 i i i 节点不为蓝线中的 u u u 节点时最大值。

对于 d p i , 0 dp_i,_{0} dpi,0 i i i 节点与子节点连线可为红线,也可为蓝线,所以:

d p i , 0 = ∑ v , f a v = u m a x ( d p v , 0 , d p v , 1 + w e i g h t u , v ) dp_i,_{0} = \sum_{v, fa_v = u} max(dp_v,_{0}, dp_v,_{1} + weight_u,_{v}) dpi,0=v,fav=umax(dpv,0,dpv,1+weightu,v)

对于 d p i , 1 dp_i,_{1} dpi,1 i i i 节点与其中一个且仅有一个子节点连线为蓝线肯定选择对答案贡献最大的一个,与其他子节点连线可为红线可为蓝线中的 f a fa fa,所以:

d p i , 1 = d p i , 0 + m a x ( d p v , 0 + w e i g h t u , v − m a x ( d p v , 0 , d p v , 1 + w e i g h t u , v ) ) dp_i,_{1} = dp_i,_{0} + max(dp_v,_{0} + weight_u,_{v} - max(dp_v,_{0}, dp_v,_{1} + weight_u,_{v})) dpi,1=dpi,0+max(dpv,0+weightu,vmax(dpv,0,dpv,1+weightu,v))

记录 m x mx mx 数组为 m a x ( d p v , 0 + w e i g h t u , v − m a x ( d p v , 0 , d p v , 1 + w e i g h t u , v ) ) max(dp_v,_{0} + weight_u,_{v} - max(dp_v,_{0}, dp_v,_{1} + weight_u,_{v})) max(dpv,0+weightu,vmax(dpv,0,dpv,1+weightu,v)) p p p 数组为 v v v 节点的位置, m n mn mn 为次大的 d p v , 0 + w e i g h t u , v − m a x ( d p v , 0 , d p v , 1 + w e i g h t u , v ) dp_v,_{0} + weight_u,_{v} - max(dp_v,_{0}, dp_v,_{1} + weight_u,_{v}) dpv,0+weightu,vmax(dpv,0,dpv,1+weightu,v)

根据状态转移方程可以轻松得出以某一结点为根时的答案,接下来就考虑换根。

考虑两种情况:

情况一:换的子节点不为根节点的 p p p 数组的值,即这个点不是根节点为蓝边中间节点时的儿子,就应该将 d p u , 0 dp_u,_{0} dpu,0 v v v 的贡献减掉,同时,对于新的根节点,同样可以与旧的根节点组成蓝边,需重新计算贡献。

情况二:换的子节点是根节点的 p p p 数组的值,即这个点是根节点为蓝边中间节点时的儿子,那么在减掉此点贡献的同时,需重新计算 d p i , 1 dp_i,_{1} dpi,1,即采用次大值更新。这也是为什么要记录次大值数组 m n mn mn

两次 DFS 代码:

void Tr_dp(ll u,ll fa) {
	mx[u] = mn[u] = -0x3f3f3f3f3f3f;
	for (int i = 0; i < Gra[u].size(); i ++){
		ll v = Gra[u][i].v, w = Gra[u][i].w;
		if (v == fa){
			continue;
		}
		Tr_dp(v, u);
		ll k = max(dp[v][0], dp[v][1] + w);
		dp[u][0] += k;
		if (dp[v][0] + w - k > mx[u]){
			mn[u] = mx[u];
			mx[u] = dp[v][0] + w - k;
			p[u] = v;
		}
		else if (dp[v][0] + w - k > mn[u]){
			mn[u] = dp[v][0] + w - k;
		}
	}
	dp[u][1] = dp[u][0] + mx[u];
}

void Tr_DP(int u, int fa, int val){
	int k = max(dp2[fa][0], dp2[fa][1] + val);
	dp[u][0] += k;
	if (fa && dp2[fa][0] + val - k > mx[u]){
		mn[u] = mx[u];
		mx[u] = dp2[fa][0] + val - k;
		p[u] = fa;
	}
	else if (fa && dp2[fa][0] + val - k > mn[u]){
		mn[u] = dp2[fa][0] + val - k;
	}
	dp[u][1] = dp[u][0] + mx[u];
	ans = max(ans, dp[u][0]);
	for (int i = 0; i < Gra[u].size(); i ++){
	 	ll v = Gra[u][i].v, w = Gra[u][i].w;
		if (v == fa){
			continue;
		}
		dp2[u][0] = dp[u][0] - max(dp[v][0], dp[v][1] + w);
		if (v == p[u]){
			dp2[u][1] = dp2[u][0] + mn[u];
		}
		else {
			dp2[u][1] = dp2[u][0] + mx[u];
		}
		Tr_DP(v, u, w);
	}
}

[SHOI2014] 概率充电器

前置知识:概率

大概题意:给你一棵树,树上每一点有直接通电的概率,每条边有导电的概率。求出期望有电节点数。

分析:

不难发现能使一个节点 u u u 导电有三种情况:

节点 u u u 自己生电。

儿子节点给节点 u u u 导电。

父亲节点给节点 u u u 导电。

难点来了:如何状态转移使 DP 没有后效性?

因为我们发现根节点不受父亲的影响,考虑先忽略父亲节点给节点 u u u 导电的概率,求出另两种情况使节点 u u u 导电的概率,也就是从叶结点出发到根节点进行递归。然后再从根节点出发考虑父节点的贡献。

d p u : dp_u: dpu: u u u 个节点在考虑自身生电的概率其子树导电的概率以后,没电的概率。

u u u 节点的一个子节点 v v v 单独拿出来考虑,发现 v v v 节点不能给 u u u 节点导电有两种情况:自身没电以及路上断电,相加即为节点 u u u 不会受到 v v v 的电流影响的概率。

d p u = ( 1 − c o u ) ∏ v , f a v = u ( d p v + ( 1 − d p v ) ∗ ( 1 − w e i g h t u , v ) ) dp_u = (1 - co_u) \prod_{v, fa_v = u} (dp_v + (1 - dp_v) * (1 - weight_u,_v)) dpu=(1cou)v,fav=u(dpv+(1dpv)(1weightu,v))

其中 c o u co_u cou 即为 u u u 能单独生电的概率。

所以我们能够发现除开节点 v v v ,节点 u u u 不能生电的概率为 d p u d p v + ( 1 − d p v ) ( 1 − w e i g h t u , v ) \frac{dp_u}{dp_v + (1 - dp_v)(1 -weight_u,_v)} dpv+(1dpv)(1weightu,v)dpu

然后考虑 u u u 节点父节点对 u u u 的影响。

d p 2 u : dp2_u: dp2u: u u u 个节点在考虑其父节点导电的概率以后,没电的概率。

所以父节点除开 u u u 节点的影响后没电的概率为 d p 2 u ∗ d p u d p u + ( 1 − d p u ) ( 1 − w e i g h t f a u , u ) dp2_u * \frac{dp_u}{dp_u + (1 - dp_u)(1 -weight_{fa_u},_u)} dp2udpu+(1dpu)(1weightfau,u)dpu, 设为 p p p

所以 d p 2 u dp2_u dp2u,即第 u u u 个节点在考虑其父节点导电的概率以后,没电的概率就为父节点本身没电的概率乘上路上断电的概率。

所以 d p 2 u = p + ( 1 − p ) ( 1 − w e i g h t f a u , u ) dp2_u = p + (1 - p)(1 - weight_{fa_u},_u) dp2u=p+(1p)(1weightfau,u)

所以 u u u 节点没电的概率就为 d p u ∗ d p 2 u dp_u * dp2_u dpudp2u,有电的概率就为 1 − d p u ∗ d p 2 u 1 - dp_u * dp2_u 1dpudp2u

最终的答案即使所有节点有电的概率和。

P S : PS: PS: 注意运算时分母不能为 0。

两次 DFS 代码:

void Tr_dp(ll u, ll fa) {
	dp[u] = 1.0 - co[u] * 1.0 / 100;
	for (int i = 0; i < Gra[u].size(); i ++){
		ll v = Gra[u][i].v;
		db w = Gra[u][i].w * 1.0 / 100;
		if (v == fa){
			continue;
		}
		Tr_dp(v, u);
		dp[u] *= (dp[v] + (1 - dp[v]) * (1 - w));
	}
}

void Tr_DP(ll u, ll fa) {
	db te;
	for (int i = 0; i < Gra[u].size(); i ++){
		ll v = Gra[u][i].v;
		db w = Gra[u][i].w * 1.0 / 100;
		if (v == fa){
			db t;
			if (dp[u] + (1 - dp[u]) * (1 - w) == 0){
				t = dp2[fa] * dp[fa];
			}
			else {
				t = dp2[fa] * dp[fa] / (dp[u] + (1 - dp[u]) * (1 - w));
			}
			dp2[u] = (t + (1 - t) * (1 - w));
			break;
		}
	}
	if (fa == 0){
		dp2[u] = 1;
	}
	for (int i = 0; i < Gra[u].size(); i ++){
		ll v = Gra[u][i].v;
		if (v == fa){
			continue;
		}
		Tr_DP(v, u);
	}
}

N o . 3 No.3 No.3 总结

二次扫描与换根的题目算 DP,需要大量的思维难度,也需要笔和草稿纸的配合。

努力ing

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值