BZOJ1491【floyd】

跑floyd的时候顺便统计一下最短路的方案数.

/* I will wait for you */

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iostream>
#include<fstream>
#include<vector>
#include<queue>
#include<deque>
#include<set>
#include<map>
#include<string>

typedef long long LL;
typedef unsigned long long ULL;

using namespace std;

const int maxn=110;
const int maxm=5010;
const int maxs=26;
const int INF=1<<29;
const int P=1000000007;
const double error=1e-9;

int n,m,cnt=0;
double wi[maxn],dis[maxn][maxn],f[maxn][maxn];
	
int main()
{
	scanf("%d%d",&n,&m);
	
	for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) f[i][j]=0.0,dis[i][j]=(double)INF;
	
	for(int i=0,u,v;i<m;i++) 
	{
		double w;scanf("%d%d%lf",&u,&v,&w);
		dis[u][v]=dis[v][u]=w;f[u][v]=f[v][u]=1.0;
	}
	
	for(int i=1;i<=n;i++) for(int s=1;s<=n;s++) for(int t=1;t<=n;t++)
	{
		if(dis[s][i]+dis[i][t]==dis[s][t]) f[s][t]+=f[s][i]*f[i][t];
		if(dis[s][i]+dis[i][t]<dis[s][t]) dis[s][t]=dis[s][i]+dis[i][t],f[s][t]=f[s][i]*f[i][t];
	}
	
	for(int i=1;i<=n;i++) f[i][i]=0;
	
	for(int i=1;i<=n;i++) for(int s=1;s<=n;s++) for(int t=1;t<=n;t++)
		if(f[s][t]&&dis[s][i]+dis[i][t]==dis[s][t]) 
			wi[i]+=(f[s][i]*f[i][t])/f[s][t];
			
	for(int i=1;i<=n;i++) printf("%.3lf\n",wi[i]);
	
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值