bzoj1486: [HNOI2009]最小圈

传送门
显然是要二分答案的。
然后爆枚起点，用dfs版的spfa判负环就ok辣。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#define N 3005
#define ll long long
using namespace std;
struct edge{
    int to,next;
    double v;
}e[10005];
int head[N],v[N];
int n,m,x,y,flag,tot;
double l,r,z,mid,dis[N];
void spfa(int x,double mid){
    if (flag) return;
    v[x]=1;
    for (int i=head[x];i;i=e[i].next)
        if (dis[x]+e[i].v-mid<dis[e[i].to]){
            if (v[e[i].to]){
                flag=1;
                return;
            }
            dis[e[i].to]=dis[x]+e[i].v-mid;
            spfa(e[i].to,mid);
        }
    v[x]=0;
}
inline bool jud(double mid){
    memset(dis,0,sizeof(dis));
    memset(v,0,sizeof(v));
    flag=0;
    for (int i=1;i<=n;i++){
        spfa(i,mid);
        if (flag) return 1;
    }
    return 0;
}
int main(){
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++){
        scanf("%d%d%lf",&x,&y,&z);
        e[++tot].to=y;
        e[tot].v=z;
        e[tot].next=head[x];
        head[x]=tot;
    }
    l=-10000000;
    r=10000000;
    while (l+1e-9<r){
        mid=(l+r)/2;
        if (jud(mid)) r=mid; else l=mid;
    }
    printf("%.8lf",l);
}