P5540-[BalkanOI2011]timeismoney|最小乘积生成树【最小生成树,凸壳】
正题
题目链接:https://www.luogu.com.cn/problem/P5540
题目大意
给出\(n\)个点\(m\)条边边权是一个二元组\((a_i,b_i)\),求出一棵生成树最小化
\]
的情况下最小化\(\sum_{e\in T}a_e\)
\(1\leq n\leq 200,1\leq m\leq 10^4\)
解题思路
这种带乘积的可以维护凸壳,对于一棵生成树\(T\)我们视为一个\((\sum_{e\in T}a_e,\sum_{e\in T}b_i)\)的点,这样我们打答案一定在下凸壳上。
可以用一种分治求凸壳的方法,我们先找出下凸壳的两个端点(\(x\)最小的和\(y\)最小的)记为\(A,B\),然后找到一个在\(A\)与\(B\)的连边下面的一个最凸的点\(C\)(可以视为最大化\(S_{\bigtriangleup ACB}\),这样\(C\)一定在凸壳上),然后分治下去做\(\vec{AC}\)和\(\vec{CB}\)。
考虑怎么求这个\(C\),就是最大化\(\vec{AC}\times \vec{CB}\)
\]
\]
然后就是相当于最小化\(x_C(y_B-y_A)+y_C(x_A-x_B)\),拿这个当边权跑就可以跑出\(C\)了。
然后时间复杂度据说是\(O(m\log m\sqrt{\ln n!})\)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int N=210,M=1e4+10;
struct node{
ll x,y,w,id;
}e[M];
struct point{
ll x,y;
point(ll xx=0,ll yy=0)
{x=xx;y=yy;return;}
}ans;
ll n,m,x[M],y[M],a[M],b[M],fa[N];
point operator-(point x,point y)
{return point(x.x-y.x,x.y-y.y);}
ll operator*(point x,point y)
{return x.x*y.y-x.y*y.x;}
bool cmp(node x,node y)
{return (x.w==y.w)?(a[x.id]<a[y.id]):(x.w<y.w);}
ll find(ll x)
{return (fa[x]==x)?x:(fa[x]=find(fa[x]));}
point Kruskal(){
ll cnt=0;point res=0;
for(ll i=1;i<=n;i++)fa[i]=i;
sort(e+1,e+1+m,cmp);
for(ll i=1;i<=m;i++){
ll x=find(e[i].x),y=find(e[i].y);
if(x==y)continue;
fa[x]=y;cnt++;
res.x+=a[e[i].id];
res.y+=b[e[i].id];
if(cnt==n-1)break;
}
if(res.x*res.y<ans.x*ans.y)ans=res;
else if(res.x*res.y==ans.x*ans.y&&res.x<ans.x)
ans=res;
return res;
}
void solve(point A,point B){
for(ll i=1;i<=m;i++)
e[i]=(node){x[i],y[i],(B.x-A.x)*b[i]+(A.y-B.y)*a[i],i};
point C=Kruskal();
if((C-A)*(B-A)<=0)return;
solve(A,C);solve(C,B);
}
signed main()
{
scanf("%lld%lld",&n,&m);
for(ll i=1;i<=m;i++){
scanf("%lld%lld%lld%lld",&x[i],&y[i],&a[i],&b[i]);
x[i]++;y[i]++;
}
ans.x=ans.y=1e9;
for(ll i=1;i<=m;i++)e[i]=(node){x[i],y[i],a[i],i};
point A=Kruskal();
for(ll i=1;i<=m;i++)e[i]=(node){x[i],y[i],b[i],i};
point B=Kruskal();
solve(A,B);
printf("%lld %lld\n",ans.x,ans.y);
return 0;
}
P5540-[BalkanOI2011]timeismoney|最小乘积生成树【最小生成树,凸壳】的更多相关文章
- 洛谷 P5540 - [BalkanOI2011] timeismoney | 最小乘积生成树(最小生成树)
洛谷题面传送门 大概是一个比较 trivial 的小 trick?学过了就不要忘了哦( 莫名奇妙地想到了 yyq 的"hot tea 不常有,做过了就不能再错过了" 首先看到这种二 ...
- bzoj2395[Balkan 2011]Timeismoney最小乘积生成树
所谓最小乘积生成树,即对于一个无向连通图的每一条边均有两个权值xi,yi,在图中找一颗生成树,使得Σxi*Σyi取最小值. 直接处理问题较为棘手,但每条边的权值可以描述为一个二元组(xi,yi),这也 ...
- 【BZOJ2395】【Balkan 2011】Timeismoney 最小乘积生成树
链接: #include <stdio.h> int main() { puts("转载请注明出处[辗转山河弋流歌 by 空灰冰魂]谢谢"); puts("网 ...
- 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 ...
- bzoj 2395 Timeismoney —— 最小乘积生成树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2395 参考博客:https://www.cnblogs.com/autsky-jadek/p ...
- Luogu5540 最小乘积生成树
Luogu5540 最小乘积生成树 题目链接:洛谷 题目描述:对于一个\(n\)个点\(m\)条边的无向连通图,每条边有两个边权\(a_i,b_i\),求使\((\sum a_i)\times (\s ...
- HDU5697 刷题计划 dp+最小乘积生成树
分析:就是不断递归寻找靠近边界的最优解 学习博客(必须先看这个): 1:http://www.cnblogs.com/autsky-jadek/p/3959446.html 2:http://blog ...
- 【算法】最小乘积生成树 & 最小乘积匹配 (HNOI2014画框)
今天考试的时候果然题目太难于是我就放弃了……转而学习了一下最小乘积生成树. 最小乘积生成树定义: (摘自网上一篇博文). 我们主要解决的问题就是当k = 2时,如何获得最小的权值乘积.我们注意到一张图 ...
随机推荐
- webGis概念
参考:https://blog.csdn.net/qq_36375770/article/details/80077533 参考:https://blog.csdn.net/BuquTianya/ar ...
- 本文详细阐述如何用C#创建COM组件,并能用VC6.0等调用。
本文详细阐述如何用C#创建COM组件,并能用VC6.0等调用. 附:本文适用任何VS系列工具. 在用C#创建COM组件时,一定要记住以下几点: 1.所要导出的类必须为公有: 2.所有属性.方法也必须为 ...
- Wpf程序显示在任务栏
后台代码如下: using System; using System.Collections.Generic; using System.Drawing; using System.IO; using ...
- win10下安装mysql5.7-建议安装-自定义位置
所需文件: 文件名字:mysql-installer-community-5.7.13.0.msi 开始安装: 01:勾选: 02:选择下一步 03:选择Custom安装 04:选择下一步 05:选择 ...
- 由struts2中配置使用servlet引发的思考和复习
Struts2拦截器到底拦截了什么? 关于struts2中的拦截器,首先再次理解其实只能过滤其中访问的action的映射!再者,因为struts中的action其实就是起到替代servlet作用的,所 ...
- Learning ROS: Packaging your ROS project as a snap
Snaps are containerised software packages that are simple to create and install. They auto-update an ...
- SQLServer数据实时同步PostgreSQL
SQLServer数据实时同步至PostgreSQL 前言: 为迎合工作需求有时候传送的数据保存在SQLServer中但由于工作需要需要保存到PostgreSQL中进行处理,本文主要通过在SQLSer ...
- 经典多级时间轮定时器(C语言版)
经典多级时间轮定时器(C语言版) 文章目录 经典多级时间轮定时器(C语言版) 1. 序言 2. 多级时间轮实现框架 2.1 多级时间轮对象 2.2 时间轮对象 2.3 定时任务对象 2.4 双向链表 ...
- 硕盟 TYPE C转HDMI+VGA+USB3.0+PD3.0四口扩展坞
硕盟SM-T54是一款USB-C 四合一多功能扩展坞,支持四口同时使用,您可以将含有USB 3.1协议的电脑主机,通过此产品连接到具有HDMI或VGA的显示器.电视机或其他显示设备.产品可以接入硬盘. ...
- 马哈鱼数据血缘分析器分析case-when语句
马哈鱼数据血缘分析器是一个分析数据血缘关系的平台,可以在线直接递交 SQL 语句进行分析,也可以选择连接指定数据库获取 metadata.从本地上传文件目录.或从指定 git 仓库获取脚本进行分析. ...