Luogu 2868 [USACO07DEC]观光奶牛Sightseeing Cows

01分数规划复习。

这东西有一个名字叫做最优比率环。

首先这个答案具有单调性,我们考虑如何检验。

设$\frac{\sum_{i = 1}^{n}F_i}{\sum_{i = 1}^{n}T_i} = e$,我们需要检验的就是$\sum_{i = 1}^{n}(F_i - mid * T_i) \geq 0$是否存在。

感觉这玩意不好算,再变形一下:$\sum_{i = 1}^{n}(e * T_i - F_i) < 0$,就变成一个负环的检验了。

$F_i$应当可以任取一条有向边的入点和出点。

注意二分时的边界问题。

时间复杂度$O(logn (spfa???))$。

Code:

#include <cstdio>
#include <cstring>
using namespace std;
typedef long double db;

const int N = 1005;
const int M = 5005;
const db inf = 1e10;
const db eps = 1e-6;

int n, m, tot = 0, head[N];
db dis[N], a[N];
bool vis[N], ex;

struct Edge {
    int to, nxt;
    db val;
} e[M];

inline void add(int from, int to, db val) {
    e[++tot].to = to;
    e[tot].val = val;
    e[tot].nxt = head[from];
    head[from] = tot;
}

template <typename T>
inline void chkMax(T &x, T y) {
    if(y > x) x = y;
}

template <typename T>
inline void chkMin(T &x, T y) {
    if(y < x) x = y;
}

void dfs(int x, db mid) {
    if(ex) return;
    vis[x] = 1;
    for(int i = head[x]; i; i = e[i].nxt) {
        int y = e[i].to;
        if(dis[y] > dis[x] + mid * e[i].val - a[y]) {
            dis[y] = dis[x] + mid * e[i].val - a[y];
            if(vis[y]) {
                ex = 1;
                return;
            }
            dfs(y, mid);
        }
    }
    vis[x] = 0;
}

inline bool chk(db mid) {
    for(int i = 1; i <= n; i++) {
        dis[i] = 0.0;
        vis[i] = 0;
    }
    ex = 0;
    
    for(int i = 1; i <= n; i++) {
        dfs(i, mid);
        if(ex) break;
    }
    
    return ex;
}

int main() {
    scanf("%d%d", &n, &m);
    db ln = 0.0, rn = 0.0, mid, res;
    for(int i = 1; i <= n; i++) scanf("%Lf", &a[i]);
    for(int i = 1; i <= m; i++) {
        int x, y; db v;
        scanf("%d%d%Lf", &x, &y, &v);
        add(x, y, v);
        rn += v;
    }
    
    for(; ln + eps <= rn; ) {
        mid = (ln + rn) * 0.5;
        if(chk(mid)) ln = mid, res = mid;
        else rn = mid;
    }
    
    printf("%.2Lf\n", res);
    return 0;
}
View Code

 

转载于:https://www.cnblogs.com/CzxingcHen/p/9872495.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值