http://acm.hdu.edu.cn/showproblem.php?pid=4035

题意:给一棵n个节点的树,每个节点有值k[i]和e[i],分别表示k[i]概率走向1号节点,e[i]概率获得胜利(即停止),如果没有进行上边任意操作,则等概率的走向与这个节点连边的点。问走过的边的期望。(n<=10000)

#include <cstdio>
#include <cstring>
using namespace std;
const int N=10005;
const double eps=1e-10;
struct dat { double a, b, c, K, E; }D[N];
struct E { int next, to; }e[N<<1];
int n, cnt, ihead[N];
void clr() {
memset(ihead, 0, sizeof(int)*(n+1));
cnt=0;
}
void add(int u, int v) {
e[++cnt].next=ihead[u]; ihead[u]=cnt; e[cnt].to=v;
e[++cnt].next=ihead[v]; ihead[v]=cnt; e[cnt].to=u;
}
double abs(double x) { return x<0?-x:x; }
void dfs(int x, int fa) {
bool flag=1;
double temp=1-D[x].K-D[x].E, suma=0, sumb=0, sumc=0;
int m=0;
for(int i=ihead[x]; i; ++m, i=e[i].next) if(e[i].to!=fa) {
flag=0;
dfs(e[i].to, x);
suma+=D[e[i].to].a;
sumb+=D[e[i].to].b;
sumc+=D[e[i].to].c;
}
if(flag) { D[x].a=D[x].K; D[x].b=D[x].c=temp; return; }
double di=m-temp*sumb;
D[x].a=(D[x].K*m+temp*suma)/di;
D[x].b=temp/di;
D[x].c=temp*(sumc+m)/di;
}
void work() {
dfs(1, -1);
if(abs(1-D[1].a)<=eps) puts("impossible");
else printf("%.6f\n", D[1].c/(1-D[1].a));
}
int main() {
int T; scanf("%d", &T);
for(int TT=1; TT<=T; ++TT) {
scanf("%d", &n);
for(int i=0; i<n-1; ++i) { int x, y; scanf("%d%d", &x, &y); add(x, y); }
for(int i=1; i<=n; ++i) scanf("%lf%lf", &D[i].K, &D[i].E), D[i].K/=100, D[i].E/=100;
printf("Case %d: ", TT);
work();
clr();
}
return 0;
}

  

好神的题...做出这题感觉对期望的认识大大加深了呢...

设状态$E[i]$表示当前在$i$节点还需要走$E[i]$的期望边即可结束,答案就是$E[1]$,容易得到:

对于叶子节点:

$$
\begin{align}
E[i] & = k[i]E[1]+e[i]*0+(1-k[i]-e[i])*(E[父亲]+1) \\
      & = k[i]E[1]+(1-k[i]-e[i])E[父亲]+(1-k[i]-e[i])
\end{align}
$$

对于非叶子节点:

$$
\begin{align}
E[i] & = k[i]E[1]+e[i]*0+\frac{(1-k[i]-e[i])(E[父亲]+1+\sum_{j是i孩子} (E[j]+1)}{m} \\
      & = k[i]E[1]+\frac{1-k[i]-e[i]}{m}E[父亲] + \frac{1-k[i]-e[i]}{m} \sum_{j是i孩子} E[j] + (1-k[i]-e[i])
\end{align}
$$

由于是树,我们要深度挖掘他们的性质:

1、我们发现非叶节点里只是多了$\sum_{j是i孩子} E[j]$

2、由于期望本来就是求极限意义,如果裸的这样求显然是递归无限的,我们需要通过数学技巧约掉式子

比如我们可以为了方便,设$E[i]$的一般式子为$E[i]=A[i]E[1]+B[i]E[父亲]+C[i]$(为什么带个父亲未知数呢?那是因为$i$的儿子需要知道自己父亲(即$i$)的期望,那么我们可以通过带入到原式中且移项即可解出$E[i]$,是不是很神奇呢0.0)

那么$\sum_{j是i孩子} E[j] = \sum_{j是i孩子} \left( A[j]E[1] + B[j]E[i] + C[j] \right)$

最终化简超强式子,变为求出这个一般方程的系数!:

对于叶子节点:

$$A[i]=k[i], B[i]=C[i]=1-k[i]-e[i]$$

对于非叶子节点:

设$di=m-(1-k[i]-e[i])\sum_{j} B[j], m=所有与i连边的节点数$

$$
\begin{align}
A[i] & = (k[i]m+(1-k[i]-e[i])\sum_{j} A[j]) / di \\
B[i] & = (1-k[i]-e[i]) / di \\
C[i] & = (1-k[i]-e[i])(m+\sum_{j} C[j]) / di
\end{align}
$$

最后答案:

$E[1]=C[i]/(1-A[1])$,判一下除数是否为0即可判断有无解。

【HDU】4035 Maze的更多相关文章

  1. 【HDU】4888 Redraw Beautiful Drawings 网络流【推断解是否唯一】

    传送门:pid=4888">[HDU]4888 Redraw Beautiful Drawings 题目分析: 比赛的时候看出是个网络流,可是没有敲出来.各种反面样例推倒自己(究其原因 ...

  2. 【HDU】2191 多重背包问题

    原题目:悼念512汶川大地震遇难同胞——珍惜现在,感恩生活 [算法]多重背包(有限背包) 动态规划 [题解]http://blog.csdn.net/acdreamers/article/detail ...

  3. 【HDU】6110 路径交(2017百度之星) 线段树+RMQ-LCA+树链的交

    [题目]2017"百度之星"程序设计大赛 - 初赛(A) [题意]给定n个点的带边权树,m条编号1~m的路径,Q次询问编号区间[L,R]所有链的交集的长度.n<=500000 ...

  4. 【HDU】6148 Valley Numer 数位DP

    [算法]数位DP [题意]定义V-number为从左到看单位数字未出现先递增后递减现象的数字,求0~N中满足条件的数字个数.T<=200,lenth(n)<=100 [题解]百度之星201 ...

  5. 【HDU】5269 ZYB loves Xor I

    [算法]trie [题解] 为了让数据有序,求lowbit无法直接排序,从而考虑倒过来排序,然后数据就会呈现出明显的规律: 法一:将数字倒着贴在字典树上,则容易发现两数的lowbit就是它们岔道结点的 ...

  6. 【HDU】3068 最长回文

    [算法]manacher [题解][算法]字符串 #include<cstdio> #include<algorithm> #include<cstring> us ...

  7. 【HDU】2222 Keywords Search

    [算法]AC自动机 [题解]本题注意题意是多少关键字能匹配而不是能匹配多少次,以及可能有重复单词. 询问时AC自动机与KMP最大的区别是因为建立了trie,所以对于目标串T与自动机串是否匹配只需要直接 ...

  8. 【HDU】6012 Lotus and Horticulture (BC#91 T2)

    [算法]离散化 [题解] 答案一定存在于区间的左右端点.与区间左右端点距离0.5的点上 于是把所有坐标扩大一倍,排序(即离散化). 让某个点的前缀和表示该点的答案. 初始sum=∑c[i] 在l[i] ...

  9. 【HDU】6146 Pokémon GO

    [题意]一个2*n的网格,再保证步数最少的情况下,求从任意格出发遍历完所有格的方案数,格子八连通.n<=10000,T<=100. [算法]递推,DP [题解]原题链接:蓝桥杯 格子刷油漆 ...

随机推荐

  1. Android屏幕旋转总结

    转自:http://www.myexception.cn/operating-system/1452058.html 1. ProjectConifg.mk中定义宏MTK_LCM_PHYSICAL_R ...

  2. Android Studio 配置

    Android配置:[转]原地址:http://www.cnblogs.com/smyhvae/p/4022844.html [开发环境] 物理机版本:Win7旗舰版(64位) Android Stu ...

  3. php array_intersect() 和 array_diff() 函数

    在PHP中,使用 array_intersect 求两个数组的交集比使用 array_diff 求同样两个数组的并集要快. 如果要求数组 $a 与数组 $b 的差集的个数,应该使用 count($a) ...

  4. jfreeChart柱状图各属性详细设置

    一. 下载与环境配置 此最新版本为 1.0.13 解压jfreechart-1.0.13.zip 将lib目录下的jfreechart-1.0.13.jar .jcommon-1.0.16.jar 复 ...

  5. jq获取后台json并解析

    参考: $(function () { $.ajax({ url: 'tsconfig.json', type: 'GET', dataType: 'json', timeout: 1000, cac ...

  6. C#复习、面向对象阶段开始

    C#复习:在控制台程序中使用结构体.集合,完成下列要求项目要求:一.连续输入5个学生的信息,每个学生都有以下4个内容:1.序号 - 根据输入的顺序自动生成,不需要手动填写,如输入第一个学生的序号是1, ...

  7. svn代码提交注意事项

    先全部add,查看是否有gen和bin,然后再忽略这两个,忽略时选择最后rescury....

  8. LoadRunner录制图片验证码

    LoadRunner录制图片验证码 LoadRunner自身是无法捕获到图片验证码的,但是我们可以帮助LoadRunner来实现验证码的捕获. 1.图片验证码 图片验证码的产生来自服务器端,由服务器生 ...

  9. position-relative 的问题

    对100%宽度的元素0001添加position-relative属性,如果再给left/right属性,可能会导致0001元素超出其父盒子的范围.如果盒子0001的父级元素是body,可能会出现滚动 ...

  10. 【转】】Android ADB命令大全

    ADB很强大,记住一些ADB命令有助于提高工作效率. 获取序列号: adb get-serialno 查看连接计算机的设备: adb devices 重启机器: adb reboot 重启到bootl ...