Luogu5540 最小乘积生成树

题目链接:洛谷

题目描述:对于一个\(n\)个点\(m\)条边的无向连通图,每条边有两个边权\(a_i,b_i\),求使\((\sum a_i)\times (\sum b_i)\)最小的生成树。

数据范围:\(n\le 200,m\le 10000,a_i,b_i\le 255\)

这题是一道非常妙的计算几何题目。

我们对于每个生成树,用\((\sum a_i,\sum b_i)\)这个二维平面上的点来表示它,那么就是求所有点中横坐标乘纵坐标的最小值。

画画图就可以发现,答案只有可能在下凸包上,为什么呢?

因为如果\(C\)在线段\(AB\)上方,其中\(x_Ay_A=x_By_B\),因为反比例函数下凸,所以\(x_Cy_C>x_Ay_A\)。

但是生成树可能有很多个,怎么得到下凸包上的点呢?

Step1 求最靠近\(x,y\)轴的两个点\(A,B\)

为什么呢?因为\(A,B\)两个点必定在下凸包上面。令\(w_i=a_i\)或\(b_i\)用最小生成树求\(A,B\)。

Step2 求\(C\)在直线\(AB\)下方且\(S_{\Delta ABC}\)最大

为什么呢?因为\(C\)点必定在下凸包上面,否则它不是最大的点。我们发现

\[\begin{aligned}
-\frac{1}{2}S_{\Delta ABC}&=\overrightarrow{AB}\times \overrightarrow{AC} \\
&=(x_B-x_A)(y_C-y_A)-(y_B-y_A)(x_C-x_A) \\
&=((x_A-x_B)y_A-(y_A-y_B)x_A)+(x_B-x_A)y_C-(y_B-y_A)x_C
\end{aligned}
\]

令\(w_i=(x_B-x_A)b_i-(y_B-y_A)a_i\)就会得到\(C\)

Step3 将(A,C)和(C,B)代入Step2递归

这样就可以求出下凸包上所有的点。那什么时候终止递归呢?当然就是\(C\)不存在,或者说求出的\(C\)在\(AB\)上方。

时间复杂度为\(O(km\log m)\),其中\(k\)为下凸包上点的个数,在随机数据下不会很大。

#include<bits/stdc++.h>
#define Rint register int
using namespace std;
typedef long long LL;
const int N = 10003;
int n, m, a[N], b[N], fa[N];
struct Point {
int x, y;
inline Point(int _x = 0, int _y = 0): x(_x), y(_y){}
inline Point operator - (const Point &o) const {return (Point){x - o.x, y - o.y};}
} A, B, ans(1e9, 1e9);
inline LL cross(Point a, Point b){return (LL) a.x * b.y - (LL) b.x * a.y;}
struct Edge {
int u, v, w, id;
inline bool operator < (const Edge &o) const {return w < o.w;}
} e[N];
inline int getfa(int x){return x == fa[x] ? x : fa[x] = getfa(fa[x]);}
inline Point Kruskal(){
Point res; sort(e + 1, e + m + 1);
for(Rint i = 1;i <= n;i ++) fa[i] = i;
for(Rint i = 1, p = 1;i <= m && p < n;i ++){
int u = getfa(e[i].u), v = getfa(e[i].v);
if(u != v){fa[u] = v; ++ p; res.x += a[e[i].id]; res.y += b[e[i].id];}
}
LL Ans = (LL) ans.x * ans.y, Res = (LL) res.x * res.y;
if(Ans > Res || Ans == Res && ans.x > res.x) ans = res;
return res;
}
inline void solve(Point A, Point B){
for(Rint i = 1;i <= m;i ++) e[i].w = (A.y - B.y) * a[e[i].id] - (A.x - B.x) * b[e[i].id]; Point C = Kruskal();
if(cross(B - A, C - A) >= 0) return;
solve(A, C); solve(C, B);
}
int main(){
scanf("%d%d", &n, &m);
for(Rint i = 1;i <= m;i ++) scanf("%d%d%d%d", &e[i].u, &e[i].v, a + i, b + i), ++ e[i].u, ++ e[i].v, e[i].id = i;
for(Rint i = 1;i <= m;i ++) e[i].w = a[e[i].id]; Point A = Kruskal();
for(Rint i = 1;i <= m;i ++) e[i].w = b[e[i].id]; Point B = Kruskal();
solve(A, B); printf("%d %d", ans.x, ans.y);
}

Luogu5540 最小乘积生成树的更多相关文章

  1. bzoj2395[Balkan 2011]Timeismoney最小乘积生成树

    所谓最小乘积生成树,即对于一个无向连通图的每一条边均有两个权值xi,yi,在图中找一颗生成树,使得Σxi*Σyi取最小值. 直接处理问题较为棘手,但每条边的权值可以描述为一个二元组(xi,yi),这也 ...

  2. HDU5697 刷题计划 dp+最小乘积生成树

    分析:就是不断递归寻找靠近边界的最优解 学习博客(必须先看这个): 1:http://www.cnblogs.com/autsky-jadek/p/3959446.html 2:http://blog ...

  3. 【BZOJ2395】【Balkan 2011】Timeismoney 最小乘积生成树

    链接: #include <stdio.h> int main() { puts("转载请注明出处[辗转山河弋流歌 by 空灰冰魂]谢谢"); puts("网 ...

  4. 【算法】最小乘积生成树 & 最小乘积匹配 (HNOI2014画框)

    今天考试的时候果然题目太难于是我就放弃了……转而学习了一下最小乘积生成树. 最小乘积生成树定义: (摘自网上一篇博文). 我们主要解决的问题就是当k = 2时,如何获得最小的权值乘积.我们注意到一张图 ...

  5. Bzoj2395: [Balkan 2011]Timeismoney(最小乘积生成树)

    问题描述 每条边两个权值 \(x,y\),求一棵 \((\sum x) \times (\sum y)\) 最小的生成树 Sol 把每一棵生成树的权值 \(\sum x\) 和 \(\sum y\) ...

  6. bzoj2395 [Balkan 2011]Timeismoney(最小乘积生成树+计算几何)

    题意 每条边有两个权值\(c,t\),请求出一颗生成树,使得\(\sum c\times \sum t\)最小 题解 为什么生成树会和计算几何扯上关系-- 对于每棵树,设\(x=c,y=t\),我们可 ...

  7. BZOJ2395 [Balkan 2011]Timeismoney 【最小乘积生成树】

    题目链接 BZOJ2395 题意:无向图中每条边有两种权值,定义一个生成树的权值为两种权值各自的和的积 求权值最小的生成树 题解 如果我们将一个生成树的权值看做坐标,那么每一个生成树就对应一个二维平面 ...

  8. P5540-[BalkanOI2011]timeismoney|最小乘积生成树【最小生成树,凸壳】

    正题 题目链接:https://www.luogu.com.cn/problem/P5540 题目大意 给出\(n\)个点\(m\)条边边权是一个二元组\((a_i,b_i)\),求出一棵生成树最小化 ...

  9. bzoj 2395 [Balkan 2011]Timeismoney——最小乘积生成树

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2395 如果把 \( \sum t \) 作为 x 坐标,\( \sum c \) 作为 y ...

随机推荐

  1. Spring Boot配置文件yml讲解--行内对象的配置方式

    yml行内对象的配置方法,一般是采取 上面的缩进方式,我只想配置在一行怎么处?——

  2. Golang高阶:Golang1.5到Golang1.12包管理

    Golang1.5到Golang1.12包管理 1. 前言 Golang 是一门到如今有十年的静态高级语言了,2009年的时候算是正式推出了,然后到最近的一两年,2017-2018年的时候,突然直线上 ...

  3. 转:common.js 常用js公共函数库

    转自其他博主,自己开发备用 var h = {}; h.get = function (url, data, ok, error) { $.ajax({ url: url, data: data, d ...

  4. iOS - Scenekit3D引擎初探之 - 导出DAE文件(3Dmax为例)

    DAE文件格式是3D交互文件格式,一般用于多个图形程序之间交换数字数据,Autodesk专有并在COLLADA(COLLAborative Design Activity)基础上改进创建的XML框架的 ...

  5. git推送代码Gogs报401错误

    1.git push 报错:RPC failed; HTTP 401 curl 22 The requested URL returned error: 401 The remote end hung ...

  6. layui 后台分页

    <!DOCTYPE html> <html> <head> <meta charset="utf-8"> <title> ...

  7. VMware网络配置三种网络模式(桥接、NAT、Host-only)

    VMware网络配置三种网络模式(桥接.NAT.Host-only) 一.虚拟安装后三种网络模式显示 当安装好后,的“虚拟网络编辑器”中也存在三种模式,分别对应:桥接-VMnet0.Host-only ...

  8. django操作mysql

    连接mysql 1.安装pymysql 操作指令 : pymsql: pip install pymysql 2.导入库 在项目目录下的__init__.py文件中导入pymysql模块 加入以下两行 ...

  9. OpenStack环境搭建

    实验环境 CentOS-7-x86_64-Minimal-1708.iso openstack_N.tar.gz 创建虚拟机 controller部署 computer网络配置 OpenStack环境 ...

  10. 用js刷剑指offer(二叉树中和为某一值的路径)

    题目描述 输入一颗二叉树的根节点和一个整数,打印出二叉树中结点值的和为输入整数的所有路径.路径定义为从树的根结点开始往下一直到叶结点所经过的结点形成一条路径.(注意: 在返回值的list中,数组长度大 ...