bzoj3036 绿豆蛙的归宿

3036: 绿豆蛙的归宿

Time Limit: 2 Sec   Memory Limit: 128 MB
Submit: 319   Solved: 224
[ Submit][ Status][ Discuss]

Description

随着新版百度空间的下线,Blog宠物绿豆蛙完成了它的使命,去寻找它新的归宿。

给出一个有向无环的连通图,起点为1终点为N,每条边都有一个长度。绿豆蛙从起点出发,走向终点。
到达每一个顶点时,如果有K条离开该点的道路,绿豆蛙可以选择任意一条道路离开该点,并且走向每条路的概率为 1/K 。
现在绿豆蛙想知道,从起点走到终点的所经过的路径总长度期望是多少?

Input

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

Output


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


Sample Input

4 4
1 2 1
1 3 2
2 3 3
3 4 4

Sample Output

7.00

HINT



对于100%的数据  N<=100000,M<=2*N

Source




概率期望DP

f[i]表示从i到n的期望长度,则状态转移方程为:

f[i]=∑(f[e[i].to]+e[i].v)/outd[i]

其中outd[i]表示i点的出度,e[i].to是与i相连的点。

一个小的总结:如果一个点的每条出边概率相同,则逆向DP更加容易求解




#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
#define ll long long
#define maxn 100005
using namespace std;
int n,m,cnt,head[maxn];
double f[maxn];
bool vst[maxn];
struct edge_type{int next,to,v;}e[maxn*2];
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
inline void add_edge(int x,int y,int v)
{
	e[++cnt]=(edge_type){head[x],y,v};
	head[x]=cnt;
}
inline void dfs(int x)
{
	if (vst[x]) return;
	vst[x]=true;
	int outd=0;
	for(int i=head[x];i;i=e[i].next)
	{
		dfs(e[i].to);
		f[x]+=f[e[i].to]+e[i].v;
		outd++;
	}
	if (!outd) return;
	f[x]/=outd;
}
int main()
{
	n=read();m=read();
	F(i,1,m)
	{
		int x=read(),y=read(),z=read();
		add_edge(x,y,z);
	}
	dfs(1);
	printf("%.2lf\n",f[1]);
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值