2019.01.14 bzoj2648: SJY摆棋子(kd-tree)
传送门
kd−treekd-treekd−tree模板题。
题意简述:支持在平面上插入一个点,求对于一个点的最近点对。
思路:cdqcdqcdq是一种很不错的分治方法 只是好像码量有点窒息
所以我用了kd−treekd-treekd−tree来做。
其实就是两个维度分别作为键值来建立二叉搜索树即可。
这道题在查询时可以通过判断下一层的有没有可能对答案有贡献来剪个枝。
由于蒟蒻暂时没写过方差划分因此时间很慢。
代码:
#include<bits/stdc++.h>
#define mid (l+r>>1)
#define ri register int
using namespace std;
inline int read(){
int ans=0;
char ch=getchar();
while(!isdigit(ch))ch=getchar();
while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
return ans;
}
typedef pair<int,int> pii;
const int N=1e6+5,inf=0x3f3f3f3f;
int n,m,mx[N][2],mn[N][2],son[N][2],ans,tot=0;
bool cmpdir;
struct Pot{
int a[2];
inline int&operator[](const int&k){return a[k];}
inline const int&operator[](const int&k)const{return a[k];}
}p[N];
inline bool cmp(const Pot&a,const Pot&b){return a[cmpdir]<b[cmpdir];}
inline void pushup(int p){
if(son[p][0]){
mn[p][0]=min(mn[p][0],mn[son[p][0]][0]);
mx[p][0]=max(mx[p][0],mx[son[p][0]][0]);
mn[p][1]=min(mn[p][1],mn[son[p][0]][1]);
mx[p][1]=max(mx[p][1],mx[son[p][0]][1]);
}
if(son[p][1]){
mn[p][0]=min(mn[p][0],mn[son[p][1]][0]);
mx[p][0]=max(mx[p][0],mx[son[p][1]][0]);
mn[p][1]=min(mn[p][1],mn[son[p][1]][1]);
mx[p][1]=max(mx[p][1],mx[son[p][1]][1]);
}
}
inline int build(int l,int r,bool d){
cmpdir=d;
nth_element(p+l,p+mid,p+r+1,cmp);
mn[mid][0]=mx[mid][0]=p[mid][0],mn[mid][1]=mx[mid][1]=p[mid][1];
if(l<mid)son[mid][0]=build(l,mid-1,d^1);
if(r>mid)son[mid][1]=build(mid+1,r,d^1);
return pushup(mid),mid;
}
inline void insert(int k,Pot v,bool d){
if(p[k][d]<=v[d]){
if(son[k][1])insert(son[k][1],v,d^1);
else{
p[son[k][1]=++tot]=v;
mx[tot][0]=mn[tot][0]=v[0];
mx[tot][1]=mn[tot][1]=v[1];
}
}
else{
if(son[k][0])insert(son[k][0],v,d^1);
else{
p[son[k][0]=++tot]=v;
mx[tot][0]=mn[tot][0]=v[0];
mx[tot][1]=mn[tot][1]=v[1];
}
}
pushup(k);
}
inline int ask(int k,Pot v){
int ret=0;
ret+=max(0,mn[k][0]-v[0]);
ret+=max(0,mn[k][1]-v[1]);
ret+=max(0,v[0]-mx[k][0]);
ret+=max(0,v[1]-mx[k][1]);
return ret;
}
inline int dist(Pot a,Pot b){return abs(a[0]-b[0])+abs(a[1]-b[1]);}
inline void query(int k,Pot v,bool d){
ans=min(ans,dist(p[k],v));
int al=inf,ar=inf;
if(son[k][0])al=ask(son[k][0],v);
if(son[k][1])ar=ask(son[k][1],v);
if(al<ar){
if(al<ans)query(son[k][0],v,d^1);
if(ar<ans)query(son[k][1],v,d^1);
}
else{
if(ar<ans)query(son[k][1],v,d^1);
if(al<ans)query(son[k][0],v,d^1);
}
}
int main(){
n=read(),m=read(),tot=0;
for(ri i=1,x,y;i<=n;++i)x=read(),y=read(),p[++tot]=(Pot){x,y};
for(ri rt=build(1,n,0),i=1,x,y,op;i<=m;++i){
op=read(),x=read(),y=read();
if(op==1)insert(rt,(Pot){x,y},0);
else ans=inf,query(rt,(Pot){x,y},0),cout<<ans<<'\n';
}
return 0;
}
2019.01.14 bzoj2648: SJY摆棋子(kd-tree)的更多相关文章
- luogu4169 [Violet]天使玩偶/SJY摆棋子 / bzoj2648 SJY摆棋子 k-d tree
k-d tree + 重构的思想,就能卡过luogu和bzoj啦orz #include <algorithm> #include <iostream> #include &l ...
- BZOJ 2648: SJY摆棋子(K-D Tree)
Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 6051 Solved: 2113[Submit][Status][Discuss] Descript ...
- BZOJ2648: SJY摆棋子&&2716: [Violet 3]天使玩偶
BZOJ2648: SJY摆棋子 BZOJ2716: [Violet 3]天使玩偶 BZOJ氪金无极限... 其实这两道是同一题. 附上2648的题面: Description 这天,SJY显得无聊. ...
- [BZOJ2648] SJY摆棋子 kd-tree
2648: SJY摆棋子 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 5421 Solved: 1910[Submit][Status][Disc ...
- Bzoj2648 SJY摆棋子
Time Limit: 20 Sec Memory Limit: 128 MB Submit: 3128 Solved: 1067 Description 这天,SJY显得无聊.在家自己玩.在一个 ...
- bzoj 2648 SJY摆棋子 kd树
题目链接 初始的时候有一些棋子, 然后给两种操作, 一种是往上面放棋子. 另一种是给出一个棋子的位置, 问你离它最近的棋子的曼哈顿距离是多少. 写了指针版本的kd树, 感觉这个版本很好理解. #inc ...
- 【kd-tree】bzoj2648 SJY摆棋子
#include<cstdio> #include<cmath> #include<algorithm> using namespace std; #define ...
- BZOJ2648 SJY摆棋子(KD-Tree)
板子题. #include<iostream> #include<cstdio> #include<cmath> #include<cstdlib> # ...
- [bzoj2648]SJY摆棋子(带插入kd-tree)
解题关键:带插入kdtree模板题. #include<iostream> #include<cstdio> #include<cstring> #include& ...
随机推荐
- [leetcode]101. Symmetric Tree对称树
Given a binary tree, check whether it is a mirror of itself (ie, symmetric around its center). For e ...
- ThreadLocal ——android消息机制handler在非主线程创建not called Looper.prepare() 错误的原因
引用自:https://www.jianshu.com/p/a8fa72e708d3 引出: 使用Handler的时候,其必须要跟一个Looper绑定.在UI线程可直接初始化Handler来使用.但是 ...
- android的图片的初步学习理解
Android支持JPEG和PNG格式.GIF和BMP格式图片的支持. 图片最终要显示在屏幕上,都会对应一个屏幕上的点,即对应一个颜色值.不同格式的图片,只是不同压缩编码和解压算法. 也就是说,我们看 ...
- android抽屉效果
所谓抽屉 是区别于侧滑菜单 他不会把内容区域挤掉 他只是覆盖在内容区域 下边一个布局文件 一个代码 可以说的就是布局文件就是 <android.support.v4.widget.Dr ...
- opencv版本的问题
opencv版本的问题(由于我安装了两个版本opencv-2.4.9 and opencv-3.1.0),搜了相关问题,可以通过修改了CMakeList.txt解决此问题. 参考了这个博客:http: ...
- 以太坊难度炸弹是什么?极大抑制矿工继续以POW方式挖矿!
以太坊的“难度炸弹”(“Difficulty Bomb”)指的是,在挖掘算法中,使用以太币在区块链上对矿工进行奖励的难度越来越大.随着游戏变得更加复杂(矿工发现以太币难挣得多),在以太坊区块链上块的生 ...
- python 之 基础
变量 变量的作用: 标识符的命名规范: 掌握常量与变n量的区别: 变量定义规范: 声明变量: name='Alex Li' 三部分:变量名 赋值运算符 变量值 变量定义规则:1.变量名只能是字母.数字 ...
- 转:css知多少(1)——我来问你来答
1. 引言 各位前端或者伪前端(比如作者本人)的同志们,css对你们来说不是很陌生.比如我,在几年之前上大学的时候,给外面做网站就用css,而且必须用css.这样算下来也得六年多了,有些功能可能轻车熟 ...
- Android.HowToDefineCustomView
Custom View Errors E1 在使用自定义CustomView时,出现以下runtime error: Android.View.InflateException: Binary XML ...
- Codeforces 799D. String Game 二分
D. String Game time limit per test:2 seconds memory limit per test:512 megabytes input:standard inpu ...