51nod 1213 二维曼哈顿距离最小生成树


第1行:1个数N,表示点的数量。(2 <= N <= 50000)
第2 - N + 1行:每行2个数,表示点的坐标(0 <= x, y <= 1000000)
输出N个点所组成的完全图的最小生成树的边权之和。
3
0 0
1 0
1 1
2
所以我们只要求一个点在其45°角的区域内离他最近的点就行了,而这可以用线段树或树状数组解决
我们以y轴正半轴往右偏45°角的区域为例:
点j在点i的这个区域要满足的条件是:
yj-xj>yi-xi
且xj>xi
那么我们将点以x为第一关键字,y为第二关键字,排序后倒序插入线段树
线段树的线段这一维是离散后的y-x,值是y+x
我们要求的是大于yi-xi的最小的y+x,而xj>xi这个条件已经由插入顺序满足了
这样我们成功的解决了这个区域的点
而其他区域的点我们可以通过坐标变换转移到这个区域
由于对称性,我们注意到其实只要求x轴或y轴正半轴所在的四个区域就行了
那么这个问题就这样解决了
#include <map>
#include <set>
#include <cmath>
#include <ctime>
#include <stack>
#include <queue>
#include <cstdio>
#include <memory>
#include <cctype>
#include <bitset>
#include <string>
#include <vector>
#include <climits>
#include <cstring>
#include <iostream>
#include <iomanip>
#include <algorithm>
#include <functional>
//#define FIN freopen("input.txt","r",stdin);
//#define FOUT freopen("output.txt","w+",stdout);
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const double eps=1e-8;
const double Pi=acos(-1.0);
const int N=50010;
struct point
{
int x,y,id;
bool operator<(const point p)const
{
return x!=p.x?x<p.x:y<p.y;
}
} p[N];
struct BIT
{
int min_val,pos;
void init()
{
min_val=INF;
pos=-1;
}
} bit[N];
int par[N];//并查集中父亲
int hight[N];//并查集树的高度
struct edge
{
int u,v,cost;
};
edge G[N<<2];//边集(边数)
int V,E;//顶点数和边数
int get_Manhadm_dis(point a,point b)
{
return abs(a.x-b.x)+abs(a.y-b.y);
}
void addedge(int u,int v,int w)
{
G[E].u=u;
G[E].v=v;
G[E++].cost=w;
}
int lowbit(int x)
{
return x&(-x);
}
void update(int x,int val,int pos)
{
for(int i=x; i>=1; i-=lowbit(i))
if(val<bit[i].min_val)
{
bit[i].min_val=val;
bit[i].pos=pos;
}
}
int ask(int x,int m)
{
int min_val=INF;
int pos=-1;
for(int i=x; i<=m; i+=lowbit(i))
if(bit[i].min_val<min_val)
{
min_val=bit[i].min_val;
pos=bit[i].pos;
}
return pos;
}
void make_edge()
{
int a[N],b[N];
for(int dir=0; dir<4; dir++)
{
if(dir==1||dir==3)
for(int i=0; i<V; i++)
swap(p[i].x,p[i].y);
else if(dir==2)
for(int i=0; i<V; i++)
p[i].x=-p[i].x;
sort(p,p+V);
for(int i=0; i<V; i++)
a[i]=b[i]=p[i].y-p[i].x;
sort(b,b+V);
int m=unique(b,b+V)-b;
for(int i=1; i<=m; i++)
bit[i].init();
for(int i=V-1;i>=0; i--)
{
int pos=lower_bound(b,b+m,a[i])-b+1;
int ans=ask(pos,m);
if(ans!=-1)
addedge(p[i].id,p[ans].id,get_Manhadm_dis(p[i],p[ans]));
update(pos,p[i].x+p[i].y,i);
}
}
}
//并查集初始化
void Init_union_find(int n)
{
for(int i=0; i<n; i++)
{
par[i]=i;
hight[i]=0;
}
}
//查询树的根
int find(int x)
{
if(par[x]==x)
return x;
else
return par[x]=find(par[x]);
}
//合并x和y所属的集合
void unite(int x,int y)
{
x=find(x);
y=find(y);
if(x==y)
return ;
if(hight[x]<hight[y])
par[x]=y;
else
{
par[y]=x;
if(hight[x]==hight[y])
hight[x]++;
}
}
//判断x和y是否属于同一个集合
bool same(int x,int y)
{
return find(x)==find(y);
}
bool cmp(const edge& a,const edge& b)
{
return a.cost<b.cost;
}
int kruskal()
{
sort(G,G+E,cmp);//按照edge.cost的顺序从小到大排列
Init_union_find(V);//并查集初始化
int ans=0;
for(int i=0; i<E; i++)
{
edge e=G[i];
if(!same(e.u,e.v))
{
unite(e.u,e.v);
ans+=e.cost;
}
}
return ans;
}
int main()
{
scanf("%d",&V);
for(int i=0; i<V; i++)
{
scanf("%d %d",&p[i].x,&p[i].y);
p[i].id=i+1;
}
E=0;
make_edge();
printf("%d\n",kruskal());
}
51nod 1213 二维曼哈顿距离最小生成树的更多相关文章
- [51nod1213]二维曼哈顿距离最小生成树
二维平面上有N个坐标为整数的点,点x1 y1同点x2 y2之间的距离为:横纵坐标的差的绝对值之和,即:Abs(x1 - x2) + Abs(y1 - y2)(也称曼哈顿距离).求这N个点所组成的完全图 ...
- 【POJ 3241】Object Clustering 曼哈顿距离最小生成树
http://poj.org/problem?id=3241 曼哈顿距离最小生成树模板题. 核心思想是把坐标系转3次,以及以横坐标为第一关键字,纵坐标为第二关键字排序后,从后往前扫.扫完一个点就把它插 ...
- 曼哈顿距离最小生成树 codechef Dragonstone
曼哈顿距离最小生成树 codechef Dragonstone 首先,对于每一个点来说有用的边只有它向它通过 x=0,y=0,y=x,y=-x 切出来的八个平面的最近点. 证明 我不会 反正当结论记住 ...
- LA 3662 Another Minimum Spanning Tree (曼哈顿距离最小生成树 模板)
题目大意: 曼哈顿最小距离生成树 算法讨论: 同上. 这回的模板真的准了. #include <iostream> #include <cstring> #include &l ...
- hdu6435 Problem J. CSGO标程讲解以及改正标程的一个错误(本来第一个样例过不了2333) 以及 poj2926 五维曼哈顿距离模板
比赛的时候抄poj2926的模板,但改不来啊orz #include <iostream> #include <cstdio> #include <cstring> ...
- POJ 3241 曼哈顿距离最小生成树 Object Clustering
先上几个资料: 百度文库有详细的分析和证明 cxlove的博客 TopCoder Algorithm Tutorials #include <cstdio> #include <cs ...
- POJ 3241Object Clustering曼哈顿距离最小生成树
Object Clustering Description We have N (N ≤ 10000) objects, and wish to classify them into several ...
- BZOJ 2735: 世博会 主席树+切比雪夫距离转曼哈顿距离
2735: 世博会 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 124 Solved: 51[Submit][Status][Discuss] D ...
- Unity插件系列之二维码
1.二维码常见的生成与识别途径 1.草料二维码 https://cli.im/text 2.在软件中实现生成和扫描二维码 使用zxing实现 zxing是一个用java写的开源项目,zxing.net ...
随机推荐
- Python使用django搭建web开发环境
安装 Python 去 Python 官方网站找到 Python 3 的下载地址,根据你的系统选择 32 位或者 64 位的安装包,下载好后双击安装即可. 检测是否安装完成 C:\WINDOWS\sy ...
- Highcharts:X轴分组堆叠图
在设计一个项目中的数据展示页面时.想要设计双X轴,一个轴显示须要的项.一个轴对这些项进行分组.效果如图: Highcharts自带双X轴展示方式.可是效果不是太理想.调整起来也会麻烦些 看到Highc ...
- ASP.NET MVC 学习笔记-7.自定义配置信息 ASP.NET MVC 学习笔记-6.异步控制器 ASP.NET MVC 学习笔记-5.Controller与View的数据传递 ASP.NET MVC 学习笔记-4.ASP.NET MVC中Ajax的应用 ASP.NET MVC 学习笔记-3.面向对象设计原则
ASP.NET MVC 学习笔记-7.自定义配置信息 ASP.NET程序中的web.config文件中,在appSettings这个配置节中能够保存一些配置,比如, 1 <appSettin ...
- 小胖说事35-----Terminating app due to uncaught exception 'CALayerInvalidGeometry', reason: 'CALayer posi
2011-06-11 15:19:17.167 ***[930:707] *** Terminating app due to uncaught exception 'CALayerInvalidGe ...
- 再理解HDFS的存储机制
再理解HDFS的存储机制 1. HDFS开创性地设计出一套文件存储方式.即对文件切割后分别存放: 2. HDFS将要存储的大文件进行切割,切割后存放在既定的存储块(Block)中,并通过预先设定的优化 ...
- JMeter Web测试快速入门教程
学习前的准备 学习本教程前,你的电脑上至少应该有Apache JMeter这款软件.如果你没有,点击此处下载. 当你点进去后,你会发现它是一个依赖Java的软件 因此如果你电脑上没有Java环境,你应 ...
- c++学习笔记之基础---类内声明线程函数的调用
近日需要将线程池封装成C++类,类名为Threadpool.在类的成员函数exec_task中调用pthread_create去启动线程执行例程thread_rounter.编译之后报错如下: spf ...
- Elasticsearch - 搜索类型与搜索位置
一.搜索类型: Elasticsearch同意用户选择其所希望的处理查询的方式. 由于存在一些不同的情形,对其使用不同的搜索类型才是合适的.为了控制查询的运行方式,我们能够在请求中使用search_t ...
- mysql--DML--视图,索引,事务,字符集,函数
Mysql中,表/列可以改名,database不能改名. phpMyAdmin似乎有这功能? 他是建新库,把所有表复制到新库,再删旧库完成的. 删除一个数据库: drop database 数据库名; ...
- bzoj1560: [JSOI2009]火星藏宝图
考虑到一个性质,A到B的距离一定不小于A到C再到B的距离,因为能够到达这三个点不可能构成锐角三角形 对于当前点的更新只需要找那些无法经过其它点再到当前点的点,相当于是一个y坐标单调减的上凸包,随便维护 ...