HDOJ 7322 Circuit 最小环计数 —— 2023“钉耙编程”中国大学生算法设计超级联赛（4）（2023杭电多校第四场）

原创已于 2023-07-28 18:31:11 修改 · 423 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#算法 #图论 #c++

于 2023-07-28 15:45:33 首次发布

ACM题目专栏收录该内容

3 篇文章

订阅专栏

文章讲述了使用Floyd最短路算法来解决图中的最小环问题，包括如何计算环的长度和数量，以及处理过程中涉及到的动态规划和空间优化策略。最后，提供了具体的代码实现来演示算法的运用。

题意： $T(1≤T≤15)T(1\leq T\leq 15)$ 组测试，每组测试给 $n$ 个点 $m$ 条有向边 $(u→v,w)(u\to v,w)$ , $(1≤n≤500,0≤m≤n(n−1),1≤u,v≤n,1≤w≤109)(1\leq n \leq 500, 0\leq m \leq n(n-1), 1\leq u,v\leq n, 1\leq w\leq 10^9)$ ，（并未保证无重边、自环），你需要求最小环的长度和数量 $(mod998244353)\pmod {998244353}$ 。保证 $∑n≤5000\sum n \leq 5000$ 。

题解：

要求最小环，不难想到 $Fl oy d$ 最短路。用 $d i s$ 表示距离矩阵，将 $d i s [i] [j]$ 全部初始化为 $in f$ （不要将 $d i s [i] [i]$ 初始化为 $0$ ），将边存入 $d i s$ ，跑完 $Fl oy d$ ，得到的 $d i s [i] [i]$ 即为以 $i$ 为起点的最小环。
还要计数，设 $c n t [i] [j] = i$ 到 $j$ 的最短路的数量，在跑 $Fl oy d$ 的过程中，若 $d i s [i] [k] + d i s [k] [j] < d i s [i] [j]$ ，则根据乘法原理，i到j的最短路数量有 $c n t [i] [k] * c n t [k] [j]$ 条

跑完Floyd后求出最小环长度 $mn=mini=1ndis[i][i]mn=\mathop{\text{min}}\limits_{i=1}^n dis[i][i]$ ，再求最小环数量 $ans=∑i∣dis[i][i]=mncnt[i][i]ans=\sum\limits_{i|dis[i][i]=mn}cnt[i][i]$ ，此时就会发现数量算多了，因为根据当前这个方法，对于一个 $x$ 元环，这 $x$ 个点都会把它算上一次。

现在需要想办法，让每个环只被一个起点计数一次，而不是被每个起点计数。不难想到一种删点的方法：求出以 $i$ 为起点的最小环数量，删除 $i$ 点，不断循环直到点被删完。

这个做法看似复杂度过大，实则在 $Fl oy d$ 的过程中已经得到。我们再仔细想一下 $Fl oy d$ 算法的动态规划原理：
状态： $d i s [k] [i] [j] = i$ 到 $j$ ，除了 $ij$ 以外只经过了 $1∼k1\sim k$ 的最短路
初态： $d i s [0] [i] [j] = i$ 到 $j$ 的边权
状态转移方程： $d i s [k] [i] [j] = min (d i s [k - 1] [i] [j], d i s [k - 1] [i] [k] + d i s [k - 1] [k] [j])$
空间优化：随着 $k$ 增大， $d i s [k] [i] [j]$ 单调递减，可以优化掉第一个维度 $k$ 。

假如我们现在记录了 $d i s [k] [i] [j] = i$ 到 $j$ 除了 $ij$ 以外只经过了 $1∼k1\sim k$ 的最短路， $c n t [k] [i] [j] = i$ 到 $j$ 除了 $ij$ 以外只经过了 $1∼k1\sim k$ 的最短路数量，那么从大到小枚举 $k$ ，若 $d i s [k] [k] [k]$ 等于最小环，则答案加上 $c n t [k] [k] [k]$ ，当 $k$ 变为 $k - 1$ 时，相当于删除了第 $k$ 个点，再询问以第 $k - 1$ 个点为起点的最小环。

但是这样需要开 $O(n^3)$ 的空间，这题 $n≤500n\leq 500$ ，空间过大（PS：这题 $10$ 组 $500^3$ 简直有毒）。此时不难意识到，这个 $d i s [k] [k] [k]$ 的 $k$ 从大到小枚举，和从小到大实际上是一样的道理，可以在 $Fl oy d$ 的过程中统计答案，优化掉 $k$ 这一维的空间。 $Fl oy d$ 算法的某个 $k$ 的迭代结束后，更新最小环大小，若 $d i s [k] [k] = mn$ 则数量加上 $c n t [k] [k]$ ，若最小环大小被更新，则之前累计的数量清空。

#include<bits/stdc++.h>

using namespace std;
using ll = long long;
const int N = 505, mod = 998244353;
const ll inf = 1e15;
int n, m;
ll dis[N][N], cnt[N][N];

void work() {
    cin >> n >> m;
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= n; ++j) {
            dis[i][j] = inf;
            cnt[i][j] = 0;
        }
    }
    for (int i = 1; i <= m; ++i) {
        int x, y;
        ll z;
        cin >> x >> y >> z;
        if (z < dis[x][y]) {
            dis[x][y] = z;
            cnt[x][y] = 1;
        }
        else if (z == dis[x][y]) {
            ++cnt[x][y];
        }
    }
    ll mn = inf, ans = 0;
    for (int k = 1; k <= n; ++k) {
        for (int i = 1; i <= n; ++i) {
            for (int j = 1; j <= n; ++j) {
                ll t = dis[i][k] + dis[k][j];
                if (t < dis[i][j]) {
                    dis[i][j] = t;
                    cnt[i][j] = cnt[i][k] * cnt[k][j] % mod;
                }
                else if (t == dis[i][j]) {
                    cnt[i][j] += cnt[i][k] * cnt[k][j];
                    cnt[i][j] %= mod;
                }
            }
        }
        for (int i = 1; i <= n; ++i) {
            if (dis[i][i] < mn) {
                mn = dis[i][i];
                ans = 0;
            }
        }
        if (dis[k][k] == mn) {
            ans += cnt[k][k];
        }
    }
    ans = (ans % mod + mod) % mod;
    if (mn == inf) {
        cout << "-1 -1\n";
    }
    else {
        cout << mn << " " << ans << "\n";
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    int T = 1;
    cin >> T;
    while (T--) {
        work();
    }
}