poj 3621 Sightseeing Cows 二分+spfa判断负环

题目:https://vjudge.net/problem/POJ-3621

题意：

给定一张L个点、P条边的有向图，每个点都有一个权值f[i]，每条边都有一个权值t[i]。

求图中的一个环，使“环上各点的权值之和”除以“环上各边的权值之和”最大。

输出这个最大值。

 

一开始我是想着暴力求解每一个环的，然后想着用spfa好像比较难将不同的环分开，要么就改写stack形式的spfa，用一个pre数组记录路径？

然后就忽然想到了二分。。

要求这个东西的最大值

一般来说除法对求和的合并很不友好，变成乘法

然后

二分k，假如kval-poi能够形成负环，就lef=k ,否则 rig=k

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 1e3 + 7;
const int maxm = 5e3 + 7;
int n, m, tot;
int Head[maxn], Nxt[maxm << 1], To[maxm << 1], Val[maxm << 1]; 
int poi_Val[maxn]; //这个是点的值
void add_edge(int fro, int to, int val) {
	Nxt[++tot] = Head[fro];
	To[tot] = to;
	Val[tot] = val;
	Head[fro] = tot;
}
queue<int> Q;
int cnt[maxn];
double dis[maxn];
bool vis[maxn];
bool spfa(double mid) {
	while (!Q.empty())Q.pop();
	memset(cnt, 0, sizeof(cnt));
	memset(dis, 0, sizeof(dis));
	for (int i = 1; i <= n; i++) Q.push(i);
	fill(vis + 1, vis + n + 1, 1);
	cnt[1] = 1;
	while (!Q.empty()) {
		int now = Q.front();
		Q.pop();
		vis[now] = 0;
		for (int i = Head[now]; i; i = Nxt[i]) {
			const int &to = To[i], &val = Val[i], &poival = poi_Val[to];
			if (dis[to] > dis[now] + mid * val - poival) {
				dis[to] = dis[now] + mid * val - poival;
				if (!vis[to]) Q.push(to), vis[to] = 1;
				if (++cnt[to] == n) return true;
			}
		}
	}
	return false;
}
int main() {
	cin >> n >> m;
	for (int i = 1; i <= n; i++) scanf("%d", poi_Val + i);
	register int fro, to, val;
	for (int i = 1; i <= m; i++) {
		scanf("%d %d %d", &fro, &to, &val);
		add_edge(fro, to, val);
	}
	double l = 0, r = 1e6, mid;
	while (r - l > 1e-4) {
		mid = (l + r) / 2;
		if (spfa(mid)) l = mid;
		else r = mid;
	}
	printf("%.2lf\n", l); //这个要用c++交，假如用g++交的话不能用lf
	return 0;
}