《算法竞赛进阶指南》观光奶牛

给定一张 L 个点、PP 条边的有向图，每个点都有一个权值 f[i]，每条边都有一个权值 t[i]。

求图中的一个环，使“环上各点的权值之和”除以“环上各边的权值之和”最大。

输出这个最大值。

注意：数据保证至少存在一个环。

输入格式

第一行包含两个整数 L 和 P。

接下来 L 行每行一个整数，表示 f[i]。

再接下来 P 行，每行三个整数 a，b，t[i]，表示点 a 和 b 之间存在一条边，边的权值为 t[i]。

输出格式

输出一个数表示结果，保留两位小数。

数据范围

2≤L≤1000,
2≤P≤5000,
1≤f[i],t[i]≤1000

输入样例：

5 7
30
10
10
5
10
1 2 3
2 3 2
3 4 5
3 5 2
4 5 5
5 1 3
5 2 2

输出样例：

6.00

假设:用sum_f，sum_t表示分别表示环内点权之和，边权之和：

一般像sum_f / sum_t这类题属于01分数规划问题，通用的方法是用二分解答

假设：mid为二分分界点，则令:

1:sum_f / sum_t > mid

2:变形后：sum_f > sum_t * mid

3：将右边移到左边sum_f - sum_t * mid > 0;（等价于求spfa中是否存在正环）

而我们要求的是最大的正环，所以mid的左半边表示spfa中不存在正环所以缩小范围r = mid;

右半边表示sum_f - sum_t * mid > 0成立，即spfa中存在正环所以扩大范围即l = mid;

        

注意：若要求存在正环

方法一：可以把边反向再求负环：

方法二：直接转化为求最长路判断是否存在正环

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

const int N = 1010, M = 5010;

int n, m;
int wf[N];
int cnt[N];
bool st[N];
double dist[N];
int h[N], wt[M], ne[M], e[M], idx;

void add(int a, int b, int c)
{
    e[idx] = b;
    wt[idx] = c;
    ne[idx] = h[a];
    h[a] = idx;
    idx ++ ;
}

bool check(double mid)
{
    memset(st, 0, sizeof st);
    memset(cnt, 0, sizeof cnt);
    
    queue<int> q;
    for (int i = 1; i <= n; i ++ )
    {
        st[i] = true;
        q.push(i);
    }
    
    while (q.size())
    {
        int t = q.front();
        q.pop();
        
        st[t] = false;
        
        for (int i = h[t]; i != -1; i = ne[i])
        {
            int j = e[i];
            
            if (dist[j] < dist[t] + wf[t] - mid * wt[i])//转化为最长路求正环
            {
                cnt[j] = cnt[t] + 1;
                if (cnt[j] >= n) return true;//说明存在正环
                
                dist[j] = dist[t] + wf[t] - mid * wt[i];
                if (!st[j])
                {
                    st[j] = true;
                    q.push(j);
                }
            }
        }
    }
    
    return false;
}

int main()
{
    cin >> n >> m;
    for (int i = 1; i <= n; i ++ ) scanf("%d", &wf[i]);
    
    memset(h, -1, sizeof h);
    
    while (m -- )
    {
        int a, b, c;
        scanf("%d %d %d", &a, &b, &c);
        add(a, b, c);
    }
    
    double l = 0, r = 1010;
    while (r - l > 1e-4)//保留两位小数，经验值：比答案多两位即可
    {
        double mid = (l + r) / 2;
        if (check(mid)) l = mid;//说明存在正环，扩大范围
        else r = mid;//不存在正环，缩小范围
    }
    
    printf("%.2lf", r);
    
    return 0;
}