HDU 4035 Maze 概率DP 搜索

解题报告链接：

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

先推公式，设计状态，令DP[i]表示在房间i退出要走步数的期望值，然后推导出关系式，亮点来了，不会搜索···

so，等我学会了搜索后明天再写····

哈哈，已经学会了，自己设计的，和解题报告的不太一样，一个DFS解决···

看代码吧···我不知道怎么说····

贴代码：

 //#define debug
 #include <cstdio>
 #include <vector>
 #include <cmath>
 #define N 10005
 #define eps 1e-10
 using namespace std;
 vector<int> ve[N];
 struct node
 {
     double ki,wal;
 } p[N];
 struct para
 {
     double a,b,c;
 };
 int n;
 bool vis[N];
 para dfs(int x)
 {
     vis[x] = ;
     bool flag = true;
     para t,t1,t2;
     //t1用来求sum(j),其中j为x的子节点
     t1.a = ;
     t1.b =;
     t1.c = ;
     for(int i=; i<ve[x].size(); ++i)
     {
         if(!vis[ve[x][i]])
         {
             flag =false;
             t2 = dfs(ve[x][i]);
             t1.a += t2.a;
             t1.b += t2.b;
             t1.c += t2.c;
         }
     }
 //      对于非叶子节点：j为i的子节点
 //      Ai =      (ki+(1-ki-ei)/m*∑Aj)  / (1 - (1-ki-ei)/m*∑Bj);
 //      Bi =        (1-ki-ei)/m  /   (1 - (1-ki-ei)/m*∑Bj);
 //      Ci = ( (1-ki-ei)+(1-ki-ei)/m*∑Cj ) / (1 - (1-ki-ei)/m*∑Bj);
 //      对于叶子节点：
 //      Ai = ki;
 //      Bi = 1 - ki - ei;
 //      Ci = 1 - ki - ei;
 //      wal = 1-ki-exi
     if(flag)//叶子节点
     {
         t.a = p[x].ki;
         t.b = p[x].wal;
         t.c = p[x].wal;
     }
     else//非叶子节点
     {
         int m = ve[x].size();
         double f =   - p[x].wal / m*t1.b;
         t.a =(p[x].ki + p[x].wal/m*t1.a) /f ;
         t.b = p[x].wal/m/f ;
         t.c = (p[x].wal + p[x].wal/m*t1.c) / f;
     }
     return t;
 }
 int main()
 {
 #ifdef debug
     freopen("in.c","r",stdin);
 #endif
     int t;
     scanf("%d",&t);
     for(int d =; d<=t ; ++d)
     {
         for(int i=; i<=n; ++i) ve[i].clear();
         scanf("%d",&n);
         for(int i=; i<n; ++i)
         {
             int u,v;
             scanf("%d%d",&u,&v);
             ve[u].push_back(v);
             ve[v].push_back(u);
         }
         for(int i=; i<=n; ++i)
         {
             int k,e;
             scanf("%d%d",&k,&e);
             p[i].ki = (double)k/100.0;
             p[i].wal =  - p[i].ki - (double)e/100.0;
         }
         memset(vis,,sizeof(vis));
         para ans = dfs();
         printf("Case %d: ",d);
         if(fabs(ans.a-) < eps)
             printf("impossible\n");
         else
             printf("%.6f\n",ans.c/(-ans.a));
     }
     return ;
 }