HDU4035:Maze

题目大意

传送门

思路

我们先考虑最朴素无华的dp，我们假设$dp[i]$为从i点走出迷宫的期望步数，那么，我们可以费劲九牛二虎之力 轻松地列出dp式，如果i是叶子节点，那么

$dp[i]=k_i\times dp[1]+e_i\times 0+(1-e_i-k_i)\times (dp[fa]+1)$$\Rightarrow dp[i]=k_i\times dp[1]+(1-e_i-k_i)\times (dp[fa]+1)$,

这句话什么意思呢？就是说它有$k_i$的概率被杀死回到1点，它也有$e_i$的可能走出迷宫，但是它是从i点走出去，所以它的贡献就是$e_i\times 0$不就是没有么？还在这里废话然后这个点还有$(1-e_i-k_i)$的概率继续走下去，贡献是$(1-e_i-k_i)\times (dp[fa]+1)$,所以这就是叶子结点的情况。对于不是叶子节点的情况也差不多，不过某些部分需要除以它的度数(它走到每个相邻的点的概率都是一样的),所以大概就长成这个亚子:

$dp[i]=k_i\times dp[1]+\frac{(1-k_i-e_i)\times (dp[fa]+1)}{deg[i]}+\frac{{\sum }_{son}dp[son]+1}{deg[i]}$

$\Rightarrow$$dp[i]=k_i\times dp[1]+\frac{(1-e_i-k_i)\times (dp[fa]+1+{\sum }_{son}(dp[son]+1))}{deg[i]}$,

因为跟上面差不多，所以我就不详说了，但是我们会发现一个BUG，这个东西TMD 一直在来回地推，完全就是有递无归，我们难道要用高斯消元了么？不！这个TMD可能么？拜托，$10000^3$,我们考虑一下树形动规，我们在后序遍历中，我们的$dp[son]$是一个已知量，而$dp[fa]$和$dp[1]$是未知量,所以，我们就不做了我们考虑把$dp[i]$表示为另一种形式，就长成

$dp[i]=A[i]\times dp[fa]+B[i]\times dp[1]+C[i]$,
其中$A[i],B[i],C[i]$都是我们自己要算出来的系数项。我们将

$dp[son]=A[son]\times dp[i]+B[son]\times dp[1]+C[son]$带入$dp[i]$的表达式中，我们设$P[i]=\frac{1-e_i-k_i}{deg[i]}$
因为

$dp[i]=k_i\times dp[1]+P[i]\times (dp[fa]+1)+P[i]\times {\sum }_{son}(dp[son]+1)$,

也就是说

$dp[i]=k_i\times dp[1]+P[i]\times (dp[fa]+1)+P[i]\times {\sum }_{son}(A[son]\times dp[i]+B[son]\times dp[1]+C[son])$
在我们合并同类项之后就长成这个亚子:

$(1-{\sum }_{son}A[son]\times P[i])\times dp[i]=dp[fa]\times P[i]+({\sum }_{son}B[son]\times P[i]+k_i)\times dp[1]+{\sum }_{son}(C[son]+1)\times P[i]+P[i]$
我们看着这个东西感觉非常的不爽，那我们就设$T[i]=1-{\sum }_{son}(A[son]\times P[i])$,那么根据之前$dp[i]=A[i]\times dp[fa]+B[i]\times dp[1]+C[i]$,所以我们两边同时除以$T[i]$,就得出:

$A[i]=\frac{P[i]}{T[i]},B[i]=\frac{{\sum }_{son}(B[son]\times P[i]+K[i])}{T[i]},C[i]=\frac{{\sum }_{son}(C[son]\times P[i]+deg[i]\times P[i])}{T[i]}$
所以，我们就可以跑一边树形动规就可以求出每个点的$A[i],B[i],C[i]$,我们发现点1它没有爸爸(它的妈妈是被强奸了么？)，那我们就可以直接算出来，点1的dp值就是$\frac{C[1]}{1.0-B[1]}$,好，上面说了一大啪啦，但是，我们怎么判断有无解呢？实际上，我们只有除以0的时候才无解，所以在除以的时候判断一下就好了

代码理论AC

#include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;

#define Int register int
#define MAXN 10005

int n;
int deg[MAXN];
double E[MAXN],K[MAXN],P[MAXN];
double A[MAXN],B[MAXN],C[MAXN];

vector <int> G[MAXN];

bool pd;

double Abs (double x){return x > 0 ? x : -x;}

void dfs (int u,int fa)
{
	double tmp = 0;
	double Au = P[u],Bu = K[u],Cu = 1 - K[u] - E[u];
	for (Int i = 0;i < G[u].size();++ i)
	{
		int v = G[u][i];
		if (v == fa) continue;
		dfs (v,u);
		if (pd) return ;
		Bu += B[v] * P[u];
		Cu += C[v] * P[u];
		tmp += A[v] * P[u];
	} 
	tmp = 1.0 - tmp;
	if (Abs (tmp) < 1e-10) 
	{
		pd = 1;
		return ;
	}
	A[u] = Au / tmp,B[u] = Bu / tmp,C[u] = Cu / tmp;
}

void read (int &x)
{
	x = 0;char c = getchar();int f = 1;
	while (c < '0' || c > '9'){if (c == '-') f = -f;c = getchar();}
	while (c >= '0' && c <= '9'){x = (x << 3) + (x << 1) + c - '0';c = getchar();}
	x *= f;return ;
}

void write (int x)
{
	if (x < 0){x = -x;putchar ('-');}
	if (x > 9) write (x / 10);
	putchar (x % 10 + '0');
}

signed main()
{
	int times,tot = 0;
	read (times);
	while (times --)
	{	
		pd = 0;
		memset (G,0,sizeof (G));
		memset (E,0,sizeof (E));
		memset (K,0,sizeof (K));
		memset (P,0,sizeof (P));
		memset (A,0,sizeof (A));
		memset (B,0,sizeof (B));
		memset (C,0,sizeof (C));
		memset (deg,0,sizeof (deg));
		read (n);
		for (Int i = 1;i < n;++ i)
		{
			int u,v;
			read (u),read (v);
			G[u].push_back(v);
			G[v].push_back(u);
			deg[u] ++,deg[v] ++;
		}
		for (Int i = 1;i <= n;++ i)	 
		{
			scanf ("%lf%lf",&K[i],&E[i]);
			K[i] /= 100,E[i] /= 100;
		}
		for (Int i = 1;i <= n;++ i)
			P[i] = (1 - E[i] - K[i]) * 1.0 / deg[i];
		dfs (1,0);
		printf ("Case %d: ",++ tot);
		if (pd || Abs (1.0 - B[1]) < 1e-10) puts ("impossible");
		else printf ("%.6lf\n",C[1] / (1.0 - B[1]));
	}
}
/*
3
1 2
1 3
0 0
100 0
0 100
*/
//建议零的值设小一点，误差就小一点