BZOJ 1486: [HNOI2009]最小圈（01分数规划：最优比例环）

1486: [HNOI2009]最小圈

Time Limit: 10 Sec   Memory Limit: 64 MB
Submit: 2675   Solved: 1279
[ Submit][ Status][ Discuss]

Description

Input

Output

Sample Input

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

Sample Output

3.66666667

HINT

Source

题解：之前已经学习了01分数规划问题，这里不再回顾，对于所有的01分数规划问题，我们无非就是找到

最优的x满足题目要求的最大最小值，对于x，可证明出其单调性，因此我们大多考虑二分查找，对于本题，因为

要求的是最小值，对于题上的公式，我们转化为以下形式：

                                            

因为证明出ans的单调性，我们考虑二分枚举ans，然后通过dfs查找图中的环，对于当前所有边的权值

我们将其更新为val-ans，因此只要我们在图中找到负环，说明会有更小的ans符合答案。。。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define eps 1e-10
#define inf 1000000000
struct node
{
	int u,v;
	double c;
}a[30005],e[100005];
int n,m,tot,flag;
int p[100005],nxt[100005];
double dis[3005];
bool used[3005];
void add(int x,int y,double z)
{
	tot++;nxt[tot]=p[x];p[x]=tot;
	a[tot].u=x;a[tot].v=y;a[tot].c=z;
}
void dfs(int x)
{
	used[x]=0;
	for(int i=p[x];i;i=nxt[i])
		if(dis[a[i].v]>dis[x]+a[i].c)   
			if(!used[a[i].v])   
			{   
				flag=1;   
				break;      
			}
			else
			{
				dis[a[i].v]=dis[x]+a[i].c;
				dfs(a[i].v);
			}
	used[x]=1;
}
bool check(double x)
{
	tot=0;
	memset(p,0,sizeof(p));
	memset(nxt,0,sizeof(nxt));
	memset(dis,0,sizeof(dis));
	memset(used,1,sizeof(used));
	for(int i=1;i<=m;i++)
		add(e[i].u,e[i].v,e[i].c-x);
	flag=0;
	for(int i=1;i<=n;i++)
	{
		dfs(i);
		if(flag) return 1;
	}
	return 0;
}
int main(void)
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
		scanf("%d%d%lf",&e[i].u,&e[i].v,&e[i].c);
	double l=-inf,r=inf;
	while(l+eps<r)
	{
		double mid=(l+r)/2.0;
		if(check(mid))
			r=mid;
		else
			l=mid;
	}
	printf("%.8f\n",l);
	return 0;
}