bzoj2716/2648 / P4169 [Violet]天使玩偶/SJY摆棋子
k-d tree 模板
找了好几天才发现输出优化错了....真是zz......
当子树非常不平衡时,就用替罪羊树的思想,拍扁重建。
luogu有个点开了O2才过.......
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
using namespace std;
char c;void read(int &x){
c=getchar();x=;int f=;
while(!isdigit(c)) c=getchar(),f=(f&&c!='-');
while(isdigit(c)) x=(x<<)+(x<<)+(c^),c=getchar();
x=f?x:-x;
}
void print(int x){
if(x<) putchar('-'),x=-x;
if(x>) print(x/);
putchar(x%+);
}
inline int max(int a,int b){return a>b?a:b;}
inline int min(int a,int b){return a<b?a:b;}
inline int abs(int a){return a<?-a:a;}
#define N 1000010
#define alph (0.75)
#define inf 1e9+7
int n,m,rt,u,ans;
int lit[N],tp1,tp2;
int d;struct node{
int p[];
void init(){read(p[]);read(p[]);}
bool operator < (const node &t) const{
return p[d]<t.p[d];
}
}b[N];
struct data{int mx[],my[],ch[],sz,fa;node ad;}a[N];
#define lc a[o].ch[0]
#define rc a[o].ch[1]
inline void up(int o){
a[o].mx[]=max(a[o].mx[],max(a[lc].mx[],a[rc].mx[]));
a[o].mx[]=min(a[o].mx[],min(a[lc].mx[],a[rc].mx[]));
a[o].my[]=max(a[o].my[],max(a[lc].my[],a[rc].my[]));
a[o].my[]=min(a[o].my[],min(a[lc].my[],a[rc].my[]));
a[o].sz=a[lc].sz+a[rc].sz+;
}
inline void news(int &o,node w){
o=tp1?lit[tp1--]:++u; a[o].ad=w;
lc=rc=; a[o].sz=;
a[o].mx[]=a[o].mx[]=w.p[];
a[o].my[]=a[o].my[]=w.p[];
}
void build(int &o,int l,int r,int wd){
if(l>r) return ;
int mid=l+((r-l)>>);
d=wd;nth_element(b+l,b+mid,b+r+);
news(o,b[mid]);
build(lc,l,mid-,wd^),a[lc].fa=o;
build(rc,mid+,r,wd^),a[rc].fa=o;
up(o);
}
void pia(int o){
if(lc) pia(lc);
lit[++tp1]=o; b[++tp2]=a[o].ad;
if(rc) pia(rc);
}
inline void check(int &o,int wd){
if(alph*a[o].sz<max(a[lc].sz,a[rc].sz))
tp2=,pia(o),build(o,,tp2,wd);
}//是否不平衡需要重建
void ins(int &o,int wd,node w){
if(!o) {news(o,w); return ;}
d=wd;ins(a[o].ch[a[o].ad<w],wd^,w);
up(o);check(o,wd);
}
inline int gdis(int o,node w){
if(!o) return inf;
int res=;
if(w.p[]>a[o].mx[]) res+=w.p[]-a[o].mx[];
if(w.p[]<a[o].mx[]) res+=a[o].mx[]-w.p[];
if(w.p[]>a[o].my[]) res+=w.p[]-a[o].my[];
if(w.p[]<a[o].my[]) res+=a[o].my[]-w.p[];
return res;
}//点与矩形的最小距离
inline int dis(node A,node B){
return abs(A.p[]-B.p[])+abs(A.p[]-B.p[]);
}
void query(int o,node w){
if(!o) return ;
ans=min(ans,dis(a[o].ad,w));
int D[]={gdis(lc,w),gdis(rc,w)},r=D[]>D[];
if(D[r]<ans) query(a[o].ch[r],w);
if(D[r^]<ans) query(a[o].ch[r^],w);
}
int main(){
read(n);read(m);
a[].mx[]=a[].my[]=-inf;
a[].mx[]=a[].my[]=inf;
for(register int i=;i<=n;++i) b[i].init();
build(rt,,n,); node q1;
for(register int i=,q2;i<=m;++i){
read(q2); q1.init();
if(q2==) ins(rt,,q1);
else{
ans=inf,query(rt,q1);
print(ans); putchar('\n');
}
}return ;
}
bzoj2716/2648 / P4169 [Violet]天使玩偶/SJY摆棋子的更多相关文章
- 洛谷 P4169 [Violet]天使玩偶/SJY摆棋子 解题报告
P4169 [Violet]天使玩偶/SJY摆棋子 题目描述 \(Ayu\)在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下.而七年后 的今天,\(Ayu\) 却忘了她把天使玩偶埋在了哪 ...
- 洛谷P4169 [Violet]天使玩偶/SJY摆棋子(CDQ分治)
[Violet]天使玩偶/SJY摆棋子 题目传送门 解题思路 用CDQ分治开了氧气跑过. 将输入给的顺序作为第一维的时间,x为第二维,y为第三维.对于距离一个询问(ax,ay),将询问分为四块,左上, ...
- P4169 [Violet]天使玩偶/SJY摆棋子
题目背景 感谢@浮尘ii 提供的一组hack数据 题目描述 Ayu 在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下.而七年后 的今天,Ayu 却忘了她把天使玩偶埋在了哪里,所以她决定仅 ...
- 洛谷P4169 [Violet]天使玩偶/SJY摆棋子
%%%神仙\(SJY\) 题目大意: 一个二维平面,有两种操作: \(1.\)增加一个点\((x,y)\) \(2.\)询问距离\((x,y)\)曼哈顿最近的一个点有多远 \(n,m\le 300 0 ...
- Luogu P4169 [Violet]天使玩偶/SJY摆棋子
传送门 二维平面修改+查询,cdq分治可以解决. 求关于某个点曼哈顿距离(x,y坐标)最近的点——dis(A,B) = |Ax-Bx|+|Ay-By| 但是如何去掉绝对值呢? 查看题解发现假设所有的点 ...
- luoguP4169 [Violet]天使玩偶/SJY摆棋子 K-Dtree
P4169 [Violet]天使玩偶/SJY摆棋子 链接 luogu 思路 luogu以前用CDQ一直过不去. bzoj还是卡时过去的. 今天终于用k-dtree给过了. 代码 #include &l ...
- [Violet]天使玩偶/SJY摆棋子 [cdq分治]
P4169 [Violet]天使玩偶/SJY摆棋子 求离 \((x,y)\) 最近点的距离 距离的定义是 \(|x1-x2|+|y1-y2|\) 直接cdq 4次 考虑左上右上左下右下就可以了-略微卡 ...
- 【LG4169】[Violet]天使玩偶/SJY摆棋子
[LG4169][Violet]天使玩偶/SJY摆棋子 题面 洛谷 题解 至于\(cdq\)分治的解法,以前写过 \(kdTree\)的解法好像还\(sb\)一些 就是记一下子树的横.纵坐标最值然后求 ...
- LG4169 [Violet]天使玩偶/SJY摆棋子
题意 Ayu 在七年前曾经收到过一个天使玩偶,当时她把它当作时间囊埋在了地下.而七年后 的今天,Ayu 却忘了她把天使玩偶埋在了哪里,所以她决定仅凭一点模糊的记忆来寻找它. 我们把 Ayu 生活的小镇 ...
随机推荐
- python3 判断字符串是否为IP
#!/usr/bin/python3 # -*- coding: utf-8 -*- import re ip = "192.168.1.1" ip = re.findall(&q ...
- golang编程之我见
golang编程之我见 学习了两个月的golang,语法算是基本掌握了,从一个C++程序员的角度,提出自己的几个看法吧. 1,没有一个好的包管理工具. 我在公司用的是glide的包管理,这个工具的好处 ...
- iOS 新浪微博-5.0 首页微博列表
首页显示微博列表,是微博的核心部分,这一章节,我们主要是显示出微博的列表. 导入第三方类库 pod 'SDWebImage', '~> 3.7.3' pod 'MJRefresh', '~> ...
- Windows服务器时间不同步问题
一台域内的服务器时间不停地被修改,我先向用户收集了一些信息 只有这一台出现此问题,其他服务器均为正常(补充一下,问题快解决完的时候用户告诉我一个重要的消息,就是时间被修改了一段时间后自动会被修改回去) ...
- JavaScript-模拟收银台小程序
<!DOCTYPE html> <html> <head lang="en"> <meta charset="UTF-8&quo ...
- django 网站的搭建(2)
这里使用nginx+uwsgi的方法来搭建生产环境 1,pip3.5 install uwsgi 下载uwsgi ,这里就不做测试了,一般不会出错 2,将django与uwsgi连接在一起 毕竟ru ...
- 关于git上的一些错误信息
如果输入$ Git remote add origin git@github.com:djqiang(github帐号名)/gitdemo(项目名).git 提示出错信息:fatal: remote ...
- mac下编译cpu only caffe并用xCode建caffe工程
mac编译caffe 好像又变容易了,直接git clone下载blvc源码,make.config里去掉了CPU_ONLY前面的注释,并没有安装任何依赖,也可能是自己mac上本来有, xCode里调 ...
- 做一次面向对象的体操:将JSON字符串转换为嵌套对象的一种方法
背景与问题 在 <一个略复杂的数据映射聚合例子及代码重构> 一文中,将一个JSON字符串转成了所需要的订单信息Map.尽管做了代码重构和配置化,过程式的代码仍然显得晦涩难懂,并且客户端使用 ...
- flask 在模板中渲染表单
在模板中渲染表单 为了能够在模板中渲染表单,我们需要把表单类实例传入模板.首先在视图函数里实例化表单类LoginForm,然后再render_template()函数中使用关键脑子参数form将表单实 ...