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;
}