Usaco 地震:
Description
一场地震把约翰家的牧场摧毁了，坚强的约翰决心重建家园。约翰已经重建了N个牧场， 现在他希望能修建一些道路把它们连接起来。研究地形之后，约翰发现可供修建的道路有M 条。碰巧的是，奶牛们最近也成立一个工程队，专门从事修复道路。而然，奶牛们很有经济 头脑，如果无利可图，它们是不会干的。
奶牛们关注的是挣钱速度，即总利润和总施工时间的比值。约翰和奶牛达成了协议，奶 牛负责修建道路，将所有牧场连通，而约翰需要支付F元。每条道路都有自己的施工时间和 建造成本。连接两个相同的牧场的道路可能有多条。保证所有的牧场必定是可连通的，不过 也有可能一些道路的建造成本之和会超过F。
请帮助奶牛们选择修复哪些道路，才能使单位时间的利润最大？
Input Format
第一行：三个整数： N，M和F，1≤N≤400，1≤M≤10000，1≤F≤2×10^9
第二行到M+1行：第i +1行表示第i条道路的信息，有四个整数：Ui，Vi，Ci和Ti，Ui和Vi表示这条道路连接的牧场编号，Ci表示这条路的施工时间，Ti表示建造成本，1≤Ui≤N，1≤Vi≤N，1≤Ci≤2×10^9，1≤Ti≤2×10^9
Output Format
第一行：一个保留四位小数的浮点数，表示奶牛们能挣到的最大单位时间利润，如果奶 牛们无钱可赚，则输出0.0000
-------------------------------------------------------------------------------------
显然这题是求最优比率生成树，
对最优比率 k 而言 (∑t)*k=F-∑v
F=(∑t)*k+∑v=∑(t*k+v)
显然上面这个式子具有单调性，随 k 而单调递增
因此可以对 k 进行二分，
以 ( k*t +v ) 作为关键字做生成树
若∑(t*k+v)>F 下调 k 否则上调 k
Ps.注意精度问题，图图大神说 ：二分的精度至少要比保留的小数位数多2位。例如这题要求保留到0.0001，那么二分的精度应该为0.000001
代码如下：
1 #include<cstring> 2 #include<algorithm> 3 #include<cstdio> 4 #include<string> 5 #include<iostream> 6 #define O 0.00000001 7 #define LL long long 8 #define INF 2000000000 9 using namespace std; 10 struct Point{ 11 int x,y; 12 double c,t; 13 }a[10011]; 14 int f[500],n,m; 15 double l,r,ans,F,mid; 16 bool cmp(const Point&X,const Point&Y){ 17 return X.c+mid*X.t<Y.c+mid*Y.t; 18 } 19 int find(int u){ 20 return f[u]==u ? f[u] : f[u]=find(f[u]); 21 } 22 bool check(long double lim){ 23 sort(a+1,a+1+m,cmp); 24 for(int i=1;i<=n;i++) f[i]=i; 25 int tot=0; 26 double temp=0; 27 for(int i=1;i<=m;i++) 28 if(find(a[i].x)!=find(a[i].y)){ 29 f[find(a[i].x)]=f[find(a[i].y)]; 30 temp+=a[i].c+mid*a[i].t; 31 ++tot; 32 if(temp>F) return false ; 33 } 34 if(tot<n-1) { 35 printf("0.000"); 36 exit(0); 37 } 38 return true; 39 } 40 int main(){ 41 scanf("%d%d%lf",&n,&m,&F); 42 for(int i=1;i<=m;i++) 43 scanf("%d%d%lf%lf",&a[i].x,&a[i].y,&a[i].c,&a[i].t); 44 l=0,r=F; 45 while(r-l>=O){ 46 mid=(l+r)/2; 47 if(check(mid)){ 48 ans=mid; 49 l=mid+O; 50 } else r=mid-O; 51 } 52 printf("%.4lf",ans); 53 }
-------------------------------------------------------------------------------------
Usaco 奶牛观光：
Description
作为对奶牛辛勤工作的回报，约翰决定带她们去附近的大城市玩一天。这个城市有L个景点，参观第i个景点会给奶牛带来Fi点欢乐度。第二天一早，奶牛可以自由选择从一个景点出发，约翰会负责开车把她们送到那里，但她们晚上必须回到这个景点和约翰汇合。
大城市里都是单行道，第i条道路从第L1i个建筑通向第L2i个建筑，走完需要Ti的时间。奶牛讨厌走路，定义一条游览线路的“欢乐指数”为参观这条线路上所有景点的欢乐度之和与花在路上的时间之和的比值。欢乐指数越大的线路越受欢迎。当然，参观同一景点两次不会带来双倍的欢乐。假设奶牛们至少参观两个景点，请帮她们找到一条欢乐指数最大的线路。
Input Format
第一行：两个用空格分开的整数：L和P，1 ≤ L ≤ 1000，2 ≤ P ≤ 5000
第二行到L + 1行：第i + 1行有一个整数Fi，1 ≤ Fi ≤ 1000
第L + 2行到L + P + 1行：第i + L + 1行有三个用空格分开的整数：L1i，L2j和Ti，1 ≤ L1i ≤ L，1 ≤ L2i ≤ L，1 ≤ Ti ≤ 1000
Output Format
第一行：输出一个实数，表示最大的欢乐指数，保留两位小数
-----------------------------------------------------------------
好把，这题是求最优比率最短路(环)( Usaco真是一站式服务- =)
最优比率 (∑t)*k=∑f => ∑(t*k-f)=0
和上一题一样，k 具有单调性，也显然可以二分
那如何判断这个环是否小于 0(负权环)
利用SPFA：
图图大神的方法是：
用edge[i]表示到i的最短路经过的边数，若边数 = N,则存在负权环
我用了种比较傻的：
time[i]表示该点被更新了几次，如果大于 n 显然存在负权环
图图大神还说：SPFA的使用十分灵活，经常配合动规使用，所以往往队列里的初始状态是不唯一的。
代码如下:
1 #include<cstring> 2 #include<algorithm> 3 #include<cstdio> 4 #include<string> 5 #include<iostream> 6 #include<queue> 7 #define INF 999999999 8 #define N 5011 9 const double E=1e-4 ; 10 using namespace std; 11 queue <int> q; 12 int n,P; 13 int T,last[N],next[N],s[N]; 14 double sum[N],t[N],l,r,mid,ans; 15 bool check(){ 16 int time[N]; memset(time,0,sizeof(time)); 17 double d[N]; for(int i=1;i<=n;i++) d[i]=INF; 18 bool in[N]; memset(in,0,sizeof(in)); 19 time[1]=in[1]=1; d[1]=0; 20 while(!q.empty()) q.pop(); 21 for(q.push(1);!q.empty();q.pop()){ 22 int now=q.front(); 23 for(int i=last[now];i;i=next[i]){ 24 double w=t[i]*mid-sum[now]; 25 if(d[s[i]]+E>d[now]+w){ 26 d[s[i]]=d[now]+w; 27 if(!in[s[i]]){ 28 in[s[i]]=1; 29 q.push(s[i]); 30 if(++time[s[i]]>n) return 1; 31 } 32 } 33 } 34 in[now]=0; 35 } 36 return 0; 37 } 38 int main(){ 39 scanf("%d%d",&n,&P); 40 for(int i=1;i<=n;i++) scanf("%lf",&sum[i]); 41 for(int i=1;i<=P;i++){ 42 int x,y; 43 double v; 44 scanf("%d%d%lf",&x,&y,&v); 45 next[++T]=last[x]; last[x]=T; s[T]=y; t[T]=v; 46 } 47 l=E,r=1000; 48 while(r-l>=E){ 49 mid=(l+r)/2; 50 if(check()){ 51 ans=mid; 52 l=mid+E; 53 } else r=mid-E; 54 } 55 printf("%.2lf",ans); 56 }