poj3621 Sightseeing Cows（最优比率环【学英语）

最新推荐文章于 2021-08-04 15:15:41 发布

Coco_T_

最新推荐文章于 2021-08-04 15:15:41 发布

阅读量277

点赞数

分类专栏： 01分数规划图论文章标签： 01分数规划

本文链接：https://blog.csdn.net/wu_tongtong/article/details/79437652

版权

图论同时被 2 个专栏收录

64 篇文章 0 订阅

订阅专栏

01分数规划

7 篇文章 0 订阅

订阅专栏

题目链接

分析：
01分数规划：最优比率环

求 $∑val[i] \over ∑e[i]$
其中 $val$ 是环中点的价值， $e$ 是环中边的价值

设答案为 $ans$
则对于所有环都有： ${∑val[i] \over ∑e[i]}<=ans$
$ans∗∑e[i]−∑val[i]>=0$
对于答案 $k$ ：
当 $k<ans$ 时，就存在至少一个环 $k∗∑e[i]−∑val[i]<0$ ，即有负权回路（边权为 $k∗e[i]−val[i]$ ）
当 $k>=ans$ 时，就对于所有的环 $k∗∑e[i]−∑val[i]>=0$ ，即没有负权回路

二分答案，把边权设为 $mid*w(x,y)-val[x]$ ，判断是否有负环
有则放大范围，无则缩小范围

tip

$dis[y]=min(dis[x]+mid*w(x,y)-val[x])$

还是要分清楚这几个单词的意思：
这里写图片描述

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

using namespace std;

const double eps=1e-6;
const int N=5002;
struct node{
    int y,nxt;
    double v;
};
node way[N];
int n,m,st[N],tot=0,cnt[N];
double val[N],dis[N];
bool in[N];

void add(int u,int w,double z) {
    tot++;way[tot].y=w;way[tot].v=z;way[tot].nxt=st[u];st[u]=tot;
}

int spfa(double x) {
    queue<int> Q;

    for (int i=1;i<=n;i++) {    //多源
        in[i]=1; dis[i]=0; Q.push(i); cnt[i]=0;
    }
    while (!Q.empty()) {
        int now=Q.front(); Q.pop();
        in[now]=0;
        for (int i=st[now];i;i=way[i].nxt) {
            int y=way[i].y;
            if (dis[now]+x*way[i].v-val[now]<dis[y]) {
                dis[y]=dis[now]+x*way[i].v-val[now];
                if (!in[y]) {
                    in[y]=1;
                    Q.push(y);
                    if (++cnt[y]>n) return 1;
                }
            }
        }
    }
    return 0;
}

int main()
{
    double l=0,r=0,ans=0;
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%lf",&val[i]),r=max(r,val[i]);
    for (int i=1;i<=m;i++) {
        int u,w; double z;
        scanf("%d%d%lf",&u,&w,&z);
        add(u,w,z);
    }

    while (r-l>=eps) {
        double mid=(l+r)/2;
        if (spfa(mid)) ans=mid,l=mid;
        else r=mid;
    }
    printf("%0.2lf\n",ans);

    return 0;
}