观光奶牛(01规划+负环+二分)

题目

01分数规划问题
这类问题首先要写出表达式，然后将表达式整理位加减法一端>=0或者<=0，此类形式，然后转为图论问题，对图的边权进行调整，最后转为最短路、最小生成树、负环等问题。

求正环那么求最长路，让d一直增加，最后最长路边数最终会大于等于n，即可找到正环。

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define x first
#define y second
typedef pair<int,int> PII;
const int N=1e3+10,M=5e3+10;
int n,m,k,ans,ecnt;

int head[N],vis[N],cnt[N];
int f[N];
double d[N]; //别忘记
struct edge {
	int u,v,w,next;
} E[M];


void add(int u,int v,int w) {
	E[++ecnt].u=u;
	E[ecnt].v=v;
	E[ecnt].w=w;
	E[ecnt].next=head[u];
	head[u]=ecnt;
}



//spfa 
bool check(double mid) {
	memset(vis,0,sizeof vis);
	memset(cnt,0,sizeof cnt);
	memset(d,0,sizeof d);
	queue<int> q;
	for(int i=1; i<=n; i++) {
		q.push(i);
		vis[i]=1;
	}
	while(!q.empty()) {
		int u=q.front();
		q.pop();
		vis[u]=0;
		for(int i=head[u]; i; i=E[i].next) {
			int v=E[i].v;
			double w= f[u]-mid*E[i].w;//01规划推导出来的表达式
			if(d[v]<d[u]+w) { //寻找正环，那么求最长路
				d[v]=d[u]+w;
				cnt[v]=cnt[u]+1;//最长路边数
				if(cnt[v]>=n)return true;
				if(!vis[v]) {
					q.push(v);
					vis[v]=1;
				}
			}
		}
	}
	return false; 
}


int main() {

	cin>>n>>m;

	for(int i=1; i<=n; i++) {
		cin>>f[i];
	}

	for(int i=1; i<=m; i++) {
		int u,v,w;
		cin>>u>>v>>w;
		add(u,v,w);
	}

	double l=0,r=1000,esp=1e-4;//保留两位小数精确到-4即可 3位到-5 加二即可

	while(r-l>esp) {
		double mid=(l+r)/2;
		if(check(mid))l=mid;
		else r=mid;
	}

	printf("%.2f\n",r);



	return 0;
}