题目大意
有n个城市(编号从0..n-1),m条公路(双向的),
从中选择n-1条边,使得任意的两个城市能够连通,一条边需要的c的费用和t的时间,
定义一个方案的权值v=n-1条边的费用和*n-1条边的时间和,你的任务是求一个方案使得v最小
Input
第一行两个整数n,m,接下来每行四个整数a,b,c,t,表示有一条公路从城市a到城市b需要t时间和费用c
Output
【output】timeismoney.out
仅一行两个整数sumc,sumt,(sumc表示使得v最小时的费用和,sumc表示最小的时间和) 如果存在多个解使得sumc*sumt相等,输出sumc最小的
Sample Input
5 7
0 1 161 79
0 2 161 15
0 3 13 153
1 4 142 183
2 4 236 80
3 4 40 241
2 1 65 92
Sample Output
279 501
HINT
【数据规模】
1<=N<=200
1<=m<=10000
0<=a,b<=n-1
0<=t,c<=255
有5%的数据m=n-1
有40%的数据有t=c
对于100%的数据如上所述
我们把每个生成树的解(∑t,∑c)=>设为坐标(t,c),k=tc,构成反比例函数
要使k最小,c=k/t尽可能贴近t轴
所以从下凸包中找一个最小的
那么首先最极端的是单纯按t排序或按c排序得到的点A,B
先找到离(A,B)最远的C,在递归(A,C)和(C,B),边界为AC×BC>=0
这样就可以求出所有下凸包上的值
离AB最远的C必定导致SΔABC最大
即AB×AC最大,把常数去掉得到一个关于C的公式,把它带入最小生成树,算出答案
丢个有图的链接,不过边界好像有点问题
http://blog.csdn.net/regina8023/article/details/45246933
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long lol;
struct Node
{
int u,v;
lol t,c,w;
}E[];
struct ZYYS
{
lol x,y;
}A,B,ans;
int n,m,set[];
bool cmp(Node a,Node b)
{
return a.w<b.w;
}
int find(int x)
{
if (set[x]!=x) set[x]=find(set[x]);
return set[x];
}
ZYYS kruskal()
{lol i,cnt;
lol anst,ansc;
ZYYS C;
anst=;ansc=;
for (i=;i<=n;i++)
set[i]=i;
cnt=;
for (i=;i<=m;i++)
{
lol p=find(E[i].u),q=find(E[i].v);
if (p!=q)
{
set[p]=q;
cnt++;
anst+=E[i].t;
ansc+=E[i].c;
if (cnt==n-) break;
}
}
return (ZYYS){ansc,anst};
}
lol cross(ZYYS p,ZYYS a,ZYYS b)
{
return (a.x-p.x)*(b.y-p.y)-(b.x-p.x)*(a.y-p.y);
}
void solve(ZYYS a,ZYYS b)
{ZYYS C;
lol i;
lol y=a.y-b.y,x=b.x-a.x;
for (i=;i<=m;i++)
E[i].w=E[i].t*x+E[i].c*y;
sort(E+,E+m+,cmp);
C=kruskal();
if (C.x*C.y<ans.x*ans.y||(ans.x*ans.y==C.x*C.y&&C.x<ans.x))
ans=C;
if (cross(C,a,b)>=) return;
solve(a,C);
solve(C,b);
}
int main()
{
lol i;
cin>>n>>m;
ans.x=ans.y=1e9;
for (i=;i<=m;i++)
{
scanf("%d%d%lld%lld",&E[i].u,&E[i].v,&E[i].c,&E[i].t);
E[i].u++;
E[i].v++;
}
for (i=;i<=m;i++)
E[i].w=E[i].c;
sort(E+,E+m+,cmp);
A=kruskal();
if (A.x*A.y<ans.x*ans.y)
ans=A;
for (i=;i<=m;i++)
E[i].w=E[i].t;
sort(E+,E+m+,cmp);
B=kruskal();
if (B.x*B.y<ans.x*ans.y)
ans=B;
solve(A,B);
cout<<ans.x<<' '<<ans.y;
}