【BZOJ1491】社交网络

题目链接：https://www.lydsy.com/JudgeOnline/problem.php?id=1491

---

 

算是Floyd的扩展吧，在求最短路的同时，记录最短路的条数。

一旦获得的任意两点间最短路的条数，就可以在O(n^3)的时间内求出每个点的l了。

提交一直出错，看了题解还是改了半天，最后才发现，题目保证任意两点间最短路条数小于10^10，但运算过程中却会爆，所以直接用double就好了。

具体的细节看代码就好了。

#include <cstdio>
#include <cstring>

const int maxn = 105, inf = 0x3f3f3f3f;

int dis[maxn][maxn];
double cnt[maxn][maxn];

int main() {
    int n, m, a, b, c;
    scanf("%d%d", &n, &m);
    memset(dis, inf, sizeof(dis));
    for (int i = 1; i <= n; ++i) dis[i][i] = 1, cnt[i][i] = 1;
    for (int i = 1; i <= m; ++i) {
        scanf("%d%d%d", &a, &b, &c);
        dis[a][b] = dis[b][a] = c;
        cnt[a][b] = cnt[b][a] = 1;
    }
    for (int k = 1; k <= n; ++k)
        for (int i = 1; i <= n; ++i)
            for (int j = 1; j <= n; ++j) {
                if (dis[i][j] > dis[i][k] + dis[k][j]) {
                    dis[i][j] = dis[i][k] + dis[k][j];
                    cnt[i][j] = cnt[i][k] * cnt[k][j];
                } else if (dis[i][j] == dis[i][k] + dis[k][j])
                    cnt[i][j] += cnt[i][k] * cnt[k][j];
            }
    for (int k = 1; k <= n; ++k) {
        double ans = 0;
        for (int i = 1; i <= n; ++i)
            for (int j = 1; j <= n; ++j) {
                if (dis[i][j] == dis[i][k] + dis[k][j])
                    ans += cnt[i][k] * cnt[k][j] / cnt[i][j];
            }
        printf("%.3f\n", ans);
    }
    return 0;
}