bzoj2395 [Balkan 2011]Timeismoney(最小乘积生成树+计算几何)
题意
每条边有两个权值\(c,t\),请求出一颗生成树,使得\(\sum c\times \sum t\)最小
题解
为什么生成树会和计算几何扯上关系……
对于每棵树,设\(x=c,y=t\),我们可以把它看成平面上的一个点,其中\(\sum c\)为横坐标,\(\sum t\)为纵坐标。那么题目就可以转化成求反比例函数图像上的点\(k=xy\)满足\(k\)最小
我们先求出一棵\(\sum c\)最小的生成树,设这个点为\(A\),和一棵\(\sum t\)最小的生成树,设为\(B\),那么如果一个点\(C\)它们更优它肯定在直线\(AB\)的下方
证明:因为\(C\)比\(A,B\)更优,所以\(A,B\)肯定在这个反比例函数图像的上方。如果点\(A\)不在直线下方,画个图就会发现它必定有纵坐标小于\(B\)或者横坐标小于\(A\),与\(A,B\)的定义矛盾。故如果该点存在必定在直线\(AB\)下方
那么我们需要寻找一个在直线\(AB\)下方且离这条直线最远的一个点\(C\)来更新答案,也就是说要满足\(AB\)叉乘\(AC\)最小,即最小化
Ans
&=(B.x-A.x)(C.y-A.y)-(C.x-A.x)(B.y-A.y)\\
&=C.y(B.x-A.x)-C.x(B.y-A.y)-A.y(B.x-A.x)+A.x(B.y-A.y)\\
\end{aligned}
\]
后面两项是常数不用考虑,那么就是要最小化\(C.y(B.x-A.x)-C.x(B.y-A.y)\)。我们把每条边的权值设为\(y(B.x-A.x)-x(B.y-A.y)\),跑一遍最小生成树就能得到\(C\)了。之后再对\((A,C),(C,B)\)递归下去即可
//minamoto
#include<bits/stdc++.h>
#define R register
#define ll long long
#define inf 1e18
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
R int res,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
const int N=205,M=10005;
struct eg{
int u,v,c,t,w;
inline bool operator <(const eg &b)const{return w<b.w;}
inline void rd(){u=read()+1,v=read()+1,c=read(),t=read();}
}e[M];
struct node{
int x,y;
node(){}
node(R int xx,R int yy):x(xx),y(yy){}
inline node operator -(const node &b)const{return node(x-b.x,y-b.y);}
inline ll operator *(const node &b)const{return 1ll*x*b.y-1ll*y*b.x;}
inline bool operator >(const node &b)const{return 1ll*x*y==1ll*b.x*b.y?x>b.x:1ll*x*y>1ll*b.x*b.y;}
}mc,mt,res(inf,inf);
int n,m,u,v,fa[N];
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
node kruskal(){
int tot=0;node P(0,0);
fp(i,1,n)fa[i]=i;
fp(i,1,m){
u=find(e[i].u),v=find(e[i].v);
if(u!=v){
fa[u]=v,++tot;
P.x+=e[i].c,P.y+=e[i].t;
if(tot==n-1)break;
}
}
return cmin(res,P),P;
}
void solve(node A,node B){
fp(i,1,m)e[i].w=1ll*e[i].t*(B.x-A.x)-1ll*e[i].c*(B.y-A.y);
sort(e+1,e+1+m);
node C=kruskal();
if((B-A)*(C-A)>=0)return;
solve(A,C),solve(C,B);
}
int main(){
// freopen("testdata.in","r",stdin);
n=read(),m=read();
fp(i,1,m)e[i].rd();
fp(i,1,m)e[i].w=e[i].c;
sort(e+1,e+1+m),mc=kruskal();
fp(i,1,m)e[i].w=e[i].t;
sort(e+1,e+1+m),mt=kruskal();
solve(mc,mt);
printf("%d %d\n",res.x,res.y);
return 0;
}
bzoj2395 [Balkan 2011]Timeismoney(最小乘积生成树+计算几何)的更多相关文章
- bzoj2395[Balkan 2011]Timeismoney最小乘积生成树
所谓最小乘积生成树,即对于一个无向连通图的每一条边均有两个权值xi,yi,在图中找一颗生成树,使得Σxi*Σyi取最小值. 直接处理问题较为棘手,但每条边的权值可以描述为一个二元组(xi,yi),这也 ...
- Bzoj2395: [Balkan 2011]Timeismoney(最小乘积生成树)
问题描述 每条边两个权值 \(x,y\),求一棵 \((\sum x) \times (\sum y)\) 最小的生成树 Sol 把每一棵生成树的权值 \(\sum x\) 和 \(\sum y\) ...
- bzoj 2395 [Balkan 2011]Timeismoney——最小乘积生成树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2395 如果把 \( \sum t \) 作为 x 坐标,\( \sum c \) 作为 y ...
- 【BZOJ2395】【Balkan 2011】Timeismoney 最小乘积生成树
链接: #include <stdio.h> int main() { puts("转载请注明出处[辗转山河弋流歌 by 空灰冰魂]谢谢"); puts("网 ...
- bzoj2395: [Balkan 2011]Timeismoney
Description 有n个城市(编号从0..n-1),m条公路(双向的),从中选择n-1条边,使得任意的两个城市能够连通,一条边需要的c的费用和t的时间,定义一个方案的权值v=n-1条边 ...
- BZOJ2395 [Balkan 2011]Timeismoney 【最小乘积生成树】
题目链接 BZOJ2395 题意:无向图中每条边有两种权值,定义一个生成树的权值为两种权值各自的和的积 求权值最小的生成树 题解 如果我们将一个生成树的权值看做坐标,那么每一个生成树就对应一个二维平面 ...
- 【最小乘积生成树】bzoj2395[Balkan 2011]Timeismoney
设每个点有x,y两个权值,求一棵生成树,使得sigma(x[i])*sigma(y[i])最小. 设每棵生成树为坐标系上的一个点,sigma(x[i])为横坐标,sigma(y[i])为纵坐标.则问题 ...
- P5540-[BalkanOI2011]timeismoney|最小乘积生成树【最小生成树,凸壳】
正题 题目链接:https://www.luogu.com.cn/problem/P5540 题目大意 给出\(n\)个点\(m\)条边边权是一个二元组\((a_i,b_i)\),求出一棵生成树最小化 ...
- 洛谷 P5540 - [BalkanOI2011] timeismoney | 最小乘积生成树(最小生成树)
洛谷题面传送门 大概是一个比较 trivial 的小 trick?学过了就不要忘了哦( 莫名奇妙地想到了 yyq 的"hot tea 不常有,做过了就不能再错过了" 首先看到这种二 ...
随机推荐
- 优化深度神经网络(一) dropout 初始化
Coursera吴恩达<优化深度神经网络>课程笔记(1)-- 深度学习的实用层面 1. Train/Dev/Test sets 训练集(Training sets).验证集(Develo ...
- Java面向对象-代码块
Java面向对象-代码块 代码块主要就是通过{}花括号 括起来的代码: 主要分为 普通代码块 构造块 静态代码块三类.后面学到线程还有一个同步代码块,到时候再说: 普通代码块:仅仅是花括号括起来的代码 ...
- leetcode897
这道题用C++来写,在本地执行正常,但是使用OJ判断输出结果是空,暂时不清楚原因.代码如下: class Solution { public: vector<int> V; //中序遍历 ...
- 【知识碎片】JavaScript篇
40.选择器组合 逗号是多选择器空格 是子子孙孙尖括号 只找儿子 39.失去焦点事件blur $("input").blur(function(){ $("input& ...
- 【286】◀▶ Python 内置函数说明
参考: Python 内置函数 01 abs() 返回数字的绝对值. 02 all() 用于判断给定的可迭代参数 iterable 中的所有元素是否不为 0.''.False 或者 itera ...
- SonarQube在CentOS上的安装
1 简介 SonarQube 是一个用于代码质量管理的开放平台.通过插件机制,Sonar 可以集成不同的测试工具,代码分析工具,以及持续集成工具.与持续集成工具(例如 Hudson/Jenkins 等 ...
- S3C6410移植u-boot
1.首先下载u-boot(ftp://ftp.denx.de/pub/u-boot) wget ftp://ftp.denx.de/pub/u-boot/u-boot-latest.tar.bz2 2 ...
- c语言实践 统计输入的一串正整数里面奇数和偶数的个数
怎么考虑这个问题. 首先先确定肯定是需要一个变量保存输入的数据的,我们叫它input,最后结果要的是个数,所以需要另外两个变量来保存奇数的个数和偶数的个数. int input int countJ ...
- ZIP压缩格式与RAR压缩格式
早已习惯了安装系统之后必须安装winrar,压缩文件也已经习惯了rar格式,这种习惯的力量真的挺可怕的.在工作中你的同事可能没有安装winrar,或者他们不喜欢安装盗版软件,这时候你给他们发送过去的是 ...
- CentOS 安装mongodb3.0 二进制包
1.下载mongodb因为64位系统CentOS,所以下载64位的安装包: wget https://fastdl.mongodb.org/linux/mongodb-linux-x86_64-3.0 ...