HDU4035 Maze 【树形DP】【期望DP】

题目分析：

以前一直不会这个方法， 我好菜啊。

转移分为三个部分，一个是直接成功，一个是转移到E1，还有一个是转移到自己周围的一圈儿点。 如果是叶子那么只能转移到父亲，如果不是叶子可以把非叶子的转移代换，这样也只转移到父亲，判一下无解就行了。

代码：

 #include<bits/stdc++.h>
 using namespace std;

 const int maxn = ;
 const double eps = 1e-;

 int n,Tmp,cas;
 double a[maxn],b[maxn];
 vector<int> g[maxn];

 double A[maxn],B[maxn],C[maxn];

 void dfs(int now,int f){
     int cnt = ;
     for(int i=;i<g[now].size();i++){
     if(g[now][i] == f) continue;
     cnt++; dfs(g[now][i],now);
     }
     if(cnt == ){A[now] = a[now]; B[now] = -a[now]-b[now]; C[now] = B[now];}
     else{
     double alpha = ,beta = ,gamma = ;
     for(int i=;i<g[now].size();i++){
         if(g[now][i] == f) continue;
         alpha += A[g[now][i]];beta += B[g[now][i]];gamma += C[g[now][i]];
     }
     if(f != ){
         gamma += cnt+;
         alpha *= (-a[now]-b[now])/(cnt+);
         beta *= (-a[now]-b[now])/(cnt+);
         gamma *= (-a[now]-b[now])/(cnt+);
         B[now] = (-a[now]-b[now])/(cnt+);
         beta = -beta; C[now] = gamma; A[now] = alpha+a[now];
         C[now] /= beta; A[now] /= beta; B[now] /= beta;
     }else{
         gamma += cnt;
         alpha *= (-a[now]-b[now])/cnt;
         beta *= (-a[now]-b[now])/cnt;
         gamma *= (-a[now]-b[now])/cnt;
         B[now] = ; C[now] = gamma; beta = -(alpha+beta+a[now]);
         if(beta <= eps){C[now] = -;}
         else C[now] /= beta;
     }
     }
 }

 void read(){
     scanf("%d",&n);
     memset(A,,sizeof(A));
     memset(B,,sizeof(B));
     memset(C,,sizeof(C));
     memset(a,,sizeof(a));
     memset(b,,sizeof(b));
     for(int i=;i<=n;i++) g[i].clear();
     for(int i=;i<n;i++){
     int u,v; scanf("%d%d",&u,&v);
     g[u].push_back(v); g[v].push_back(u);
     }
     for(int i=;i<=n;i++) {
     int u,v; scanf("%d%d",&u,&v);
     a[i] = u/100.0,b[i] = v/100.0;
     }
 }

 void work(){
     dfs(,);
     printf("Case %d: ",cas);
     if(C[] == -){puts("impossible");}
     else{printf("%.6lf\n",C[]);}
 }

 int main(){
     scanf("%d",&Tmp);
     while(cas++,Tmp--){
     read();
     work();
     }
     return ;
 }