【HNOI 2009】最小圈

传送门

---

Problem

给定一个有 $n$ 个点 $m$ 条边的有向图，边带权。在图中找到一个环，这个环上的边权平均值最小，求出这个最小值。

注意边权可以为负。

数据范围：$n\le 3000,m\le 10000,\mid w_{i,j}\mid \le 10^7$。

---

Solution

很显然这是一道分数规划的题。

假设环上有 $k$ 个点（也即有 $k$ 条边），且边权分别 $w_1,w_2,...,w_k$。那我们就是要让 $\frac{{\sum }_{i=1}^k{w}_i}{k}$ 最小。

设这个最小值为 $\lambda$，那对于图中任意一个环，都有 $\frac{{\sum }_{i=1}^{k^{\prime }}{w}_i}{k^{\prime }}\ge \lambda$。

那稍微化简一下就是 ${\sum }_{i=1}^{k^{\prime }}{w}_i-k^{\prime }\lambda \ge 0$，也即 ${\sum }_{i=1}^{k^{\prime }}(w_i-\lambda )\ge 0$。

于是我们二分答案，令 $mid$ 为当前的二分值，将所有的边权变为 $w_i-mid$。可以发现，当图中出现负环时，$mid>\lambda$，应当将 $mid$ 调小；反之应将 $mid$ 调大。

**注意：**这道题卡 $bfs$ 版的写法，要写成 $dfs$ 版才能过。

---

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 50005
#define eps 1e-10
using namespace std;
int n,m,t;
int first[N],v[N],nxt[N],vis[N];
double w[N],d[N];
void add(int x,int y,double z){
	nxt[++t]=first[x],first[x]=t,v[t]=y,w[t]=z;
}
bool spfa(int x,double mid){
	vis[x]=1;
	for(int i=first[x];i;i=nxt[i]){
		int to=v[i];
		double len=w[i]-mid;
		if(d[to]>d[x]+len){
			d[to]=d[x]+len;
			if(vis[to]||spfa(to,mid))  return true;
		}
	}
	vis[x]=0;
	return false;
}
bool check(double mid){
	int i,flag=0;
	for(i=1;i<=n;++i)
		d[i]=0,vis[i]=0;
	for(i=1;i<=n;++i)
		if(spfa(i,mid))
			flag=1;
	return flag;
}
int main(){
	int x,y,i;double z;
	scanf("%d%d",&n,&m);
	for(i=1;i<=m;++i){
		scanf("%d%d%lf",&x,&y,&z),add(x,y,z);
	}
	double l=-1e6,r=1e6;
	while(r-l>eps){
		double mid=(l+r)/2;
		if(check(mid))  r=mid;
		else  l=mid;
	}
	printf("%.8lf",l);
	return 0;
}