(LightOJ - 1321)Sending Packets(数学期望+最短路)

题目链接：Sending Packets - LightOJ 1321 - Virtual Judge

为了方便起见，我下面内容下标均从1开始

题意：有n个路由器，我们现在要从1号路由器向n号路由器传送数据，两个路由器之间的边代表数据沿这条边传输数据成功的概率。问传输完sKB的数据所需的最小期望时间，传输数据遵循以下规则：

（1）每次只能传输1KB的数据

（2）每次成功传输数据后我们在2*k时间后都能收到信息

（3）如果在传输当前数据2*k时间后我们没有收到信息，那我们就需要重新传输这条数据

分析：为什么期望时间还有最小一说呢？这是因为我们沿不同的路线传输数据成功的概率是不相同的，而每次传输数据成功的时间是固定的，所以按照不同的路径走求出来的时间期望是不同的，那按照什么样的路径求出来的时间期望才是最小的呢？显然按照传输数据成功率最高的路径所花费的时间期望是最小的。所以我们就需要先跑一边dijkstra求出起点到终点的最大概率。不妨设为p，则有   E=p*2K+（1-p）*(E+2K)，即成功期望是当前成功的概率乘以成功所花费的时间加上当前不成功的概率乘以当前不成功所要花费的时间期望。

化简得 E=2k/p，这是一次传输1KB数据的时间期望，最后答案就是传输1KB的时间期望乘以所要传输的数据量。

下面是代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<vector>
using namespace std;
const int N=20003;
typedef pair<double,int> PII;
double f[N],w[N];
int h[N],e[N],ne[N],idx;
bool vis[N];
void add(int x,int y,double z)
{
	e[idx]=y;
	w[idx]=z;
	ne[idx]=h[x];
	h[x]=idx++;
}
int T,n,m,s,k;
double dijkstra()
{
	memset(vis,false,sizeof vis);
	memset(f,0,sizeof f);
	f[1]=1;
	priority_queue<PII> q;
	q.push({f[1],1});
	while(q.size())
	{
		int begin=q.top().second;
		q.pop();
		if(vis[begin]) continue;
		vis[begin]=true;
		for(int i=h[begin];i!=-1;i=ne[i])
		{
			int j=e[i];
			if(f[j]<f[begin]*w[i])
			{
				f[j]=f[begin]*w[i];
				q.push({f[j],j});
			}
		}
	}
	return f[n];
}
int main()
{
	cin>>T;
	for(int o=1;o<=T;o++)
	{
		scanf("%d%d%d%d",&n,&m,&s,&k);
		memset(h,-1,sizeof h);
		idx=0;
		int u,v;
		double t;
		for(int i=1;i<=m;i++)
		{
			scanf("%d%d%lf",&u,&v,&t);
			add(u+1,v+1,0.01*t);add(v+1,u+1,0.01*t);//让下标从1开始 
		}
		double p=dijkstra();
		double ans=2*k/p;
		printf("Case %d: %lf\n",o,ans*s);
	}
	return 0;
}