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\) 由这两个性 ...
随机推荐
- UML和模式应用5:细化阶段(8)---逻辑架构和UML包图
1.前言 本章是从面向分析的工作过度到软件设计 典型的OO系统设计的基础是若干架构层,如UI层.应用逻辑(领域)层 本章简要考察逻辑分层架构和相关UML表示法 2.逻辑架构和层 逻辑架构 逻辑架构是软 ...
- c# webbrowser控件内核版本强制修改
int BrowserVer, RegVal; // get the installed IE version using (WebBrowser Wb = new WebBrowser()) Bro ...
- #ifndef详解
#ifndef 是"if not defined"的简写,是预处理功能(宏定义.文件包含.条件编译)当中的条件编译,可以根据是否已经定义了一个变量来进行分支选择,其作用是: 1.防 ...
- Shell 中test 单中括号[] 双中括号[[]] 的区别
Shell test 单中括号[] 双中括号[[]] 的区别 在写Shell脚本的时候,经常在写条件判断语句时不知道该用[] 还是 [[]],首先我们来看他们的类别: $type [ [[ test ...
- pom.xml的第一行报错
第一行:<project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.or ...
- Tomcat:在centos中做成自启动服务
# 创建一个启动脚本文件,脚本内容见下 vi /etc/init.d/tomcat #!/bin/bash # /etc/rc.d/init.d/tomcat # init script for to ...
- 统一异常处理@ExceptionHandler
异常处理功能中用到的注解是:@ExceptionHandler(异常类型.class). 这个注解的功能是:自动捕获controller层出现的指定类型异常,并对该异常进行相应的异常处理. 比如我要在 ...
- 升级 php composer 版本
在执行 composer update 时,报错 You made a reference to a non-existent script @php artisan package:discover ...
- hdu4812 逆元+树分治
逆元链接:https://www.cnblogs.com/zzqc/p/7192436.html 经典的树分治题 #pragma comment("linker,"/STACK:1 ...
- c++ primer 学习杂记2【派生类到基类转换的可访问性】
参考: http://blog.csdn.net/rehongchen/article/details/7930853 http://blog.csdn.net/ming_road/article/d ...