Position:


List

Description

  • 大意:求曼哈顿距离最小生成树上第k大(第n-k小)的边

Knowledge

参考曼哈顿距离最小生成树
http://blog.csdn.net/yjpyjp2014/article/details/52180707

Solution

分析:
1.由上面资料可知,每45度的范围只算一个最近点
2.看下面一种情况:
      
对于下面那个点,上面的点在它的右上方。而它又在上面的点的左下方,故我们只需计算左上和右上,即角[0,180]的范围,每45度一个,每个点4个,故边数不超过4×n。
3-1.为了写代码的方便,我们每次将点翻转,只计算[45,90]范围的点。翻转:①第一次不翻②第二次将一个点的x与y交换(即沿y=x翻转,原本在[0,45]的点翻到了[45,90])③将x=-x,沿x=0翻转,左上翻到右上④
3-2.ZYT的写法(可能容易理解一点):我们只需考虑在一块区域内的点,其他区域内的点可以通过坐标变换“移动”到这个区域内。为了方便处理,我们考虑在y轴向右45度的区域。在某个点A(x0,y0)的这个区域内的点B(x1,y1)满足x1≥x0且y1-x1>y0-x0。这里对于边界我们只取一边,但是操作中两边都取也无所谓。那么|AB|=y1-y0+x1-x0=(x1+y1)-(x0+y0)。在A的区域内距离A最近的点也即满足条件的点中x+y最小的点。因此我们可以将所有点按x坐标排序,再按y-x离散,用线段树或者树状数组维护大于当前点的y-x的最小的x+y对应的点(也就是维护区间最小值)。时间复杂度O(NlogN)。至于坐标变换,一个比较好处理的方法是第一次直接做(R1==R5);第二次沿直线y=x翻转,即交换x和y坐标(R2==R6);第三次沿直线x=0翻转,即将x坐标取相反数(R7==R3);第四次再沿直线y=x翻转(R8==R4)。注意只需要做4次,因为边是双向的。
             
4.对于点(x,y)如果有一个点(x1,y1)在它的(45,90),只要满足y1-x1>y-x(证明:即用一条斜率为45的直线,过这个点,比截距的大小y=x+b,y-x=b,平面图知识qaq~)
5.按x坐标排序,从后往前扫,值域(y-x)树状数组,w记录的是曼哈顿距离,pos记录第几个点。

Notice

1.这是一个反的树状数组,c[i]记录i~i+lowbit(i)-1
2.因为只要求x~m的最小值,用树状数组可以

Code

// <ObjectClustering.cpp> - 08/11/16 16:35:38
// This file is made by YJinpeng,created by XuYike's black technology automatically.
// Copyright (C) 2016 ChangJun High School, Inc.
// I don't know what this program is. #include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#define MOD 1000000007
#define INF 1e9
#define EPS 1e-10
using namespace std;
typedef long long LL;
const int MAXN=10010;
const int MAXM=100010;
inline int max(int &x,int &y) {return x>y?x:y;}
inline int min(int &x,int &y) {return x<y?x:y;}
inline int getint() {
register int w=0,q=0;register char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')q=1,ch=getchar();
while(ch>='0'&&ch<='9')w=w*10+ch-'0',ch=getchar();
return q?-w:w;
}
int n,m,k,tot,cnt,x,pos;
int f[MAXN],hs[MAXN],t[MAXN];
struct Point{
int x,y,id;
void read(){x=getint();y=getint();}
friend bool operator < (const Point &a,const Point &b) {
return a.x==b.x?a.y<b.y:a.x<b.x;
}
}p[MAXN];
struct Edge{
int u,v,w;
friend bool operator < (const Edge &a,const Edge &b){
return a.w<b.w;
}
}e[MAXN*4];
struct Bit{
int w,pos;
void Pre(){pos=-1,w=INF;}
}b[MAXN];
inline int Find(int x){return x==f[x]?x:f[x]=Find(f[x]);}
inline int Lowbit(int x){return x&(-x);}
int Dis(int a,int b){
return abs(p[a].x-p[b].x)+abs(p[a].y-p[b].y);
}
int Query(int x){
int M=INF,pos=-1;
for(int i=x;i<=m;i+=Lowbit(i))
if(b[i].w<M)M=b[i].w,pos=b[i].pos;
return pos;
}
void Update(int x,int M,int pos){
for(int i=x;i;i-=Lowbit(i))
if(M<b[i].w)b[i].w=M,b[i].pos=pos;
}
void Insert(int u,int v,int w){
e[++cnt].u=u;e[cnt].v=v;e[cnt].w=w;
}
void Caledge(){
sort(p+1,p+1+n);
for(int i=1;i<=n;i++)hs[i]=t[i]=p[i].y-p[i].x;
sort(hs+1,hs+1+n);
m=unique(hs+1,hs+1+n)-hs;
for(int i=1;i<=m;i++)b[i].Pre();
for(int i=n;i>=1;i--){
x=lower_bound(hs+1,hs+1+m,t[i])-hs+1;
pos=Query(x);
if(pos!=-1)Insert(p[i].id,p[pos].id,Dis(i,pos));
Update(x,p[i].x+p[i].y,i);
}
}
void Kruskal(){
sort(e+1,e+1+cnt);tot=0;
for(int i=1;i<=n;i++)f[i]=i;
for(int i=1;i<=cnt;i++){
f[e[i].u]=Find(f[e[i].u]);
f[e[i].v]=Find(f[e[i].v]);
if(f[e[i].u]!=f[e[i].v])f[f[e[i].u]]=f[e[i].v],tot++;
if(tot==k){printf("%d",e[i].w);break;}
}
}
int main()
{
freopen("ObjectClustering.in","r",stdin);
freopen("ObjectClustering.out","w",stdout);
n=getint();k=n-getint();
for(int i=1;i<=n;i++)p[i].read(),p[i].id=i;
for(int j=1;j<=4;j++){
if(j==2||j==4)
for(int i=1;i<=n;i++)swap(p[i].x,p[i].y);
else if(j==3)
for(int i=1;i<=n;i++)p[i].x=-p[i].x;
Caledge();
}
Kruskal();
return 0;
}

感谢

Friend ZYT的大力支持

【Poj3241】Object Clustering的更多相关文章

  1. 【poj3241】 Object Clustering

    http://poj.org/problem?id=3241 (题目链接) MD被坑了,看到博客里面说莫队要写曼哈顿最小生成树,我就写了一个下午..结果根本没什么关系.不过还是把博客写了吧. 转自:h ...

  2. 【POJ 3241】Object Clustering 曼哈顿距离最小生成树

    http://poj.org/problem?id=3241 曼哈顿距离最小生成树模板题. 核心思想是把坐标系转3次,以及以横坐标为第一关键字,纵坐标为第二关键字排序后,从后往前扫.扫完一个点就把它插 ...

  3. 前端开发者进阶之ECMAScript新特性【一】--Object.create

    Object.create(prototype, descriptors) :创建一个具有指定原型且可选择性地包含指定属性的对象 参数:prototype 必需.  要用作原型的对象. 可以为 nul ...

  4. 【转】object标签和embed标签

    我们现在大部分人做网页,都是直接用DW插入flash,而且DW也是所见即所得,直接生成了相应的flash显示代码.可是我们又有多少人了解这些直接由DW生成的代码呢?其实我接触flash player标 ...

  5. 【JavaScript】Object.observe()带来的数据绑定变革

    Object.observe()带来的数据绑定变革 引言 一场变革即将到来.一项Javascript中的新特性将会改变你对于数据绑定的所有认识.它也将改变你所使用的MVC库观察模型中发生的修改以及更新 ...

  6. 【原】Object 异常静态

    //所有的类都继承Object类: Object a=; Object b="ddfasfda"; //正常情况下,都会省略掉:Object:但实际上是存在的: class tes ...

  7. 【JavaScript】Object.prototype.toString.call()进行类型判断

    权声明:本文为博主原创文章,未经博主允许不得转载. op = Object.prototype, ostring = op.toString, ... function isFunction(it)  ...

  8. 为什么要用lock 【readonly】object?为什么不要lock(this)?

    一. 为什么要用lock,lock了什么? 当我们使用线程的时候,效率最高的方式当然是异步,即各个线程同时运行,其间不相互依赖和等待.但当不同的线程都需要访问某个资源的时候,就需要同步机制了.也就是说 ...

  9. 【vue】---Object.defineProperty基本使用---【巷子】

    1.object.defineProperty 给一个对象定义一个新的属性或者在修改一个对象现有的属性,并返回这个对象 语法: Object.defineProperty(参数1,参数2,参数3) 参 ...

随机推荐

  1. SpringBoot中如何使用jpa和jpa的相关知识总结

    jpa常用的注解: 注解 解释 @Entity 声明类为实体或表. @Table 声明表名. @Basic 指定非约束明确的各个字段. @Embedded 指定类或它的值是一个可嵌入的类的实例的实体的 ...

  2. JS函数assign

    Object函数提供了一个叫做assign的函数,用来合并多个对象. Object.assign(...): 你可以传递多个对象给该函数,这些对象中的自有且可枚举的属性,会被拷贝给第一个对象. var ...

  3. node事件循环和消息队列简单分析

    node的好处毋庸置疑,事件驱动,异步非阻塞I/O,以及处理高并发的能力深入人心,因此大家喜欢用node做一些小型后台服务或者作为中间层和其他服务配合完成一些大型应用场景. 什么是异步? 异步和同步应 ...

  4. 洛谷——P2171 Hz吐泡泡

    P2171 Hz吐泡泡 题目描述 这天,Hz大大心血来潮,吐了n个不同的泡泡玩(保证没有重复的泡泡).因为他还要写作业,所以他请你帮他把这些泡泡排序成树(左子树<=根<右子树).输出它的后 ...

  5. keepalived(nginx的高可用)安装文档

    1. 安装环境 su - root yum -y install kernel-devel* yum -y install openssl-* yum -y install popt-devel yu ...

  6. python下载网络文件

    python下载网络文件 制作人:全心全意 下载图片 #!/usr/bin/python #-*- coding: utf-8 -*- import requests url = "http ...

  7. msdn的原版windows下载地址链接

    http://msdn.itellyou.cn/ 所有版本的下载地址 进去点左边操作系统

  8. npm 使用教程

    链接----------------------------------npm官网npm淘宝镜像 安装包----------------------------------npm install -g ...

  9. Ubuntu挂载硬盘,修改卷标

    Ubuntu挂载硬盘,修改卷标转载2016-03-06 17:03:21标签:ubuntu Ubuntu不像windows,硬盘插入电脑不会自动读取硬盘 数据,需要把硬盘挂载到文件夹上,然后才能访问硬 ...

  10. idea 编译级别的设置

    File->Settings Project Structure