题目地址http://acm.hdu.edu.cn/showproblem.php?pid=6714

题目:


Time Limit: 6000/4000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)

Problem Description

小 A 是社团里的工具人,有一天他的朋友给了他一个 n 个点,m 条边的正权连通无向图,要他计算所有点两两之间的最短路。

作为一个工具人,小 A 熟练掌握着 floyd 算法,设 w[i][j] 为原图中 (i,j) 之间的权值最小的边的权值,若没有边则 w[i][j]=无穷大。特别地,若 i=j,则 w[i][j]=0。

Floyd 的 C++ 实现如下:

```c++
for(int k=1;k<=p;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
  w[i][j]=min(w[i][j],w[i][k]+w[k][j]);
```

当 p=n 时,该代码就是我们所熟知的 floyd,然而小 A 为了让代码跑的更快点,所以想减少 p 的值。

令 Di,j 为最小的非负整数 x,满足当 p=x 时,点 i 与点 j 之间的最短路被正确计算了。

现在你需要求 ∑ni=1∑nj=1Di,j,虽然答案不会很大,但为了显得本题像个计数题,你还是需要将答案对 998244353 取模后输出。

Input

第一行一个正整数 T(T≤30) 表示数据组数

对于每组数据:

第一行两个正整数 n,m(1≤n≤1000,m≤2000),表示点数和边数。

保证最多只有 5 组数据满足 max(n,m)>200 

接下来 m 行,每行三个正整数 u,v,w 描述一条边权为 w 的边 (u,v),其中 1≤w≤109

Output

输出 T 行,第 i 行一个非负整数表示第 i 组数据的答案

Sample Input

1

4 4

1 2 1

2 3 1

3 4 1

4 1 1

Sample Output

6

 

解题思路:


(1)对于最短路直接是i->j的,=0,故我们不需要特殊考虑这种情况,直接把都初始化为0即可。

(2)对于≠0, floyd算法最外层k是从小到大的, 所以对于i-->j的路径上一直在增加节点k,随着k的增大,i-->j的最短路逐渐确定。所以k一定是i-->j最短路上除端点i,j之外的最大节点值。而i-->j的最短路有多条,那么对应的k值也就有多个,取最小的k值即为

最短路比较优的算法是spfa(O(m)的时间复杂度,m是边数),在寻找最短路的过程中做些处理即可,注意一定不包括端点i,j

 

ac代码:


嘿嘿!!

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 1005;
const ll inf = 1e18+10;
const int mod = 998244353;
struct node{
    int nxt, val;
};
vector<node> edges[maxn];
int vis[maxn], D[maxn];
ll dis[maxn];
int t, n, m, x, y, z;
queue<int> q;
void spfa(int s)
{
    for(int i = 1; i <= n; i++) dis[i] = inf;
    memset(vis, 0, sizeof vis); // 避免重复入队
    memset(D, 0, sizeof D);
    dis[s] = 0, vis[s] = 1;
    q.push(s);
    while(!q.empty())
    {
        int u = q.front(); q.pop();
        vis[u] = 0;
        for(int i = 0; i < edges[u].size(); i++)
        {
            int v = edges[u][i].nxt;
            ll w = edges[u][i].val;
            if(dis[v] > dis[u] + w)// s-->u->v,u!=v,但u可能等于s,就变成了s->v
            {
                dis[v] = dis[u] + w;
                if(u != s)D[v] = max(D[u], u);
                if(!vis[v])  vis[v] = 1, q.push(v);
            }
            else if(dis[v] == dis[u] + w)
            {
                int tmp = 0;
                if(u != s) tmp = max(D[u], u);
                if(D[v] > tmp) // 优化:多条最短路时只有当这条路的D【v】更小时才入队
                {
                    D[v] = tmp;
                    if(!vis[v]) vis[v] = 1, q.push(v);
                }
            }
        }
    }
}
int main()
{
    //freopen("/Users/zhangkanqi/Desktop/11.txt","r",stdin);
    scanf("%d", &t);
    while(t--)
    {
        scanf("%d %d", &n, &m);
        for(int i = 1; i <= n; i++)
            edges[i].clear();
        for(int i = 0; i < m; i++)
        {
            scanf("%d %d %d", &x, &y, &z);
            edges[x].push_back((node){y, z});
            edges[y].push_back((node){x, z});
        }
        ll ans = 0;
        for(int i = 1; i <= n; i++)
        {
            spfa(i);
            for(int j = 1; j <= n; j++)
                ans = (ans + D[j]) % mod;
        }
        printf("%lld\n", ans);
    }
    return 0;
}