POJ 3621 Sightseeing Cows 最优比例环

题目链接：http://poj.org/problem?id=3621

 

题意：

       给你L个点，P条边，每个点和边都有各自的权值，现在要你求一个环，使得这个环的点权/边权最大。

 

做法：

       自己做当然是。。。不会做啦。。这个坑已经放在这里很久，之前做了最优比率生成树的，这样三个01分数规划的专题就算完成了。因为我们要求一个最大的,那么我们设这个比例为x，那么我们会有一个公式,所以我们要二分一个最优值x，假设由这个公式出来的值是正的，那么就说明还有一个更大的x值满足要求，那么就把l=mid继续二分，否则就把r=mid。可是如果我们要判断是否存在一个环满足大于0有点困难，因为spfa可以判断负环，那么我们就可以把公式改一下，令，如果存在一个负环，那么就是令l=mid；

       这里我在写之前一直有个疑问，不是说所有的地点只有第一次走的时候有效吗，那为什么每条路里面都要减掉happy[i]呢，想了一会儿，然后发现，因为每个点只能在一个环中，所以每次如果松弛的时候更新到这个点的话，那么就说明还有一个更符合要求的环需要这个点了，而如果我们要这个环，那么这个点就只会被遍历到一次，所以才要每一条边都要减掉happy[i]。

---

#include<queue>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;
const int maxn=3005;
const int inf=0x3f3f3f3f;
const double eps=1e-7;
int m,n,head[maxn],now,flag,vis[maxn],ct[maxn];
double dist[maxn],happy[maxn],w,maxhp,minhp,maxti,minti;
struct node{
	int to,next;
	double w;
}e[maxn*120];
void add(int u,int v,double w){
	e[now].to=v,e[now].next=head[u];
	e[now].w=w; head[u]=now++;
}
int spfa(int x,double mid){
	queue<int> q;
	dist[x]=0;
	vis[x]=1; ct[1]=1;
	q.push(x);
	while(!q.empty()){
		int u=q.front(); q.pop();
		if(++ct[u]>n) return 1;
		vis[u]=0;
		for(int i=head[u];~i;i=e[i].next){
			int t=e[i].to;
			double w=mid*e[i].w-happy[t];;
			if(dist[t]-w>dist[u]){
				dist[t]=w+dist[u];
				if(!vis[t]){
					vis[t]=1;
					q.push(t);
				}
			}
		}
	}
	return 0;
}
int ck(double mid){
    memset(ct,0,sizeof(ct));
    memset(vis,0,sizeof(vis));
    memset(dist,125,sizeof(dist));
    dist[1]=0;
    return spfa(1,mid);
}
int main(){
    int x,y;
    memset(head,-1,sizeof(head));
    maxhp=-1,maxti=-1,minhp=inf,minti=inf;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%lf",&happy[i]);
        maxhp=max(maxhp,happy[i]);
        minhp=min(minhp,happy[i]);
    }
    for(int i=0;i<m;i++){
        scanf("%d%d%lf",&x,&y,&w);
        add(x,y,w);
        maxti=max(maxti,w);
        minti=min(minti,w);
    }
    double l=minhp/maxti,r=maxhp/minti,mid;
    while(r-l>eps){
        mid=(l+r)/2;
        if(ck(mid)) l=mid;
        else r=mid;
    }
    printf("%.2f\n",mid);
    return 0;
}