【BZOJ3470】Freda’s Walk

Description

雨后的Poetic Island空气格外清新，于是Freda和Rainbow出来散步。 Poetic Island的交通可以看作一张n个点、m 边的有向无环图。由于刚下过雨，每条边都有一个积水深度，而恰好Freda 和Rainbow都喜欢踩水玩儿，于是Ta们从某个点出发，选择走向哪条边的概率与该边的积水深度是成正比的。即：如果Freda和Rainbow现在在点u，点u 出发的所有边的积水深度之和为s，从u到v的边积水深度为w，那么Ta们选择走向v的概率就是 w/s。
Ta们会一直走下去，直到到达一个没有出边的点，那么散步的路程长度就是走过的边的数量。更特殊的是，Freda和Rainbow在出发之前还可以选择一条边，在散步过程中无视这条边的存在（当然也可以不选择）。请你帮忙计算一下，Ta 们从0号点出发，散步的路程长度的期望值最大是多少？

Input

第一行两个正整数 n、m。
接下来m行每行三个整数u、v、w，表示从u到v有一条无向边，积水深度为w。

Output

输出Freda和Rainbow散步的路程长度的最大期望值，四舍五入保留六位小数。

Sample Input

4 5
0 1 2
0 2 1
0 3 3
1 3 1
2 3 4

Sample Output

2.000000

HINT

对于 100% 的数据，2<=n<=10000，1<=m<=100000，0<=u,v<n，1<=w<=1000。

题解：由于是DAG，所以我们先反着跑拓扑排序，求出从每个点开始走的期望步数f[i]，再正着跑拓扑排序，求出从0号点走到这个点的概率p[i]。

然后枚举删除那条边<a,b>。首先删除这条边会使答案减去p[a]*(f[b]+1)，其次a的其他出边的概率都会增加。算一下贡献即可。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
using namespace std;
const int maxn=10010;
const int maxm=100010;
int n,m,cnt;
int to[maxm],next[maxm],val[maxm],head[maxn],d[maxn],s[maxn],pa[maxm],pb[maxm],pc[maxm];
double ans;
double p[maxn],f[maxn];
queue<int> q;
inline void add(int a,int b,int c)
{
	to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
}
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<'0'||gc>'9')	{if(gc=='-')	f=-f;	gc=getchar();}
	while(gc>='0'&&gc<='9')	ret=ret*10+(gc^'0'),gc=getchar();
	return ret*f;
}
int main()
{
	n=rd(),m=rd();
	int i,u,a,b,c;
	memset(head,-1,sizeof(head)),cnt=0;
	for(i=1;i<=m;i++)	a=pa[i]=rd()+1,b=pb[i]=rd()+1,c=pc[i]=rd(),add(b,a,c),s[a]+=c,d[a]++;
	for(i=1;i<=n;i++)	if(!d[i])	q.push(i);
	while(!q.empty())
	{
		u=q.front(),q.pop();
		for(i=head[u];i!=-1;i=next[i])
		{
			d[to[i]]--,f[to[i]]+=(f[u]+1)*val[i]/s[to[i]];
			if(!d[to[i]])	q.push(to[i]);
		}
	}
	memset(head,-1,sizeof(head)),cnt=0;
	for(i=1;i<=m;i++)	add(pa[i],pb[i],pc[i]),d[pb[i]]++;
	p[1]=1;
	for(i=1;i<=n;i++)	if(!d[i])	q.push(i);
	while(!q.empty())
	{
		u=q.front(),q.pop();
		for(i=head[u];i!=-1;i=next[i])
		{
			d[to[i]]--,p[to[i]]+=p[u]*val[i]/s[u];
			if(!d[to[i]])	q.push(to[i]);
		}
	}
	ans=f[1];
	for(i=1;i<=m;i++)
	{
		a=pa[i],b=pb[i],c=pc[i];
		double g=(f[a]-(f[b]+1)*c/s[a])*s[a]/(s[a]-c)-f[a];
		ans=max(ans,f[1]+g*p[a]);
	}
	printf("%.6lf",ans);
	return 0;
}