SHOI2014充電確率（DP木）（DPの確率）（および除外）

ポータル

良い質問ああ！それは彼の父から再びカウント、その後、再びサブツリーの数から明らかです

メイク $F（U）$確率のサブツリーからの電力で電気または電荷に直接電荷を表す
合併の確率はどの包含排除を必要と $P（A | B）= P（A）+ P（B）-P（A）* P（B）$

それから $F（U）$拡張子はに父から渡される
、それがダウンして行くために自分自身に着手していない彼の父親から来た場合には、明らかに

その確率は父と父の他の息子を考慮した確率であります
$P（A）= \ FRAC {P（A | B）-P（B）} {1-P（B）}$ それは息子の父親との確率で他の電荷から父であります

そして、父合併静止および除外の合併から下ります

木の複合効果 $DP$ルーチンと確率 $DP$コンテンツ転送撥

#include<bits/stdc++.h>
#define cs const
using namespace std;
cs int N = 5e5 + 5;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt * f;
}
int first[N], nxt[N << 1], to[N << 1], tot; double w[N << 1];
int n; double a[N], f[N];
void add(int x, int y, double z){nxt[++tot] = first[x], first[x] = tot, to[tot] = y, w[tot] = z; }
void dfs1(int u, int fa){
	f[u] = a[u];
	for(int i = first[u]; i; i = nxt[i]){
		int t = to[i]; if(t == fa) continue;
		dfs1(t, u); double p = w[i] * f[t];
		f[u] = f[u] + p - f[u] * p;
	}
}
void dfs2(int u, int fa){
	for(int i = first[u]; i; i = nxt[i]){
		int t = to[i]; if(t == fa) continue;
		double pa = f[t] * w[i];
		if(pa < 1){
			double pc = (f[u] - pa) / (1 - pa);
			double p = pc * w[i];
			f[t] = f[t] + p - p * f[t];
		} dfs2(t, u);
	}
}
int main(){
	n = read();
	for(int i = 1; i < n; i++){
		int x = read(), y = read(); double z = read(); z = z/100.0;
		add(x, y, z); add(y, x, z);
	} 
	for(int i = 1; i <= n; i++) a[i] = 1.0 * read() / 100.0;
	dfs1(1, 0); dfs2(1, 0); double ans = 0;
	for(int i = 1; i <= n; i++) ans += f[i];
	printf("%.6lf", ans); return 0;
}