“蔚来杯“2022牛客暑期多校训练营2 D Link with Game Glitch

传送门：D-Link with Game Glitch_"蔚来杯"2022牛客暑期多校训练营2 (nowcoder.com)

看了官方题解以后深感难想，一个不算明确的图论题除了判负环还需要二分+取对数。在此记录一种相对简单的方法，参考补题时观摩大佬的代码：代码查看 (nowcoder.com)，本题解会简单说明该做法的一些正确性？（本蒟蒻没有能力严谨的证明呜呜呜）

还是先交代题意，可直接抽象为：n个点，m条边，边按a, b, c, d输入，b到d有向，权值为a/c；要求给所有边权加权重w，使得该图内所有简单环的权值之积小于等于1，求w最大值。

做法：

part1. 首先当然是按题解所说的取对数，不对的话无论是向上还是向下都可能溢出。

part2. 使用dp[i][u]表示走i步到u的权重之积，每一步操作我们都枚举每一条边，对于边{u, v, w}, dp[i][u] = min( dp[i][u], dp[i - 1][v] + w)，如此我们可以去除重边。当然不止于此，具体地说，相同步数i到达u的所有方法中我们只取权值之积最小的即可。因为我们只需求简单环（非简单环可以分解为简单环），所以最多走n步，即所有顶点相连。（哎~感觉讲的很烂，建议大家系统学习Bellman-Ford算法）

part3. 计算完dp[N][N]后，对于每个材料或者称顶点u，要计算以该店为首尾的简单环权值就很简单了：dp[n][u] - dp[i][u]，i为0到n-1（当然其中n-1显然不会构成环），表示步数或者更直接的我们称为深度。显然环的大小为n-i，所以我们枚举步数计算w：logw = (dp[n][u] - dp[i][u]) / (n - i)，因为之前我们已经做了很好的去重工作了，所以在此我们只要取max即可。

part4. 最后我们给每个顶点的简单环最大权值max{w1, w2, ... , wn-2}取min得到最终的logw，输出exp(ans)即可

#include<bits/stdc++.h>
using namespace std;
const int N = 1005;
const double inf = -1e8;
int n, m;
double dp[N][N];
struct node{
    int v,u;
    double w;
    node(int a, int b, double c){
        v = a, u = b, w = c;
    }
};
vector<struct node> v;
int main(){
    cin >> n >> m;
    for(int i = 1; i <= m;i ++){
        int b, d;
        double a, c;
        cin >> a >> b >> c >> d;
        node tmp(b, d, log(a/c));
        v.push_back(tmp);
    }
    // 所有简单环大小一定在n步之内
    for(int i = 1; i <= n; i ++){
        for(auto [v, u, w] : v){
            // dp[i][u] 表示第i步走到u的权值，取min操作对于重边等情况可以把权值较小的更新掉
//             cout << v << " " << u << " " << w << endl;
            dp[i][u] = min(dp[i][u], dp[i - 1][v] + w);
        }
    }
    double res = 0;
    for(int i = 1; i <= n; i ++){
        double ans = inf;
//         cout << ans << endl;
        for(int j = 0; j < n - 1;j ++){
            // 如何证明每种材料可以换取infinite的最小环消除以后就不会有材料infinite（局部->全局）
            ans = max(ans, (dp[n][i] - dp[j][i]) / (double)(n - j));
//             cout << exp(ans) << endl;
        }
        res = min(res , ans);
    }
    printf("%.8lf", exp(res));
    return 0;
}