poj-3621(01分数规划)

题目:Sightseeing Cows
思路:题目意思是值给出L个点和P条有向边,给个点都有点权fi和边权ti,要求从起点沿着某条路径转一圈再回到起点,求这个环上所有点权之和/所有边权之和的最大值。
在图论里,对于求像是在这里插入图片描述,上面一堆东西除下面一堆东西的最值,一般采用01分数规划。
求01分数规划一般用二分来求:

在这里插入图片描述

因此,可以用spfa来判断是否存在正环,来不断缩小答案所在的区间即可。
代码:

#include <iostream>
#include <cstring>

using namespace std;

const int N = 1010, M = 5210;

int n,m;
int wf[N];
int h[N],wt[M],e[M],ne[M],idx;
double dist[N];
int q[N],cnt[N];
bool st[N];

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){
	memset(cnt,0,sizeof cnt);
	memset(st,0,sizeof st);
	int hh = 0, tt = 0;
	for(int i = 1; i <= n; i++){
		q[tt++] = i;
		st[i] = true;
	}
	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(int argc, char** argv) {
	while(~scanf("%d%d",&n,&m)){
		memset(h,-1,sizeof h);
		idx = 0;
		for(int i = 1; i <= n; i++) scanf("%d",&wf[i]);
		for(int i = 0; i < m; i++){
			int a,b,c;
			scanf("%d%d%d",&a,&b,&c);
			add(a,b,c);
		}
		if(!check(0)) puts("0");
		else{
			
		}
		double l = 0, r = 1010;
		while(r - l > 1e-4){
			double mid = (l + r) / 2.0;
			if(check(mid)) l = mid;
			else r = mid;
		}
		printf("%.2lf\n",l);
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值