2018.09.12 poj3621Sightseeing Cows（01分数规划+spfa判环）

最新推荐文章于 2019-09-20 22:53:34 发布

SC.ldxcaicai

最新推荐文章于 2019-09-20 22:53:34 发布

阅读量163

点赞数 1

分类专栏： # 01分数规划 # spfa 文章标签： 01分数规划 spfa

本文链接：https://blog.csdn.net/dreaming__ldx/article/details/82667606

版权

图论专题同时被 3 个专栏收录

63 篇文章 1 订阅

订阅专栏

01分数规划

11 篇文章 0 订阅

订阅专栏

spfa

7 篇文章 0 订阅

订阅专栏

传送门
01分数规划板题啊。
发现就是一个最优比率环。
这个直接二分+spfa判负环就行了。
代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 1005
#define M 5005
using namespace std;
int n,m,first[N],cnt=0;
bool in[N];
double l,r,dis[N],d[N];
struct Node{int v,next;double w;}e[M<<1];
inline void add(int u,int v,double w){e[++cnt].v=v,e[cnt].next=first[u],e[cnt].w=w,first[u]=cnt;}
inline bool spfa(int p,double mid){
    in[p]=true;
    for(int i=first[p];i;i=e[i].next){
        int v=e[i].v;
        double tmp=e[i].w*mid-d[p];
        if(dis[v]>dis[p]+tmp){
            if(dis[v]=dis[p]+tmp,in[v])return true;
            if(spfa(v,mid))return true;
        }
    }
    return in[p]=false;
}
inline bool check(double tmp){
    memset(dis,0,sizeof(dis)),memset(in,false,sizeof(in));
    for(int i=1;i<=n;++i)if(spfa(i,tmp))return true;
    return false;
}
int main(){
    scanf("%d%d",&n,&m),r=1000.0,l=0.0;
    for(int i=1;i<=n;++i)scanf("%lf",&d[i]);
    for(int i=1;i<=m;++i){
        int u,v;
        double w;
        scanf("%d%d%lf",&u,&v,&w),add(u,v,w);
    }
    while(r-l>1e-4){
        double mid=(l+r)/2.0;
        if(check(mid))l=mid;
        else r=mid;
    }
    printf("%.2f",r);
    return 0;
}