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

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

题解：

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

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

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

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

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

但是这样需要开$O(n^3)$的空间，这题$n\le 500$，空间过大（PS：这题$10$组$500^3$简直有毒）。此时不难意识到，这个$dis[k][k][k]$的$k$从大到小枚举，和从小到大实际上是一样的道理，可以在$Floyd$的过程中统计答案，优化掉$k$这一维的空间。$Floyd$算法的某个$k$的迭代结束后，更新最小环大小，若$dis[k][k]=mn$则数量加上$cnt[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();
    }
}