题目链接
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;
}