Luogu P1768 天路 0/1分数规划+dfs spfa

“那是一条神奇的天路诶~~把第一个神犇送上天堂” 怕不是某大佬早就A了这题,然鹅我又调了很久很久。。。


好吧就是0/1分数规划,但是跑的dfs的spfa(好像题解说bfs过不了????不知) 发现把spfa写成bool的很难调。。。于是重构了一遍代码。。。

#include<cstdio>
#include<iostream>
#include<cstring>
#define R register int
using namespace std;
const int M=20010,N=7010;
const int eps=1E-2;
inline int g() {
    R ret=0; register char ch; while(!isdigit(ch=getchar())) ;
    do ret=ret*10+(ch^48); while(isdigit(ch=getchar())); return ret;
}
int n,m,cnt;
int vr[M],nxt[M],c[M],p[M],fir[N];
double d[N],md;
bool vis[N],flg;
inline void add(int u,int v,int cc,int pp) {vr[++cnt]=v,nxt[cnt]=fir[u],fir[u]=cnt,c[cnt]=cc,p[cnt]=pp;}
inline void spafa(int u) { vis[u]=true;
    for(R i=fir[u];i;i=nxt[i]) { R v=vr[i]; register double w=md*p[i]-c[i];
        if(d[v]>d[u]+w) { d[v]=d[u]+w;
            if(vis[v]) {flg=true; return ;}
            else spafa(v);    
        }
    } vis[u]=false;
}
signed main() {
    n=g(),m=g();
    for(R i=1,u,v,c,p;i<=m;++i) u=g(),v=g(),c=g(),p=g(),add(u,v,c,p);
    register double l=0,r=200;
    while(l+eps<r) { flg=false; memset(d,0,sizeof(d)); memset(vis,0,sizeof(vis));
        md=(l+r)/2; for(R i=1;i<=n;++i) if(!flg) spafa(i); else break; 
        if(flg) l=md+0.01; else r=md-0.01;
    } if(l) printf("%.1lf\n",l); else printf("-1\n");
}

2019.04.20 ...已经快5月了。。。

 

转载于:https://www.cnblogs.com/Jackpei/p/10740669.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值