[省选前题目整理][BZOJ 1486][HNOI 2009]最小圈(01分数规划)

题目链接

http://www.lydsy.com/JudgeOnline/problem.php?id=1486

思路

首先我们知道，将一个圈中所有的边权全部减去它们之和的平均数的话，处理后的圈的边权和为0。因此我们可以二分答案，每次二分出答案后，将图上所有的边权全部减去这个答案，如果这个答案可行的话(但是这个答案不一定是最大的)，图上一定会存在负环(环的边权之和是小于等于0)，因此我们用SPFA判负环，如果图上存在负环，则说明当前二分的答案大了，降低上界。否则说明当前二分的最大答案小了，增大下界。

代码

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
#include <cmath>

#define MAXE 11000
#define MAXV 4100
#define INF 0x3f3f3f3f
#define EPS 1e-10

using namespace std;

int n,m;

struct edge
{
    int u,v,next;
    double w,w0; //w0是原来的边权
}edges[MAXE];

int head[MAXV],nCount=0;

void AddEdge(int U,int V,double W)
{
    edges[++nCount].u=U;
    edges[nCount].v=V;
    edges[nCount].w0=W;
    edges[nCount].next=head[U];
    head[U]=nCount;
}

bool hasFind,vis[MAXV]; //hasFind=true表示SPFA已经找到了负环
double dist[MAXV];

void SPFA(int u)
{
    vis[u]=true;
    for(int p=head[u];p!=-1;p=edges[p].next)
    {
        int v=edges[p].v;
        if(dist[u]+edges[p].w<dist[v])
        {
            if(vis[v])
            {
                hasFind=true;
                return;
            }
            dist[v]=dist[u]+edges[p].w;
            SPFA(v);
            if(hasFind) return;
        }
    }
    vis[u]=false; //!!!!
}

bool check() //判断当前的图是否有负环
{
    hasFind=false;
    memset(vis,false,sizeof(vis));
    memset(dist,0,sizeof(dist));
    for(int i=1;i<=n;i++)
    {
        SPFA(i);
        if(hasFind) return true;
    }
    return false;
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int u,v;
        double w;
        scanf("%d%d%lf",&u,&v,&w);
        AddEdge(u,v,w);
    }
    double ans;
    double lowerBound=-1e9,upperBound=1e9;
    while(fabs(upperBound-lowerBound)>EPS)
    {
        double mid=(upperBound+lowerBound)/2;
        for(int i=1;i<=nCount;i++) edges[i].w=edges[i].w0-mid;
        bool tmp=check();
        if(tmp) upperBound=mid;
        else lowerBound=mid,ans=mid; //!!!!
    }
    printf("%.8lf\n",ans);
    return 0;
}