最小比例生成树

小澳的葫芦
【 题目描述】
小澳最喜欢的歌曲就是《葫芦娃》。
一日表演唱歌，他尽了洪荒之力，唱响心中圣歌。
随之，小澳进入了葫芦世界。
葫芦世界有 n 个葫芦，标号为 1~ n。 n 个葫芦由 m 条藤连接， 每条藤连接了
两个葫芦， 这些藤构成了一张有向无环图。小澳爬过每条藤都会消耗一定的能量。
小澳站在 1 号葫芦上（你可以认为葫芦非常大，可以承受小澳的体重），他
想沿着藤爬到 n 号葫芦上， 其中每个葫芦只经过一次。
小澳找到一条路径，使得消耗的能量与经过的葫芦数的比值最小。
【 输入格式】
输入文件名为 calabash.in。
输入文件第一行两个正整数 n,m， 分别表示葫芦的个数和藤数。
接下来 m 行，每行三个正整数 u,v,w，描述一条藤，表示这条藤由 u 连向 v，
小澳爬过这条藤需要消耗 w 点能量。
【 输出格式】
输出文件名为 calabash.out。
一行一个实数， 表示答案（ 误差不超过 10^-3）。
【 输入输出样例】
calabash.in calabash.out
4 6
1 2 1
2 4 6
1 3 2
3 4 4
2 3 3
1 4 8
2.000【输入输出样例说明】
有 4 种爬法：
1->4，消耗能量 8，经过 2 个葫芦， 比值为 8/2=4。
1->2->4， 消耗能量 1+6=7， 经过 3 个葫芦， 比值为 7/3≈2.33。
1->3->4， 消耗能量 2+4=6， 经过 3 个葫芦， 比值为 6/3=2。
1->2->3->4， 消耗能量 1+3+4=8， 经过 4 个葫芦， 比值为 8/4=2。
所以选第三种或第四种方案， 答案为 2。
【数据规模与约定】
测试点编号 n m 特殊说明
1 2 1
2 100 99 除 1 外，所有葫芦的入度均为 1
3 100 105 所有从 1到 n的路径经过的葫芦数相等
4 100 1000
5 100 1000
6 199 198 除 1 外，所有葫芦的入度均为 1
7 200 231 所有从 1到 n的路径经过的葫芦数相等
8 200 2000
9 200 2000
10 200 2000
对于所有数据，小澳爬过每条藤消耗的能量不会超过 10^3， 且一定存在一条从 1
到 n 的路径。

根据题解来说，这是一道最小比例生成树的题。
另建一个起点 0，连接一条 0 到 1 长度为 0 的边，就此将问题转化为长度和
边数最小比值。 这个问题的求解需要分数规划。
假设答案为 ans，对于任意一条由 k 条边组成的路径，有：
(w1+w2+w3+…+wk)/k>=ans；
转化一下：
(w1+w2+w3+…+wk) >=ans*k；
即(w1-ans)+(w2-ans)+(w3-ans)+…+(wk-ans)>=0。
于是就得到了这样一个算法：
二分答案 x，每次将每一条边的权值减去 x 求最短路，判断 1~n 的最短路是
否大于 0： 若大于 0，则说明答案 ans>x；否则说明 ans

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define inf 0x7fffffff
#define lim 1e-4 
#define M 2010
#define N 210
using namespace std;

queue<int>q;double dis[N];int n,m;
bool et[N];
struct e{
    int next,to;double w,x;
}a[M];int cnt,last[N];

void bud(int u,int v,double w)
{
    a[++cnt].to=v;a[cnt].x=w;a[cnt].next=last[u];
    last[u]=cnt;
}

bool spfa(double x)
{
    for(int i=1;i<=cnt;++i)a[i].w=a[i].x-x;
    for(int i=1;i<=n;++i)dis[i]=inf;
    dis[0]=0;et[0]=1;q.push(0);
    while(!q.empty()){
        int x=q.front();q.pop();et[x]=0;
        for(int i=last[x];i;i=a[i].next){
            if(dis[a[i].to]-lim>dis[x]+a[i].w)
            dis[a[i].to]=dis[x]+a[i].w;
            if(!et[a[i].to])q.push(a[i].to),et[a[i].to]=1;      
        }

    }
    return dis[n]<-lim;
}
int main()
{
    freopen("calabash.in","r",stdin);
    freopen("calabash.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;++i)
    {
        int u,v;double w;
        scanf("%d%d%lf",&u,&v,&w);
        bud(u,v,w);
    }
    bud(0,1,0);
    double l=0,r=1e3;
    while(r-l>=lim){
        double mid=(l+r)/2.0;
        if(spfa(mid))r=mid;
        else     l=mid;
    }
    printf("%.3lf",l);
    return 0;
}