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

题意:
有n个房间,由n-1条隧道连通起来,实际上就形成了一棵树,
从结点1出发,开始走,在每个结点i都有3种可能:
1.被杀死,回到结点1处(概率为ki)
2.找到出口,走出迷宫 (概率为ei)
3.和该点相连有m条边,随机走一条
求:走出迷宫所要走的边数的期望值。

设 E[i]表示在结点i处,要走出迷宫所要走的边数的期望。E[1]即为所求。

叶子结点:
E[i] = ki*E[1] + ei*0 + (1-ki-ei)*(E[father[i]] + 1);
= ki*E[1] + (1-ki-ei)*E[father[i]] + (1-ki-ei);

非叶子结点:(m为与结点相连的边数)
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);

设对每个结点:E[i] = Ai*E[1] + Bi*E[father[i]] + Ci;

对于非叶子结点i,设j为i的孩子结点,则
∑(E[child[i]]) = ∑E[j]
= ∑(Aj*E[1] + Bj*E[father[j]] + Cj)
= ∑(Aj*E[1] + Bj*E[i] + Cj)
带入上面的式子得
(1 - (1-ki-ei)/m*∑Bj)*E[i] = (ki+(1-ki-ei)/m*∑Aj)*E[1] + (1-ki-ei)/m*E[father[i]] + (1-ki-ei) + (1-ki-ei)/m*∑Cj;
由此可得
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;

从叶子结点开始,直到算出 A1,B1,C1;

E[1] = A1*E[1] + B1*0 + C1;
所以
E[1] = C1 / (1 - A1);
若 A1趋近于1则无解...+

/*
HDU 4035 dp求期望的题。
题意:
有n个房间,由n-1条隧道连通起来,实际上就形成了一棵树,
从结点1出发,开始走,在每个结点i都有3种可能:
1.被杀死,回到结点1处(概率为ki)
2.找到出口,走出迷宫 (概率为ei)
3.和该点相连有m条边,随机走一条
求:走出迷宫所要走的边数的期望值。 设 E[i]表示在结点i处,要走出迷宫所要走的边数的期望。E[1]即为所求。 叶子结点:
E[i] = ki*E[1] + ei*0 + (1-ki-ei)*(E[father[i]] + 1);
= ki*E[1] + (1-ki-ei)*E[father[i]] + (1-ki-ei); 非叶子结点:(m为与结点相连的边数)
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); 设对每个结点:E[i] = Ai*E[1] + Bi*E[father[i]] + Ci; 对于非叶子结点i,设j为i的孩子结点,则
∑(E[child[i]]) = ∑E[j]
= ∑(Aj*E[1] + Bj*E[father[j]] + Cj)
= ∑(Aj*E[1] + Bj*E[i] + Cj)
带入上面的式子得
(1 - (1-ki-ei)/m*∑Bj)*E[i] = (ki+(1-ki-ei)/m*∑Aj)*E[1] + (1-ki-ei)/m*E[father[i]] + (1-ki-ei) + (1-ki-ei)/m*∑Cj;
由此可得
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; 从叶子结点开始,直到算出 A1,B1,C1; E[1] = A1*E[1] + B1*0 + C1;
所以
E[1] = C1 / (1 - A1);
若 A1趋近于1则无解... */
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<math.h>
#include<vector>
using namespace std;
const int MAXN=;
const double eps=1e-;//这里1e-8会WA。设为1e-9和1e-10可以
double k[MAXN],e[MAXN];
double A[MAXN],B[MAXN],C[MAXN]; vector<int>vec[MAXN];//存树 bool dfs(int t,int pre)//t的根结点是pre
{
int m=vec[t].size();//点t的度
A[t]=k[t];
B[t]=(-k[t]-e[t])/m;
C[t]=-k[t]-e[t];
double tmp=;
for(int i=;i<m;i++)
{
int v=vec[t][i];
if(v==pre)continue;
if(!dfs(v,t))return false;
A[t]+=(-k[t]-e[t])/m*A[v];
C[t]+=(-k[t]-e[t])/m*C[v];
tmp+=(-k[t]-e[t])/m*B[v];
}
if(fabs(tmp-)<eps)return false;
A[t]/=(-tmp);
B[t]/=(-tmp);
C[t]/=(-tmp);
return true;
}
int main()
{
int T;
int n;
int u,v;
int iCase=;
scanf("%d",&T);
while(T--)
{
iCase++;
scanf("%d",&n);
for(int i=;i<=n;i++)vec[i].clear();
for(int i=;i<n;i++)
{
scanf("%d%d",&u,&v);
vec[u].push_back(v);
vec[v].push_back(u);
}
for(int i=;i<=n;i++)
{
scanf("%lf%lf",&k[i],&e[i]);
k[i]/=;
e[i]/=;
}
printf("Case %d: ",iCase);
if(dfs(,-)&&fabs(-A[])>eps)
{
printf("%.6lf\n",C[]/(-A[]));
}
else printf("impossible\n");
}
}

最新文章

  1. Spark运行模式与Standalone模式部署
  2. libevent源码分析:listener
  3. [转]mvc3 使用session来存储类来存储用户登陆信息
  4. jquery 隐藏表单元素
  5. CSS选择器的权重与优先规则
  6. (十)stm32中FSMC的使用(用于LCD)
  7. 【和我一起学Python吧】Python3.0与2.X版本的区别
  8. Cannot connect to (local) sql server 2008
  9. c#使用XSLT将xml文档转换为html文档
  10. ORA-04021
  11. Ch06 验证
  12. OC 获取城市首字母
  13. Java 7之基础 - 强引用、弱引用、软引用、虚引用
  14. 关于vue-axios的post方式,后台无法解析传参问题
  15. 图文详解互联网根基之HTTP
  16. java8_api_集合框架
  17. centos7下zookeeper集群安装部署
  18. centos 7 mysql启动失败--学会看错误日志
  19. compact 创建一个包含变量名为数组的键和它们的值为数组的值的数组
  20. PHP json_encode函数中需要注意的地方

热门文章

  1. Log4Net使用详解(续)
  2. SpringBoot 使用logback
  3. 【转载】研发应该懂的binlog知识(下)
  4. Spring MVC-学习笔记(2)DispatcherServlet、@Controller、@RequestMapping、处理方法参数类型、可返回类型、Model、ModelMap、ModelAndView
  5. CodeForces.1174D.EhabandtheExpectedXORProblem(构造前缀异或和数组)
  6. 7、 正则化(Regularization)
  7. 使用 js 修饰器封装 axios
  8. js如何获取到select的option值???
  9. k8s结合helm部署
  10. vue.js(16)--vue的组件