poj3621 最优比例生成环

//给出一个有向图，问求一个回路，使得回路上的点权之和/边权之和 最大
//01分数规划
//首先可以证明最优解一定是一个单独的环。
//设大环套小环最优，大环为花费x1，收益y1，小环为花费x2，收益y2，重叠部分花费x3，收益y3
//那么大环套小环花费x1+x2,收益y1+y2-y3,可以推出（y1+y2-y3）/(x1+x2)要么小于y1/x1，要么小于y2/x2
//则最优解一定是一个单独环。
//一个环有n个点，n条边，我们可以把起点点权合并到边上。
//我们想求的是这样一个函数 f = max{ sigma(v[i]*x[i])/sigmax(e[i]*x[i]) } x[i]属于{0，1}
//那么设 g(y) = max{sigma(v[i]*x[i])-y*sigmax(e[i]*x[i]}
//此函数有两个性质 1.非严格单调递减   2.仅当y为最优解y*的时候，g(y)=0;
//当g(y)=0时，y=y*；
//当g(y)>0时, y<y*;
//当g(y)<0时, y>y*;
//所以只要图中存在正环，则g(y)>0
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace::std;
const int N=1010,M=100100;
double F[N],dis[N];//一定为double
struct Edge{
    int to,next;
    double cap;
}edge[M];
int cnt=0,head[N];
void addedge(int u,int v,double w){
    edge[cnt].to=v;edge[cnt].cap=w;edge[cnt].next=head[u];head[u]=cnt++;
}
int vis[N],c[N];    int Q[M*10];
bool spfa(double mid,int n){
    int front,rear;
    front=rear=0;
    memset(vis, 0, sizeof(vis));
    for(int i=1;i<=n;i++){
        vis[i]=c[i]=1;
        Q[rear++]=i;
        dis[i]=0;
    }
    while(front!=rear){
        int top=Q[front++];vis[top]=0;
        for(int i=head[top];i!=-1;i=edge[i].next){
            int v=edge[i].to;
            double w=F[top]-mid*edge[i].cap;
            if(dis[v]<dis[top]+w){
                dis[v]=dis[top]+w;
                if(!vis[v]){
                    Q[rear++]=v;
                    vis[v]=1;
                    c[v]++;
                    if(c[v]>n)return 1;
                }
            }
        }
    }
    return 0;
}
int main(){
    cnt=0;
    memset(head, -1, sizeof(head));
    int L,P;
    scanf("%d%d",&L,&P);
    for(int i=1;i<=L;i++)
        scanf("%lf",&F[i]);
    for(int i=0;i<P;i++){
        int u,v;double w;
        scanf("%d%d%lf",&u,&v,&w);
        addedge(u,v,w);
    }
    double l=0,r=10000;
    double eps=1e-6;
   while(r-l>eps){
        double mid=(l+r)/2;
        if(spfa(mid,L)){//有正环
            l=mid;//增大l
        }
        else{
            r=mid;
        }
    }
    printf("%.2lf\n",l);
    return 0;
}