BZOJ.2177.曼哈顿最小生成树(Kruskal)
\(Solution\)
参考
对于每个点,向唯一有可能与它形成MST的8个点连边,由于是双向单边,所以每个点最多连出4条边(证明见blog)
怎么找到一个区域内最近的点?
只考虑y轴右侧45°的区域,其余部分可以通过坐标变换移动到这一区域
设当前点P(x0,y0),这一区域一点P1(x1,y1),满足x1>x0 && y1-x1>y0-x0
那么dis(AB)= y1-y0+x1-x0 = x1+y1-(x0+y0)
这样这一区域内dis(AB)最小的点即 在满足之前条件的点中 x1+y1最小的点
先对x排序;then由于是只查询某个后缀,可以对y-x离散后用树状数组找到最小的x+y
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define lb(x) (x&-x)
#define gc() getchar()
#define dis(i,j) (std::abs(p[i].x-p[j].x)+std::abs(p[i].y-p[j].y))
const int N=5e4+5;
int n,Enum,fa[N],A[N],ref[N];
struct Point
{
int x,y,id;
inline bool operator <(const Point &a)const {return x==a.x?y<a.y:x<a.x;}
}p[N];
struct Bit_Tree
{
int val[N],pos[N];
inline void Init() {memset(val,0x7f,sizeof val);}
void Update(int p,int v,int vp)
{
for(;p;p-=lb(p))
if(val[p]>v) val[p]=v,pos[p]=vp;
}
int Query(int p,int n)
{
int mn=1<<30,res=-1;
for(;p<=n;p+=lb(p))
if(mn>val[p]) mn=val[p],res=pos[p];
return res;
}
}t;
struct Edge
{
int to,fr,val;
inline bool operator <(const Edge &a)const {return val<a.val;}
}e[N<<2];
inline void AddEdge(int u,int v,int w)
{
e[++Enum].to=v, e[Enum].fr=u, e[Enum].val=w;
}
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
int Getf(int x)
{
return x==fa[x]?x:fa[x]=Getf(fa[x]);
}
int MST()
{
for(int dir=0;dir<4;++dir)//坐标变换
{//Time1.不变
if(dir==1||dir==3)
for(int i=1;i<=n;++i) std::swap(p[i].x,p[i].y);//T2.T4.关于y=x对称
else if(dir==2)
for(int i=1;i<=n;++i) p[i].x=-p[i].x;//T3.关于x=0对称
std::sort(p+1,p+1+n);
t.Init();
for(int i=1;i<=n;++i) A[i]=ref[i]=p[i].y-p[i].x;
std::sort(A+1,A+1+n);
int m=1; ref[1]=A[1];
for(int i=2;i<=n;++i)
if(A[i]!=A[i-1]) ref[++m]=A[i];
for(int pos,res,i=n;i;--i)
{
pos=std::lower_bound(ref+1,ref+1+m,p[i].y-p[i].x)-ref;
if((res=t.Query(pos,m))!=-1) AddEdge(p[i].id,p[res].id,dis(i,res));
t.Update(pos,p[i].x+p[i].y,i);
}
}
std::sort(e+1,e+1+Enum);
for(int i=1;i<=n;++i) fa[i]=i;
int res=0;
for(int r1,r2,i=1,k=0;i<=Enum;++i)
{
r1=Getf(e[i].fr),r2=Getf(e[i].to);
if(r1!=r2)
{
res+=e[i].val, fa[r1]=r2;
if(++k==n-1) break;
}
}
return res;
}
int main()
{
n=read();
for(int i=1;i<=n;++i) p[i].x=read(),p[i].y=read(),p[i].id=i;
printf("%d",MST());
return 0;
}
BZOJ.2177.曼哈顿最小生成树(Kruskal)的更多相关文章
- BZOJ 2177: 曼哈顿最小生成树
Sol 考了好几次曼哈顿最小生成树,然而一直不会打...这次终于打出来了...神tm调试了2h...好蛋疼... 首先曼哈顿最小生成树有个结论就是讲它每45度分出一个象限,对于每个点,只与每个象限中离 ...
- 【BZOJ-2177】曼哈顿最小生成树 Kruskal + 树状数组
2177: 曼哈顿最小生成树 Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 190 Solved: 77[Submit][Status][Discu ...
- poj 3241 Object Clustering (曼哈顿最小生成树)
Object Clustering Time Limit: 2000MS Memory Limit: 131072K Total Submissions: 2640 Accepted: 806 ...
- POJ3241 Object Clustering 曼哈顿最小生成树
题意:转换一下就是求曼哈顿最小生成树的第n-k条边 参考:莫涛大神的论文<平面点曼哈顿最小生成树> /* Problem: 3241 User: 96655 Memory: 920K Ti ...
- 模板——最小生成树kruskal算法+并查集数据结构
并查集:找祖先并更新,注意路径压缩,不然会时间复杂度巨大导致出错/超时 合并:(我的祖先是的你的祖先的父亲) 找父亲:(初始化祖先是自己的,自己就是祖先) 查询:(我们是不是同一祖先) 路径压缩:(每 ...
- 最小生成树——Kruskal与Prim算法
最小生成树——Kruskal与Prim算法 序: 首先: 啥是最小生成树??? 咳咳... 如图: 在一个有n个点的无向连通图中,选取n-1条边使得这个图变成一棵树.这就叫“生成树”.(如下图) 每个 ...
- 【转】最小生成树——Kruskal算法
[转]最小生成树--Kruskal算法 标签(空格分隔): 算法 本文是转载,原文在最小生成树-Prim算法和Kruskal算法,因为复试的时候只用到Kruskal算法即可,故这里不再涉及Prim算法 ...
- BZOJ 1016: [JSOI2008]最小生成树计数( kruskal + dfs )
不同最小生成树中权值相同的边数量是一定的, 而且他们对连通性的贡献是一样的.对权值相同的边放在一起(至多10), 暴搜他们有多少种方案, 然后乘法原理. ----------------------- ...
- BZOJ.1016.[JSOI2008]最小生成树计数(Matrix Tree定理 Kruskal)
题目链接 最小生成树有两个性质: 1.在不同的MST中某种权值的边出现的次数是一定的. 2.在不同的MST中,连接完某种权值的边后,形成的连通块的状态是一样的. \(Solution1\) 由这两个性 ...
随机推荐
- Linux设备驱动之Ioctl控制【转】
转自:http://www.cnblogs.com/geneil/archive/2011/12/04/2275372.html 大部分驱动除了需要具备读写设备的能力之外,还需要具备对硬件控制的能力. ...
- 【bzoj1901】dynamic ranking(带修改主席树/树套树)
题面地址(权限题) 不用权限题的地址 首先说说怎么搞带修改主席树? 回忆一般的kth问题,我们的主席树求的是前缀和,这样我们在目标区间的左右端点的主席树差分下就能求出kth. 那么我们如何支持修改操作 ...
- adb启动和停止android app方法
一.启动app adb shell am start 包名/MainActivity 上面涉及到的包名及mainactivity获取办法 方法一: 1.adb shell 2.dumpsys act ...
- oracle 数据库链路
Oracle数据库链路的建立和使用 一.数据库链路的建立: 1.数据库链路的建立语法一般是:CREATE DATABASE LINK[db_link_name] CONNECT TO [user_na ...
- 使用JDBC连接数据库报“找不到驱动程序”错误解决
1.jre安装不成功: 2.jre中没有数据库驱动(D:\jre\lib\ext 中查看)
- [工具/PC]计算机中丢失libiconv-2.dll,丢失libintl-8.dll,无法定位程序输入点libiconv于动态链接库libiconv-2.dll上问题解决方法
CodeBlocks 1. 背景,为了学习C语言,在win系统上下载了codeBlock,先简单介绍下:Code::Blocks 是一个开放源码的全功能的跨平台C/C++集成开发环境. Code::B ...
- LeetCode(55): 跳跃游戏
Medium! 题目描述: 给定一个非负整数数组,你最初位于数组的第一个位置. 数组中的每个元素代表你在该位置可以跳跃的最大长度. 判断你是否能够到达最后一个位置. 示例 1: 输入: [2,3,1, ...
- LiteQuery MAX(Integer)、MAX(String) 判断是否返回值
unit Unit6; interface uses Winapi.Windows, Winapi.Messages, System.SysUtils, System.Variants, System ...
- pyhon----模块导入
正常情况 报错(两个导入都报错)
- python接口自动化测试八:更新Cookies、session保持会话
s = requests.session() # 此方法只适用于网站是cookies这种,网站是token的没用 # 这样做的好处就是可以保存cookies并保持会话,不用每次都去获取.传参 Toke ...