BZOJ1486 HNOI2009 最小圈 【01分数规划】

BZOJ1486 HNOI2009 最小圈

---

Description

考虑带权的有向图$G=(V,E)​$以及$w:E\rightarrow R​$ ,每条边$e=(i,j)(i\neq j,i\in V,j\in V)​$的权值定义为$w_{i,j}​$，令 $n=|V|​$ 。$c=(c_1,c_2,\cdots,c_k)(c_i\in V)​$是G中的一个圈当且仅当 $(c_i,c_{i+1})(1\le i<k)​$和$(c_k,c_1)​$都在E中，这时称k 为圈c的长度同时令$c_{k+1}=c_1​$，并定义圈 $c=(c_1,c_2,\cdots,c_k)​$的平均值为$μ(c)=\sum\limits_{i=1}^{k} w_{c_i,c_{i+1}}/k​$，即c上所有边的权值的平均值。令$\mu'(c)=Min(\mu(c))​$为G中所有圈c的平均值的最小值。现在的目标是：在给定了一个图$G=(V,E)​$以及$w:E\rightarrow R​$之后，请求出G中所有圈c的平均值的最小值 $\mu'(c)=Min(\mu(c)​)$

Input

第一行2个正整数，分别为n和m ，并用一个空格隔开，只用$n=|V|,m=|E|$分别表示图中有n个点m条边。 接下来m行，每行3个数$i,j,w_{i,j}$，表示有一条边$(i,j)$且该边的权值为$w_{i,j}$​。输入数据保证图$G=(V,E)$连通，存在圈且有一个点能到达其他所有点。

Output

请输出一个实数$\mu'(c)=Min(\mu(c))$ ，要求输出到小数点后8位。

Sample Input

4 5
1 2 5
2 3 5
3 1 5
2 4 3
4 1 3

Sample Output

3.66666667

---

应该算是01分数规划的裸板题了吧。。但是第一次写还是遇到了一些困难vis数组不清零之类的
假设一个答案成立，那么一定可以找到一个环使得其边权和大于等于边数$*ans$。
可以发现答案是具有二分性的，二分出一个临时答案$ans'$，并且用$ans'$对每条边进行约束，再用深搜SPFA判断一下负环，如果有负环说明当前解可行，继续二分就好了。
注意题目要求保留到小数点后八位，多开一点二分次数防止精度不够啊

---

/**************************************************************
    Problem: 1486
    User: yangkai
    Language: C++
    Result: Accepted
    Time:2292 ms
    Memory:1648 kb
****************************************************************/

#include<bits/stdc++.h>
using namespace std;
int read(){
    int ans=0,w=1;char c=getchar();
    while(c!='-'&&!isdigit(c))c=getchar();
    if(c=='-')w=-1,c=getchar();
    while(isdigit(c))ans=ans*10+c-'0',c=getchar();
    return ans*w;
}
#define N 3010
#define M 10010
#define INFF 10000000
int n,m,tot,head[N];
struct Edge{int v,next;double w;}E[M];
void add(int u,int v,double w){
    E[++tot]=(Edge){v,head[u],w};head[u]=tot;
}
int u[M],v[M],vis[N];
double w[M],dis[N];
void build(double val){
    for(int i=1;i<=n;i++)head[i]=0;tot=0;
    for(int i=1;i<=m;i++)add(u[i],v[i],w[i]-val);
} 
bool flag=0;
void SPFA(int x){
    vis[x]=1;
    for(int i=head[x];i;i=E[i].next){
        int y=E[i].v;
        if(dis[x]+E[i].w<dis[y]){
            if(vis[y]){flag=1;break;}
            else{
                dis[y]=dis[x]+E[i].w;
                SPFA(y);
            }
        }
    }
    vis[x]=0;
}
int main(){
    n=read();m=read();
    for(int i=1;i<=m;i++){
        u[i]=read(),v[i]=read();
        scanf("%lf",&w[i]);
    }
    double l=-INFF,r=INFF;
    for(int p=1;p<=60;p++){
        double mid=(l+r)/2;
        build(mid);
        flag=0;
        for(int i=1;i<=n;i++){
            memset(vis,0,sizeof(vis));
            memset(dis,0,sizeof(dis));
            SPFA(i);
            if(flag)break;
        }
        if(flag)r=mid;
        else l=mid;
    }
    printf("%.8lf",l);
    return 0;
}