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