[bzoj1941][Sdoi2010]Hide and Seek_KD-Tree
Hide and Seek bzoj-1941 Sdoi-2010
题目大意:给出平面上n个点,选出一个点,使得距离这个点的最远点曼哈顿距离减去距离这个点的最近非己点的曼哈顿距离最小。输出最小曼哈顿距离差。
注释:$1\le n \le 5\cdot 10^5$。
想法:这个题的话处理方法就是枚举所有点,然后对于每一个点求近邻点和“远”邻点。
其实区别就是估价函数的区别。近邻点对的估价函数就是SJY摆旗子的估价函数。
远邻点的话我们就另设一个估价函数,就表示当前点到这个矩形的可能最远曼哈顿距离。也就是矩形的4个顶点中距离当前点曼哈顿距离最远的曼哈顿距离。
然后枚举更新答案即可。
最后,附上丑陋的代码... ...
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 1000010
#define inf 0x7f7f7f7f
using namespace std;
int ans_min,ans_max,root,d;
struct Node
{
int c[2],p[2],maxn[2],minn[2];
}a[N];
inline bool cmp(const Node &a,const Node &b)
{
return a.p[d]==b.p[d]?a.p[d^1]<b.p[d^1]:a.p[d]<b.p[d];
}
void pushup(int k , int s)
{
a[k].minn[0] = min(a[k].minn[0] , a[s].minn[0]);
a[k].minn[1] = min(a[k].minn[1] , a[s].minn[1]);
a[k].maxn[0] = max(a[k].maxn[0] , a[s].maxn[0]);
a[k].maxn[1] = max(a[k].maxn[1] , a[s].maxn[1]);
}
int build(int l , int r , int now)
{
int mid = (l + r) >> 1;
d = now , nth_element(a + l , a + mid , a + r + 1 , cmp);
a[mid].minn[0] = a[mid].maxn[0] = a[mid].p[0];
a[mid].minn[1] = a[mid].maxn[1] = a[mid].p[1];
if(l < mid) a[mid].c[0] = build(l , mid - 1 , now ^ 1) , pushup(mid , a[mid].c[0]);
if(r > mid) a[mid].c[1] = build(mid + 1 , r , now ^ 1) , pushup(mid , a[mid].c[1]);
return mid;
}
int getdis_min(int qx,int qy,int x)
{
int ans=0;
if(qx<a[x].minn[0]) ans+=a[x].minn[0]-qx;
if(qx>a[x].maxn[0]) ans+=qx-a[x].maxn[0];
if(qy<a[x].minn[1]) ans+=a[x].minn[1]-qy;
if(qy>a[x].maxn[1]) ans+=qy-a[x].maxn[1];
return ans;
}
int getdis_max(int qx,int qy,int x)
{
int ans=0;
ans+=max(abs(qx-a[x].minn[0]),abs(qx-a[x].maxn[0]));
ans+=max(abs(qy-a[x].minn[1]),abs(qy-a[x].maxn[1]));
return ans;
}
void query_min(int qx,int qy,int x)
{
int dis_now=abs(qx-a[x].p[0])+abs(qy-a[x].p[1]);
if(dis_now) ans_min=min(dis_now,ans_min);
int ls=a[x].c[0],rs=a[x].c[1];
int dis_l=ls?getdis_min(qx,qy,ls):inf;
int dis_r=rs?getdis_min(qx,qy,rs):inf;
if(dis_l<dis_r)
{
if(ans_min>dis_l) query_min(qx,qy,ls);
if(ans_min>dis_r) query_min(qx,qy,rs);
}
else
{
if(ans_min>dis_r) query_min(qx,qy,rs);
if(ans_min>dis_l) query_min(qx,qy,ls);
}
}
void query_max(int qx,int qy,int x)
{
int dis_now=abs(qx-a[x].p[0])+abs(qy-a[x].p[1]);
ans_max=max(dis_now,ans_max);
int ls=a[x].c[0],rs=a[x].c[1];
int dis_l=ls?getdis_max(qx,qy,ls):0;
int dis_r=rs?getdis_max(qx,qy,rs):0;
if(dis_l>dis_r)
{
if(ans_max<dis_l) query_max(qx,qy,ls);
if(ans_max<dis_r) query_max(qx,qy,rs);
}
else
{
if(ans_max<dis_r) query_max(qx,qy,rs);
if(ans_max<dis_l) query_max(qx,qy,ls);
}
}
int main()
{
int n,ans=inf; cin >> n ;
for(int i=1;i<=n;i++) scanf("%d%d",&a[i].p[0],&a[i].p[1]);
root=build(1,n,1);
for(int i=1;i<=n;i++)
{
ans_min=inf,ans_max=0;
query_min(a[i].p[0],a[i].p[1],root), query_max(a[i].p[0],a[i].p[1],root);
ans=min(ans,ans_max-ans_min);
}
printf("%d\n",ans);
return 0;
}
小结:KD-Tree如果觉得复杂度不对就不对吧,这题就是驶过的。还有,那个tm30秒我还以为是每个店周围的斜正方形之内的点,然后一顿敲就是过不去样例,还没有样例解释...差评!
[bzoj1941][Sdoi2010]Hide and Seek_KD-Tree的更多相关文章
- [BZOJ1941][Sdoi2010]Hide and Seek
[BZOJ1941][Sdoi2010]Hide and Seek 试题描述 小猪iPig在PKU刚上完了无聊的猪性代数课,天资聪慧的iPig被这门对他来说无比简单的课弄得非常寂寞,为了消除寂寞感,他 ...
- BZOJ1941:[SDOI2010]Hide and Seek(K-D Tree)
Description 小猪iPig在PKU刚上完了无聊的猪性代数课,天资聪慧的iPig被这门对他来说无比简单的课弄得非常寂寞,为了消除寂寞感,他决定和他的好朋友giPi(鸡皮)玩一个更加寂寞的游戏- ...
- 【kd-tree】bzoj1941 [Sdoi2010]Hide and Seek
枚举每个点,计算离他最近的和最远的点. #include<cstdio> #include<cmath> #include<algorithm> using nam ...
- 【BZOJ1941】[Sdoi2010]Hide and Seek KDtree
[BZOJ1941][Sdoi2010]Hide and Seek Description 小猪iPig在PKU刚上完了无聊的猪性代数课,天资聪慧的iPig被这门对他来说无比简单的课弄得非常寂寞,为了 ...
- 【BZOJ-1941】Hide and Seek KD-Tree
1941: [Sdoi2010]Hide and Seek Time Limit: 16 Sec Memory Limit: 162 MBSubmit: 830 Solved: 455[Submi ...
- bzoj:1941: [Sdoi2010]Hide and Seek
1941: [Sdoi2010]Hide and Seek Time Limit: 16 Sec Memory Limit: 162 MBSubmit: 531 Solved: 295[Submi ...
- 【BZOJ1941】Hide and Seek(KD-Tree)
[BZOJ1941]Hide and Seek(KD-Tree) 题面 BZOJ 洛谷 题解 \(KD-Tree\)对于每个点搜一下最近点和最远点就好了 #include<iostream> ...
- BZOJ_1941_[Sdoi2010]Hide and Seek_KDtree
BZOJ_1941_[Sdoi2010]Hide and Seek_KDtree Description 小猪iPig在PKU刚上完了无聊的猪性代数课,天资聪慧的iPig被这门对他来说无比简单的课弄得 ...
- BZOJ 1941: [Sdoi2010]Hide and Seek(k-d Tree)
Time Limit: 16 Sec Memory Limit: 162 MBSubmit: 1712 Solved: 932[Submit][Status][Discuss] Descripti ...
随机推荐
- 一、Linux文件权限与目录配置
行文结构如下: 用户和用户组 Linux文件权限概念 Linux目录配置 重点回顾 1.用户与用户组 Linux是个多用户.多任务的系统,可能有多人同时使用这台机器进行工作,为了考虑每个人的隐私和工作 ...
- vue-resource 拦截器的使用
园友参考 https://www.cnblogs.com/lhl66/p/8022823.html vue-resource 拦截器使用在vue项目使用vue-resource的过程中,临时增加了一 ...
- 【洛谷4770/UOJ395】[NOI2018]你的名字(后缀数组_线段树合并)
题目: 洛谷4770 UOJ395 分析: 一个很好的SAM应用题-- 一句话题意:给定一个字符串\(S\).每次询问给定字符串\(T\)和两个整数\(l\).\(r\),求\(T\)有多少个本质不同 ...
- Objective-C——关联对象
动态语言 OC是一种动态语言,它的方法,对象的类型都是到运行的时候才能够确定的.所以这就使得OC存在了关联对象这一强大的机制. 关联对象 所谓关联对象,其实就是我们在运行时对一个已存在的对象上面绑定一 ...
- java学习笔记_网络
客户端 import java.io.*; import java.net.*; public class DailyAdviceClient { public void go() { try { S ...
- linux 卸载 mongo2.6
要求:linux 卸载 mongo2.6 版本:linux系统:Ubuntu 16.04 mongo: mongo 2.6.12 1. 查看安装的mongo版本和服务 # dpkg –l | gr ...
- BAT文件如何注释
1.BAT文件中如何注释: 1.:: 注释内容(第一个冒号后也可以跟任何一个非字母数字的字符) 2.rem 注释内容(不能出现重定向符号和管道符号) 3.echo 注释内容(不能出现重定向符号和管道符 ...
- Ajax——异步基础知识(三)
封装异步请求 1.将函数作为参数进行使用 2.因为获取数据是在一个注册事件中获取的,所以只有事件触发的时候才会调用此函数 <!DOCTYPE html> <html lang=&qu ...
- Effective Java中文版
译者序 序 前言 第一章 引言 第二章 创建和销毁对象 第1条:考虑用静态工厂方法代替构造函数 第2条:使用私有构造函数强化singleton属性 第3条:通过私有构造函数强化不可实例化属性 第4条: ...
- ajax不执行success的问题
有时候经常会遇到ajax请求后台,然后后台返回数据后,不触发ajax的success函数的问题,归根到底,这与ajax的参数设置dataType和后台的返回值的类型有关,现总结如下: 一.后台返回值的 ...