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] = 1;
    bool flag = true;
    para t,t1,t2;
    //t1用来求sum(j),其中j为x的子节点
    t1.a = 0;
    t1.b =0;
    t1.c = 0;
    for(int i=0; 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 =  1 - 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 =1; d<=t ; ++d)
    {
        for(int i=1; i<=n; ++i) ve[i].clear();
        scanf("%d",&n);
        for(int i=1; i<n; ++i)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            ve[u].push_back(v);
            ve[v].push_back(u);
        }
        for(int i=1; i<=n; ++i)
        {
            int k,e;
            scanf("%d%d",&k,&e);
            p[i].ki = (double)k/100.0;
            p[i].wal = 1 - p[i].ki - (double)e/100.0;
        }
        memset(vis,0,sizeof(vis));
        para ans = dfs(1);
        printf("Case %d: ",d);
        if(fabs(ans.a-1) < eps)
            printf("impossible\n");
        else
            printf("%.6f\n",ans.c/(1-ans.a));
    }
    return 0;
}