217. 绿豆蛙的归宿 - 概率dp

给出一个有向无环的连通图,起点为 1,终点为 N,每条边都有一个长度。

数据保证从起点出发能够到达图中所有的点,图中所有的点也都能够到达终点。

绿豆蛙从起点出发,走向终点。

到达每一个顶点时,如果有 K 条离开该点的道路,绿豆蛙可以选择任意一条道路离开该点,并且走向每条路的概率为 1/K。

现在绿豆蛙想知道,从起点走到终点所经过的路径总长度的期望是多少?

输入格式

第一行: 两个整数 N,M,代表图中有 N 个点、M 条边。

第二行到第 1+M 行: 每行 3 个整数 a,b,c,代表从 a 到 b 有一条长度为 c 的有向边。

输出格式

输出从起点到终点路径总长度的期望值,结果四舍五入保留两位小数。

数据范围

1≤N≤105,

1≤M≤2N

输入样例:

复制代码
4 4
1 2 1
1 3 2
2 3 3
3 4 4

输出样例:

复制代码
7.00
cpp 复制代码
#include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
#define endl '\n'

using namespace std;

typedef pair<int, int> PII;
typedef long long ll;

const int N = 100010, M = 200010;

int h[N], e[M], w[M], ne[M], idx;
int dout[N], d[N];
double f[N];

void add(int a, int b, int c)
{
	e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx ++;
}

void dfs(int u)
{
	queue<int> q;
	q.push(u);
	
	while(q.size())
	{
		int t = q.front();
		q.pop();
		
		for(int i = h[t]; ~i; i = ne[i])
		{
			int j = e[i];
			f[j] += (w[i] + f[t]) / (double)dout[j];
			d[j] ++;
			if(d[j] == dout[j])q.push(j);
		}
	}
}

int main()
{
	IOS
	int n, m;
	cin >> n >> m;
	memset(h, -1, sizeof h);
	while(m --)
	{
		int a, b, c;
		cin >> a >> b >> c;
		add(b, a, c);
		dout[a] ++;
	}
	
	f[n] = 0;
	dfs(n);
	printf("%.2lf", f[1]);
	
	return 0;
}

E(aX + bY) = aE(X) + bE(Y);

每个点的期望是k个(下一个点的期望+到下一个点的路径长度)相加除k,从最后一个点开始算。

相关推荐
cpp_learner17 分钟前
C++ 实现责任链模式(Chain of Responsibility):从一堆 if-else 到可插拔的处理管道
c++·设计模式
residual_fan26 分钟前
深度渐进收缩学习(DPSL):面向强噪声与类内分散的机械复合故障诊断方法
人工智能·算法·数据挖掘·数据分析
一木 之林32 分钟前
五、C++ 新特性、关键字与编译原理(进阶)(一)
c语言·开发语言·c++
2601_9622982733 分钟前
交易开拓者通常使用什么编程语言?这种语言如何帮助自动化交易?
java·c++·python·编程语言·自动化交易
码匠许师傅39 分钟前
【设计模式精讲】23.备忘录模式(Memento)
c++·设计模式·软件工程·uml·备忘录模式
a187927218311 小时前
【算法】动态规划第五篇:区间 DP——从“看什么都像背包“到“最后戳谁“
算法·leetcode·动态规划·dp·区间·区间dp·算法讲解
键盘会跳舞1 小时前
【C++多线程】thread_local 深度剖析:线程局部存储的生命周期与实现
c++·多线程·thread_local
whitelbwwww1 小时前
c++ 多线程
开发语言·c++·算法
额鹅恶饿呃1 小时前
随着CentOS官方停服的时间越来越久,大量仍在使用CentOS7的企业和运维从业者
java·python·算法·c#·ruby
白狐_7981 小时前
408 数据结构|KMP算法核心:为什么不用主串回退
数据结构·算法