BZOJ 3597 [SCOI 2014] 方伯伯运椰子 （分数规划+消圈算法）

题目链接：BZOJ 3597

因为总流量不可改变，只能改变流量的分布，所以就用消圈算法吧。。。

因为要求(x-y)/k的最大值，所以用0-1分数规划吧。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;

const int maxn=5000+10;

double eps=1e-3;

int N,M,col;
bool vis[maxn];
double dis[maxn];
int head[maxn],fa[maxn];

struct node{int v,w,next;}e[100000];

inline int read(){
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

int k=0;

void adde(int u,int v,int w){
	e[k].v=v; e[k].w=w; e[k].next=head[u]; head[u]=k++;
}

void input(){
	memset(head,-1,sizeof(head));
	
	N=read(); M=read();
	for(int i=1;i<=M;i++){
		int u=read(),v=read(),a=read(),b=read(),c=read(),d=read();
		adde(u,v,b+d);  if(c)adde(v,u,a-d);
	}
}

void spfa(int x,double mid){
	if(col)return ;
	
	vis[x]=1;
	
	for(int i=head[x];i!=-1;i=e[i].next){
		int v=e[i].v;
		if(dis[x]+e[i].w+mid<dis[v]){
			dis[v]=dis[x]+e[i].w+mid;
			
			if(vis[v]){
				col=1; return ;
			}
			spfa(v,mid);
		}
	}
	vis[x]=0;
}

bool check(double mid){
	memset(dis,0,sizeof(dis));
	memset(vis,0,sizeof(vis));
	
	col=0;
	
	for(int i=1;i<=N+2;i++){
		spfa(i,mid);
		if(col)return true;
	}
	return false;
}

void solve(){
	double L=0,R=1e10,mid;
	while(R-L>eps){
		mid=(L+R)/2;
		if(check(mid))L=mid;
		else R=mid;
	}
	printf("%.2lf\n",(L+R)/2);
}

int main(){
	input();
	solve();
	return 0;
}