[HNOI2009] 最小圈

[HNOI2009] 最小圈

题目描述

考虑带权有向图 $G=(V,E)$ 以及 $w:E\to \mathbb{R}$，每条边 $e=(i,j)$（$i\ne j$，$i,j\in V$）的权值定义为 $w_{i,j}$。设 $n=|V|$。

$c=(c_1,c_2,\cdots ,c_k)$（$c_i\in V$）是 $G$ 中的一个圈当且仅当 $(c_i,c_{i+1})$（$1\le i<k$）和 $(c_k,c_1)$ 都在 $E$ 中。称 $k$ 为圈 $c$ 的长度，同时记 $c_{k+1}=c_1$，并定义圈 $c=(c_1,c_2,\cdots ,c_k)$ 的平均值为
$$\mu (c)=\frac{1}{k}\sum _{i=1}^k{w}_{c_i,c_{i+1}}$$
即 $c$ 上所有边的权值的平均值。设 ${\mu }^{\prime }(G)={\min }_c\mu (c)$ 为 $G$ 中所有圈 $c$ 的平均值的最小值。

给定图 $G=(V,E)$ 以及 $w:E\to \mathbb{R}$，求出 $G$ 中所有圈 $c$ 的平均值的最小值 ${\mu }^{\prime }(G)$。

输入格式

第一行两个正整数，分别为 $n$ 和 $m$，并用一个空格隔开。其中 $n=|V|$，$m=|E|$ 分别表示图中有 $n$ 个点 和 $m$ 条边。

接下来 $m$ 行，每行三个数 $i,j,w_{i,j}$，表示有一条边 $(i,j)$ 且该边的权值为 $w_{i,j}$，注意边权可以是实数。输入数据保证图 $G=(V,E)$ 连通，存在圈且有一个点能到达其他所有点。

输出格式

一个实数 ${\mu }^{\prime }(G)$，要求精确到小数点后 $8$ 位。

样例 #1

样例输入 #1

4 5
1 2 5
2 3 5
3 1 5
2 4 3
4 1 3

样例输出 #1

3.66666667

样例 #2

样例输入 #2

2 2
1 2 -2.9
2 1 -3.1

样例输出 #2

-3.00000000

提示

对于 $100\%$ 的数据，$2\le n\le 3000$，$1\le m\le 10000$，$|w_{i,j}|\le 10^7$，$1\le i,j\le n$ 且 $i\ne j$。

---

提示：本题存在 $O(nm)$ 的做法，但是 $O(nm\log n)$ 的做法也可以通过。

题意

我们需要在一张有向图上找到一个"平均权值最小的一个圈"。

分析

我们观察圈的平均权值的定义
$$\mu (c)=\frac{1}{k}\sum _{i=1}^k{w}_{c_i,c_{i+1}}$$
处理一下
$$\mu (c)=\sum _{i=1}^n\frac{W_i}{k}$$
这种形式的分式，常见的做法（我知道的做法 ）就是分式规划了。
我们设:
$$\sum _{i=1}^k\frac{W_i}{k}\le X$$
处理一下
$$\sum _{i=1}^k(W_i-X)\le 0$$
意味着我们以 $W_i-X$ 为 权值建图，在图中如果找到了负环则说明上式成立。我们知道SPFA算法可以寻找负环（简单来说，如果一个点到起点的距离被重复更新，说明图中存在负环）。

时间复杂度: 我们需要跑整个图 $NM$的复杂度，分数规划二分复杂度$log(W\ast W)$， 总复杂度:
$$O(NMlog(W\ast W))$$

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
//#define int long long
#define double long double
#define ios ios::sync_with_stdio(false); cin.tie(0);cout.tie(0)
const int N = 1e5+5;
const double eps = 1e-10;
const int INF = 0x3f3f3f3f;
typedef pair<int, double > PID;
int n, m, vis[N];
double dis[N];
vector<PID> e[N];

bool spfa(int u, double x){
    vis[u] = 1;
    for(auto [v, w] : e[u]){
        if(dis[v] > dis[u] + w - x){
            dis[v] = dis[u] + w - x;
            if(vis[v] || spfa(v, x)) return 1;
        }
    }
    vis[u] = 0;
    return 0;
}

bool check(double x){
    memset(vis, 0, sizeof(vis));
    for(int i = 1; i <= n; i++) dis[i] = INF;
    for(int i = 1; i <= n; i++){
        if(spfa(i, x)) return 1;
    }
    return 0;
}

void solve(){
    cin >> n >> m;
    int u, v;
    double w;
    for(int i = 1; i <= m; i++){
        cin >> u >> v >> w;
        e[u].push_back({v, w});
    }
    double l = -1e7, r = 1e7; 
    while(r - l > eps){
        double mid = (l + r)/2;
        if(check(mid)) r = mid;
        else l = mid;
    }
    //cout << fixed << setprecision(10) << l <<" " << r << endl;
    cout << fixed << setprecision(8) << (r) << endl;
}

signed main(){
    ios; int _;
    _ = 1; //cin >> _;
    while(_--){
        solve();
    }
    return 0;
}