KD-Tree


  一开始看错题了aaarticlea/jpeg;base64,/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQH/2wBDAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQEBAQH/wAARCAAYABoDASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD+u/8AbK/aO134M6JD4J8C+JfBngjx/rnw8+L3xl1T4ifEbQvEfirwD8H/ANnX4ES+BL39oD43674D8OazomufEe68GReP/BPhzwf8NdG17TpvE/jHxNo82ra3aeC7TxdPX5m6L+1bolzB8Srn4n/tf/8ABSr4SfFj4d+A/B/jmLwPqvwo/wCCfEmveM4viFqvwp8IeFfB/wAMvg98Pfgt8btEsfG2t+IfjT8GLHS/gr8abyx+NXhY/En4M6p4xMGm65H4wf6U/bd8CeLNW/a3+HGvaa8F9Yaz8Fr74cWdlDqxi1P+2b34p2+sm3vNIJ2jS9aC6JJ/aF0/lXV/pAjLCS0c18Uat4Y0HxLJ+21qPjz9nn9tP4123wXvfH/wBufix41/al/Yq8A2P7Pvw3+HHxR8P/HXwl4v+DdpZ698H/EngDTv7Z8D/CX48eDPF/i3wL448atpuj/Dd08S6ybLXfDbdlbLKOX5Lk+bzx8FVzuOOTwKlT/d4HB1asPr7lze9Kq4TUIqna0K0aspzXMfA5RxLLM+KuJMipZU/q/Dbypf2h/aHuzzPGzr1J4J5fGCX+w4GPOpLFJpzUuWUleX3d+wB+3xrXx4+K3xg/ZS+M+mahpfx0+DFi3iCCbW9M8B6Tq3xH+Gel+NLj4ceIvEOraN8NvjF8Y/A1j4w+GXxF0jV/h58TLfwZ8QNe8PWHjENo1xJoXxItPiF8KbH9ZtsP8Az7f+Rn/xr8Sv2dfCukwftT/s3fGTxJJ4rk+I/iX4IaAureKvizaeH0+McGkfEn4QaT4u1H4Y/EDXtC8N+FdE1K5XxNa+G72+i8N6Jonh1PGcTaTo3gi0lV5a/cHzZPRv++1/+JrfG5NWyOeFUcc8Rh8xy3BZjhpyS5lTruumuZcnMuaMmnZxakpQlODVR78H8U1OLKOeSw1NZZiMkz3H5HmGESjKMcXgZ8s3H2lRNXvdpd03GMm2/hPx5+zj4j+K/wC0tB4q8bWsqfDXw5beDtS0G60Lxj4i8M+IJte8OXU+uaPFBe+FNb0TXNLuNF8V21rq8phujpWpeH408LavFd6LdazpVer/ABH/AGPf2V/jT4+s/ib8WPgH8M/HvjCxg0S0HiLxR4R0nVJdSsdA1ObV/Dlv4itbiJ7XxPD4V1PF74cHiW0vRoupvJqmiQ2srzzEory8XmGJx9DAxruHLluAyvD4RU4cqhChicSoOV23Kbm3WnJ6zrSlOTblK/Rw3w9gMjqZ7PBzxUpZpnuPzrGSxFf2rnjMRODqJPkjy0FdKFF3SjGEXKUVI8N/ax+BmnXFjqPxJ8AeEdUvPHviXxp4d1DxLrOkz69qutXNlYeE7TwnZQW8DX14+gaFYf8ACP6Hu0jwxaWum2Wpvrvi2VE13UvHutzfYGj3+uvpGlvrjWdhrLadYtq9hb32nT29lqjWyHULSCczkzQ2115sMUpOZERXJJ3ZKK+hjj8TPJsswVSUatHBVcb9WVSCk6ca0MHKcI3dow5qalGEUoqUpys5Nt/HVMjwVDjLifH4eeJw9fNMBw5WxvsKqhCtWp/22o1XFU9JyUrSd2mlFJJpt//Z" alt="" />

  其实是:给定n个点,从中找一个点,使得其他所有点到它距离的最大值与最小值之差最小。

  利用KD-Tree暴力求出每个点的答案(找离它最近的点以及最远的点(当然只关心距离))

  然后……两个过程分开写……

  注意一下最近的点的距离不能是0(然而我一开始用 if (o==tmp) return INF; 就WA了……)(这里o是当前搜索到的点,tmp是枚举的起始点)

  

 /**************************************************************
Problem: 1941
User: Tunix
Language: C++
Result: Accepted
Time:1484 ms
Memory:16900 kb
****************************************************************/ //BZOJ 1941
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
int r=,v=; char ch=getchar();
for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-;
for(; isdigit(ch);ch=getchar()) v=v*-''+ch;
return r*v;
}
const int N=,INF=1e9;
/*******************template********************/ struct node{
int d[],mn[],mx[],l,r;
int& operator [] (int x){return d[x];}
void read(){d[]=getint(); d[]=getint();}
}t[N];
int n,m,D,root,tmp,ask_mx,ask_mn;
bool operator < (node a,node b){return a[D]<b[D];}
#define L t[o].l
#define R t[o].r
#define mid (l+r>>1)
void Push_up(int o){
F(i,,){
t[o].mn[i]=min(t[o].mn[i],min(t[L].mn[i],t[R].mn[i]));
t[o].mx[i]=max(t[o].mx[i],max(t[L].mx[i],t[R].mx[i]));
}
} int build(int l,int r,int dir){
D=dir;
nth_element(t+l,t+mid,t+r+);
int o=mid;
F(i,,) t[o].mn[i]=t[o].mx[i]=t[o][i];
if (l<mid) L=build(l,mid-,dir^);
if (mid<r) R=build(mid+,r,dir^);
Push_up(o);
return o;
} int dis(int a){return abs(t[a][]-t[tmp][])+abs(t[a][]-t[tmp][]);}
int calc_mn(int o){
if (!o) return INF;
int ans=;
F(i,,) ans+=max(,t[o].mn[i]-t[tmp][i]);
F(i,,) ans+=max(,t[tmp][i]-t[o].mx[i]);
return ans;
} void query_mn(int o){
if (!o) return;
int dl=calc_mn(L),dr=calc_mn(R),d0=dis(o);
if (d0) ask_mn=min(ask_mn,d0);
if (dl<dr){
if (dl<ask_mn) query_mn(L);
if (dr<ask_mn) query_mn(R);
}else{
if (dr<ask_mn) query_mn(R);
if (dl<ask_mn) query_mn(L);
}
} int calc_mx(int o){
if (!o) return -INF;
int ans=;
F(i,,) ans+=max(abs(t[o].mn[i]-t[tmp][i]),abs(t[o].mx[i]-t[tmp][i]));
return ans;
}
void query_mx(int o){
if (!o) return;
int dl=calc_mx(L),dr=calc_mx(R),d0=dis(o);
ask_mx=max(ask_mx,d0);
if (dl>dr){
if (dl>ask_mx) query_mx(L);
if (dr>ask_mx) query_mx(R);
}else{
if (dr>ask_mx) query_mx(R);
if (dl>ask_mx) query_mx(L);
}
} int main(){
#ifndef ONLINE_JUDGE
freopen("1941.in","r",stdin);
freopen("1941.out","w",stdout);
#endif
n=getint();
F(i,,n) t[i].read();
F(i,,) t[].mn[i]=INF,t[].mx[i]=-INF;
root=build(,n,);
int ans=INF;
F(i,,n){
ask_mn=INF; ask_mx=-INF;
tmp=i;
query_mn(root);
query_mx(root);
// printf("%d %d\n",ask_mx,ask_mn);
ans=min(ans,ask_mx-ask_mn);
}
printf("%d\n",ans);
return ;
}

1941: [Sdoi2010]Hide and Seek

Time Limit: 16 Sec  Memory Limit: 162 MB
Submit: 385  Solved: 213
[Submit][Status][Discuss]

Description


猪iPig在PKU刚上完了无聊的猪性代数课,天资聪慧的iPig被这门对他来说无比简单的课弄得非常寂寞,为了消除寂寞感,他决定和他的好朋友
giPi(鸡皮)玩一个更加寂寞的游戏---捉迷藏。
但是,他们觉得,玩普通的捉迷藏没什么意思,还是不够寂寞,于是,他们决定玩寂寞无比的螃蟹版捉迷藏,顾名思义,就是说他们在玩游戏的时候只能沿水平或垂
直方向走。一番寂寞的剪刀石头布后,他们决定iPig去捉giPi。由于他们都很熟悉PKU的地形了,所以giPi只会躲在PKU内n个隐秘地点,显然
iPig也只会在那n个地点内找giPi。游戏一开始,他们选定一个地点,iPig保持不动,然后giPi用30秒的时间逃离现场(显然,giPi不会呆
在原地)。然后iPig会随机地去找giPi,直到找到为止。由于iPig很懒,所以他到总是走最短的路径,而且,他选择起始点不是随便选的,他想找一个
地点,使得该地点到最远的地点和最近的地点的距离差最小。iPig现在想知道这个距离差最小是多少。
由于iPig现在手上没有电脑,所以不能编程解决这个如此简单的问题,所以他马上打了个电话,要求你帮他解决这个问题。iPig告诉了你PKU的n个隐秘
地点的坐标,请你编程求出iPig的问题。

Input

第一行输入一个整数N 第2~N+1行,每行两个整数X,Y,表示第i个地点的坐标

Output

一个整数,为距离差的最小值。

Sample Input

4
0 0
1 0
0 1
1 1

Sample Output

1

HINT

对于30%的数据,N<=1000 对于100%的数据,N<=500000,0<=X,Y<=10^8 保证数据没有重点保证N>=2

Source

[Submit][Status][Discuss]

【BZOJ】【1941】【SDOI2010】Hide and Seek的更多相关文章

  1. 【BZOJ 2754 喵星球上的点名】

    Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 2512  Solved: 1092[Submit][Status][Discuss] Descript ...

  2. 【BZOJ】3052: [wc2013]糖果公园

    http://www.lydsy.com/JudgeOnline/problem.php?id=3052 题意:n个带颜色的点(m种),q次询问,每次询问x到y的路径上sum{w[次数]*v[颜色]} ...

  3. 【BZOJ】3319: 黑白树

    http://www.lydsy.com/JudgeOnline/problem.php?id=3319 题意:给一棵n节点的树(n<=1e6),m个操作(m<=1e6),每次操作有两种: ...

  4. 【BZOJ】3319: 黑白树(并查集+特殊的技巧/-树链剖分+线段树)

    http://www.lydsy.com/JudgeOnline/problem.php?id=3319 以为是模板题就复习了下hld............................. 然后n ...

  5. 【BZOJ】【1025】【SCOI2009】游戏

    DP/整数拆分 整个映射关系可以分解成几个循环(置换群的预备知识?),那么总行数就等于各个循环长度的最小公倍数+1(因为有个第一行的1~N).那么有多少种可能的排数就等于问有多少种可能的最小公倍数. ...

  6. 【BZOJ】1013: [JSOI2008]球形空间产生器sphere

    [BZOJ]1013: [JSOI2008]球形空间产生器sphere 题意:给n+1个n维的点的坐标,要你求出一个到这n+1个点距离相等的点的坐标: 思路:高斯消元即第i个点和第i+1个点处理出一个 ...

  7. 【BZOJ】1002:轮状病毒(基尔霍夫矩阵【附公式推导】或打表)

    Description 轮状病毒有很多变种,所有轮状病毒的变种都是从一个轮状基产生的.一个N轮状基由圆环上N个不同的基原子和圆心处一个核原子构成的,2个原子之间的边表示这2个原子之间的信息通道.如下图 ...

  8. 【BZOJ】【3697】采药人的路径&【3127】【USACO2013 Open】Yin and Yang

    点分治 Orz hzwer 倒是比较好想到点分治……然而在方案统计这里,我犯了两个错误…… 1.我比较傻逼的想的是:通过儿子来更新父亲,也就是统计以x为根的子树中xxxx的路径有多少条……这样转移. ...

  9. 【BZOJ】【3083】遥远的国度

    树链剖分/dfs序 其实过了[BZOJ][4034][HAOI2015]T2以后就好搞了…… 链修改+子树查询+换根 其实静态树的换根直接树链剖分就可以搞了…… 因为其实只有一样变了:子树 如果roo ...

随机推荐

  1. LoadRunner脚本篇

    LoadRunner脚本篇     1概述 2脚本录制 3脚本编写 4脚本调试   关  键  词:LoadRunner 性能测试脚本 摘      要:编写一个准确无误的脚本对性能测试有至关重要的意 ...

  2. 8-5 Wine trading in Gergovia Gergovia的酒交易 uva11054

    等价转换思维题 题意: 直线上有n(2<=n<=100000)个等距的村庄  每个村庄要么买酒 要么卖酒  设第i个村庄对酒的需求量为ai  绝对值小于一千  其中ai大于0表示买酒   ...

  3. 同一个IP不同端口号使用session失效

    背景 我有两个工程projectA.projectB,projectA放在TomcatA中,projectB放在TomcatB中,TomcatA.TomcatB在一台server上. 工程都映射的根路 ...

  4. Java 中类的初始化过程

    先来一张 JVM 中的内存模型 . 在Java 虚拟机原理这本书中介绍了类会被初始化的 5 种情况 . 1 遇到 new getstatic putstatic 和 invokestatic 这 4 ...

  5. odoo基础数据加载

    odoo 基础数据加载 这里介绍的odoo基础数据加载分两种方式,一种是演示数据加载,一种是默认数据加载,下面就是详细介绍 首先,当然是创建一个date文件夹 项目目录,右键自定义一个文件夹 XML数 ...

  6. 容易错的try

    配对方式一:try{}catch{}配对方法二:try{}catch{}finally{}配对方法三:try{}finally{}所以可用看出 catch和finally都不是必须的,try和catc ...

  7. 装部署VMware vSphere 5.5文档 (6-2) 为IBM x3850 X5服务器安装配置VMware ESXi

    部署VMware vSphere 5.5 实施文档 ########################################################################## ...

  8. python opencv3 cornerHarris 角点检测

    git:https://github.com/linyi0604/Computer-Vision 角点也是处在一个无论框框往哪边移动 框框内像素值都会变化很大的情况而定下来的点 如果框框水平方向上移动 ...

  9. Spring的优点

    Spring的优点 1.低侵入式设计,代码污染极低: 2.独立于各种应用服务器,基于Spring框架的应用,可以真正实现Write Once,Run Anywhere的承诺: 3.Spring的DI机 ...

  10. 「HNOI2018」转盘

    「HNOI2018」转盘 现场推出了大部分结论但是只写了 \(40\) 分暴力,被贺指导踩爆,现在还有点怀念 HNOI2018 贺指导对着镜子荒野行动的日子,那几天他云球迷瞎**指点篮球,被送上指导称 ...