BZOJ 1486: [HNOI2009]最小圈|分数规划

分数规划

spfa判负环

dis数组直接清0（好神奇啊而且速度好快啊）

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define T 22222
using namespace std;
double v[T],e[T],dis[T];
int head[T],nxt[T],lst[T],mark[T];
int u[T],w[T];
int n,m,tot,flag;
void insert(int x,int y,double z)
{
	lst[++tot]=y;
	nxt[tot]=head[x];
	v[tot]=z;
	head[x]=tot;
}
void build(double x)
{
	for(int i=1;i<=m;i++)
	    insert(u[i],w[i],e[i]-x);
}
void spfa(int x)
{
	mark[x]=1;
	for(int i=head[x];i;i=nxt[i])
	    if(dis[x]+v[i]<dis[lst[i]])
	    {
	    	dis[lst[i]]=dis[x]+v[i];
	    	if(mark[lst[i]]){flag=1;break;}
	    	if(flag)break;
	    	spfa(lst[i]);
	    }
	mark[x]=0;
}
bool jud(double x)
{
	for(int i=1;i<=n;i++)head[i]=0;
	for(int i=1;i<=n;i++)dis[i]=0;
	tot=0; 	build(x);
	flag=0;	
	for(int i=1;i<=n;i++)
	{
		spfa(i);
	    if(flag)return 1;
	}
	return 0;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)scanf("%d%d%lf",&u[i],&w[i],&e[i]);
	double l=-10000000,r=10000000;
	while(r-l>1e-9)
	{
		double mid=(l+r)/2;bool x=jud(mid);
		if(x)r=mid;
		else l=mid;
	}
	printf("%.8lf",l);
}