ICode9

精准搜索请尝试: 精确搜索
首页 > 编程语言> 文章详细

算法学习---概率与期望

2021-08-30 10:33:07  阅读:125  来源: 互联网

标签:ch 期望 int MAX --- ++ 算法 dp


根据期望的线性性质,

(E[aX +bY]=aE[X] + bE[Y]) 。所以另外一种求期望的方式是分别求出每一种代价产生的期望贡献,然后相加得到答案。在本题中路径期望总长度等于每条边产生的期望贡献之和而每条边产生又等于经过这条边的期望次数乘这条边的代价所以,我们只需要算每条边的期望经过次数即可

((u,v,w)) 的期望经过次数是不好直接算的,但如果我们能算得点

(u) 的期望经过次数为

(dp(u,v)) 那么边

((u,v,w)) 的期望经过次数是

(dp(u)*\frac1{|G_u|}) 对答案的贡献就是

(w*dp(u)*\frac1{|G_u|})

如何计算点

(u) 的期望经过次数

(dp(u)) 呢?我们依然考虑 DP 的方式首先有

(dp(u) = 1) 转移采取刷表的方式

[dp(e_{to})\leftarrow dp(u)*\frac1{|G_u|},e\in G_u ]

在用边

(e) 刷表的同时

(e) 的贡献就可以计算了十分简洁。因为这种方法计算答案十分的便捷,而且适用范围广所以这种『利用期望的线性性质单独计算贡献的方法』是我们计算期望首选的方法

typedef long long ll; inline int readInt() { static int n, ch; n = 0, ch = getchar(); while (!isdigit(ch)) ch = getchar(); while (isdigit(ch)) n = n * 10 + ch - '0', ch = getchar(); return n; } const int MAX_N = 100000 + 3, MAX_M = MAX_N * 2; struct edge { edge *next; int to, cost; edge(edge *next = NULL, int to = 0, int cost = 0): next(next), to(to), cost(cost) {} } pool[MAX_M], *pit = pool, *first[MAX_N]; int n, m, deg[MAX_N], outDeg[MAX_N]; double f[MAX_N]; void solve() { static int q[MAX_N]; int l = 0, r = 0; q[r++] = 0; double ans = 0; f[0] = 1.0; while (r - l >= 1) { int u = q[l++]; for (edge *e = first[u]; e; e = e->next) { f[e->to] += f[u] / outDeg[u]; ans += f[u] * e->cost / outDeg[u]; if (--deg[e->to] == 0) q[r++] = e->to; } } printf("%.2f\n", ans); } int main() { n = readInt(), m = readInt(); for (int i = 0, u, v, w; i < m; ++i) { u = readInt() - 1, v = readInt() - 1, w = readInt(); first[u] = new (pit++) edge(first[u], v, w); deg[v]++, outDeg[u]++; } solve(); }

标签:ch,期望,int,MAX,---,++,算法,dp
来源: https://www.cnblogs.com/xiaohe0012/p/15204584.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有