[luogu1768]天路

##链接##
题目描述

“那是一条神奇的天路诶_{，把第一个神犇送上天堂}”，XDM先生唱着这首“亲切”的歌曲，一道猥琐题目的灵感在脑中出现了。

和C_SUNSHINE大神商量后，这道猥琐的题目终于出现在本次试题上了，旨在难到一帮大脑不够灵活的OIer们（JOHNKRAM真的不是说你……）。

言归正传，小X的梦中，他在西藏开了一家大型旅游公司，现在，他要为西藏的各个景点设计一组铁路线。但是，小X发现，来旅游的游客都很挑剔，他们乘火车在各个景点间游览，景点的趣味当然是不用说啦，关键是路上。试想，若是乘火车一圈转悠，却发现回到了游玩过的某个景点，花了一大堆钱却在路上看不到好的风景，那是有多么的恼火啊。

所以，小X为所有的路径定义了两个值，Vi和Pi，分别表示火车线路的风景趣味度和乘坐一次的价格。现在小X想知道，乘客从任意一个景点开始坐火车走过的一条回路上所有的V之和与P之和的比值的最大值。以便为顾客们推荐一条环绕旅游路线（路线不一定包含所有的景点，但是不可以存在重复的火车路线）。

于是，小X梦醒之后找到了你……

输入输出格式

输入格式：
第一行两个正整数N，M，表示有N个景点，M条火车路线，火车路线是单向的。

以下M行，每行4个正整数，分别表示一条路线的起点，终点，V值和P值。

注意，两个顶点间可能有多条轨道，但一次只能走其中的一条。

输出格式：
一个实数，表示一条回路上最大的比值，保留1位小数。

若没有回路，输出-1。

输入输出样例

输入样例#1：
5 6
1 2 1 1
4 1 6 2
5 4 8 1
2 3 2 2
5 2 4 1
3 5 6 4
输出样例#1：
2.3
说明

对于30%的数据，1≤N≤100,1≤M≤20;

对于60%的数据，1≤N≤3,000,1≤M≤2,000;

对于100%的数据，1≤N≤7,000,1≤M≤20,000,1≤Vi,Pi≤1,000.

保证答案在200以内.

二分+spfa判负环
*我们这样想：设ans为答案，则存在任意环E保证∑V(i)/∑P(i)<=ans，可以得到∑V(i)<=ans*∑P(i)，于是∑{P(i)ans}-∑V(i)>=0。

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

const int maxn = 100000 + 100;
int n,m;
struct edge {
    int u,v;
	dd w1,w2;
    int next;
}e[maxn * 2];
int head[maxn],tot = 0;

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;
}

void add(int u, int v, int w1,int w2) {
    e[++tot] = (edge){u,v,w1,w2,head[u]};
    head[u] = tot;
}

int vis[maxn],p[maxn];
dd d[maxn];

bool spfa(dd mid, int x) {
	vis[x] = 1;
	for(int i = head[x]; i; i = e[i].next) {
		int v = e[i].v;
		if(d[v] > d[x] - e[i].w1 + e[i].w2 * mid) {
			d[v] = d[x] - e[i].w1 + e[i].w2 * mid;
			if(!vis[v]) {
				if(spfa(mid,v)) {
					return 1;
				}
			}
			else {
				return 1;
			}
		}
	}
	vis[x] = 0;
	return 0;
}

int main() {
    n = read(), m = read();
    for(int i = 1; i <= m; i++) {
    	int u = read(), v = read(), w1 = read(), w2 = read();
    	add(u,v,w1,w2);
    }
    for(int i = 1; i <= n; i++) add(0,i,0,0);
    dd l = 0, r = 200,ans = 0;
    while( l + 0.01 < r ) {
    	dd mid = (l + r) / 2;
    	memset(d,127,sizeof(d));
    	memset(vis,0,sizeof(vis));
    	d[0] = 0;
    	if(spfa(mid,0)) {
    		ans = mid;
    		l = mid;
    	}
    	else r = mid;
    }
    if(ans == 0) cout<<"-1"<<endl;
    else printf("%.1lf",ans);
	return 0;
}