【BZOJ】1491: [NOI2007]社交網路（floyd）

weixin_34262482發表於2014-08-02

http://www.lydsy.com/JudgeOnline/problem.php?id=1491

囧囧囧。。。。。。。。。。。。。。。。。。。。。。

囧1：雖然自己想到做法了，但是操作的時候，竟然忘記之前想到的求最短路的數量。。。（充分證實我是蒟蒻，）然後稍稍看了下題解程式碼。。

囧2：初始化這個梗我就T_T...還有此題精度各種爆啊。。。後2wa就是精度有問題啊囧。這題要全部開double。。要不然就。。。。我查了下double的精度，omg。。300多位。。這太爽。

此題不用說了吧，，有點最短路基礎的都懂，求路徑數量的時候只要判斷d[i][j]==d[i][k]+d[k][j]時就根據乘法原理×上。

自環的最短路數量要=1，要不然，，你懂的。

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define read(a) a=getnum()
#define print(a) printf("%d", a)
inline int getnum() { int ret=0; char c; int k=1; for(c=getchar(); c<'0' || c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0' && c<='9'; c=getchar()) ret=ret*10+c-'0'; return k*ret; }

const int N=105;

double d[N][N], I[N], f[N][N];
int main() {
	int n=getnum(), m=getnum(), u, v;
	for1(i, 1, n) for1(j, 1, n) d[i][j]=1e15;
	for1(i, 1, m) read(u), read(v), scanf("%lf", &d[u][v]), d[v][u]=d[u][v], f[u][v]=f[v][u]=1;
	for1(k, 1, n) for1(i, 1, n) for1(j, 1, n) {
		if(d[i][j]>d[i][k]+d[k][j]) d[i][j]=d[i][k]+d[k][j], f[i][j]=0;
		if(d[i][j]==d[i][k]+d[k][j]) f[i][j]+=f[i][k]*f[k][j];
	}
	for1(i, 1, n) f[i][i]=0;
	for1(k, 1, n) 
		for1(i, 1, n) for1(j, 1, n) if(d[i][j]==d[i][k]+d[k][j] && f[i][j]>0)
			I[k]+=f[i][k]*f[k][j]/f[i][j];
	for1(i, 1, n) printf("%.3lf\n", I[i]);
	return 0;
}