01分数规划 观光奶牛

01分数规划的题目通常长这样

本质上是在求$$\frac{\sum f_i}{\sum t_i}$$的最大值，这类分数求最大、最小的问题被成为01分数规划问题。

01分数规划求最大值通常使用二分，二分的范围要看题目，对于本题 $$\frac{\sum 环点权}{\sum 环边权}$$ 的范围为$(0,1000]$

于是问题可以转化为问号处是否成立，$$\frac{\sum f_i}{\sum t_i}\stackrel{?}{>}mid\iff \sum \left(f_i-mid\cdot t_i\right)>0$$ 于是就是将所有边权变为点权 - mid * 边权，求是否有正环，而求环的存在性问题可以用spfa，这样问题就解决了

/*
 * @Author: 爱学习的图灵机
 * @Date: 2022-02-20 16:48:10
 * @LastEditTime: 2022-02-20 20:43:01
 * Bilibili:https://space.bilibili.com/7469540
 * 题目地址:https://www.acwing.com/problem/content/363/
 * @keywords:  spfa 01分数规划
 */
#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;
const int N = 1010,M = 5010;
int h[N], e[M], wt[M], ne[M], idx;
int q[N];
double dist[N];
bool st[N];
int wf[N];// 点权
int cnt[N]; 

int n, m;

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

bool check(double mid){ //check要执行多次，要初始化
//!double mid
    memset(dist, 0, sizeof dist);// 可以不要
    memset(st, 0, sizeof st);
    memset(cnt, 0, sizeof cnt);
    
    int hh = 0, tt = 0;
    for(int i = 1; i <= n; ++ i){
        st[i] = true;
        q[tt ++] = i;
    }
    
    while(hh != tt){
        int t = q[hh ++]; if(hh == N) hh = 0;
        st[t] = false;
        
        for(int i = h[t]; ~ i; i = ne[i]){
            int j = e[i];
            if(dist[j] < dist[t] + wf[t] - mid * wt[i]){
                dist[j] = dist[t] + wf[t] - mid * wt[i];
                cnt[j] = cnt[t] + 1;
                if(cnt[j] >= n)return true;
                
                if(!st[j]){
                    q[tt ++] = j;
                    if(tt == N) tt = 0;
                    st[j] = true;
                }
    
            }
        }
    }
    return false;
}

int main()
{
    memset(h, -1, sizeof h);
    cin >> n >> m;
    for(int i = 1; i <= n; ++ i){
        cin >> wf[i];
    }
    
    for(int i = 1; i <= m; ++ i){
        int a, b, c;
        cin >> 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;
}