HDU4035(概率dp)

解题思路：一开始做的时候，推出来的概率递推表达式，

 是这样的   E[i] = ki*E[1] + ei*0 + (1-ki-ei)/m*( E[father[i]]+1 + ∑( E[child[i]]+1 ) );
         = ki*E[1] + (1-ki-ei)/m*E[father[i]] + (1-ki-ei)/m*∑(E[child[i]]) + (1-ki-ei);

因为子节点和父节点是相互影响的，所以想当然就是用高斯消元法来做。但复杂度是O(n^3)，会超时。后来看了一下题解：

https://www.cnblogs.com/kuangbin/archive/2012/10/03/2711108.html

发现竟然可以这么巧妙，递推关系式只与子节点dp以及父节点的系数有关，而和父节点dp无关！！！具体推导过程详见上面那篇博客。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath> 
#include<iostream>
#define N 10009
using namespace std;

const double eps=1e-10;
vector<int> s[N];
double e[N],k[N],a[N],b[N],c[N];
bool dfs(int u,int pre)
{
	bool fg=true;
	int m=s[u].size();
    a[u]=k[u];
    b[u]=(1-k[u]-e[u])/m;
    c[u]=1-k[u]-e[u];
    double tp=0;
	for(int i=0;i<m;i++)
	{
		int v=s[u][i];
		if(v==pre) continue;
		fg=false;
		if(!dfs(v,u)) return false;
		a[u]+=(1-k[u]-e[u])/m*a[v];
		c[u]+=(1-k[u]-e[u])/m*c[v];
		tp+=(1-k[u]-e[u])/m*b[v];
	}
	if(fabs(tp-1)<eps) return false;
	a[u]/=(1-tp);
	b[u]/=(1-tp);
	c[u]/=(1-tp);
	return true;
}
int main()
{
	//freopen("t.txt","r",stdin);
	int T,p1,p2,n;
	scanf("%d",&T);
	for(int ca=1;ca<=T;ca++)
	{
		printf("Case %d: ",ca);
		scanf("%d",&n);
		for(int i=1;i<=n;i++) s[i].clear();
		for(int i=0;i<n-1;i++)
		{
			scanf("%d%d",&p1,&p2);
			s[p1].push_back(p2);
			s[p2].push_back(p1);
		}
		for(int i=1;i<=n;i++)
		{
			scanf("%lf%lf",&k[i],&e[i]);
			k[i]/=100;
			e[i]/=100;
		}
		if(dfs(1,-1)&&fabs(a[1]-1)>eps)
		{
			//cout<<a[1]<<" "<<c[1]<<" ";
			printf("%.6lf\n",c[1]/(1-a[1]));
		}
		else
		printf("impossible\n");
	}
	return 0;
}