hdu 4035 Maze

题目链接：hdu 4035 Maze

　

这个几乎是照着　这篇博客　　写的，写的很好

#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>

using namespace std;
const int maxn=10010;
const double eps=1e-10;
vector<int> g[maxn];
double A[maxn],B[maxn],C[maxn],k[maxn],e[maxn];

int zero(double a)
{
    if(a>-eps&&a<eps) return 1;
    return 0;
}
void dfs(int u,int fa)
{
    double tmp1=0,tmp2=0,tmp3=0;
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        if(v==fa) continue;
        dfs(v,u);
        tmp1+=A[v];
        tmp2+=B[v];
        tmp3+=C[v];
    }
    int m = max((int)g[u].size(),1);
    double dd=(1-k[u]-e[u])/m,dm=(1-dd*tmp2);
    A[u]=(k[u]+dd*tmp1)/dm;
    B[u]=dd/dm;
    C[u]=(dd*tmp3 + (1-k[u]-e[u]))/dm;
}
int main()
{
    int n,u,v,ca=1,T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1;i<n;i++)
        {
            scanf("%d%d",&u,&v);
            g[u].push_back(v);
            g[v].push_back(u);
        }
        for(int i=1;i<=n;i++)
           scanf("%lf%lf",&k[i],&e[i]),k[i]/=100.0,e[i]/=100.0;

        dfs(1,-1);

        if(zero(A[1]-1)) printf("Case %d: impossible\n",ca++);
        else printf("Case %d: %lf\n",ca++,C[1]/(1-A[1]));

        for(int i=1;i<=n;i++) g[i].clear();
    }
    return 0;
}