[bzoj2648]SJY摆棋子(带插入kd-tree)
解题关键:带插入kdtree模板题。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1000005
#define inf (1<<30)
using namespace std;
int n,m,dim,rt,ans;
struct node{int p[],x[],y[];}a[N];
bool cmp(node x,node y){ return x.p[dim]<y.p[dim]; }
struct kd_tree{
int c[N][];
node s[N],q;
void update(int k){//管辖范围
int l=c[k][],r=c[k][];
for(int i=;i<;i++){
if(l){ s[k].x[i]=min(s[k].x[i],s[l].x[i]); s[k].y[i]=max(s[k].y[i],s[l].y[i]); }
if(r){ s[k].x[i]=min(s[k].x[i],s[r].x[i]); s[k].y[i]=max(s[k].y[i],s[r].y[i]); }
}
}
void add(int k,node t){ for(int i=;i<;i++)s[k].x[i]=s[k].y[i]=s[k].p[i]=t.p[i]; }
int dist(node t,int k){
int tmp=;
for(int i=;i<;i++) tmp+=max(,s[k].x[i]-t.p[i]);
for(int i=;i<;i++) tmp+=max(,t.p[i]-s[k].y[i]);
return tmp;
}//?
void build(int &k,int l,int r,int now){
k=(l+r)>>; dim=now;
nth_element(a+l,a+k,a+r+,cmp);
add(k,a[k]);
if(l<k) build(c[k][],l,k-,now^);
if(k<r) build(c[k][],k+,r,now^);
update(k);
} void ins(int k,int now){
if(q.p[now]<s[k].p[now]){
if(c[k][]) ins(c[k][],now^);
else c[k][]=++n,add(n,q);
}
else{
if(c[k][]) ins(c[k][],now^);
else c[k][]=++n,add(n,q);
}
update(k);
}
void qry(int k){//曼哈顿距离,且只求最短,dis是最短距离
int tmp=;
for(int i=;i<;i++) tmp+=abs(s[k].p[i]-q.p[i]);
ans=min(ans,tmp);
int dl=c[k][]?dist(q,c[k][]):inf,dr=c[k][]?dist(q,c[k][]):inf;
if(dl<dr){
if(dl<ans) qry(c[k][]);
if(dr<ans) qry(c[k][]);
}else{
if(dr<ans) qry(c[k][]);
if(dl<ans) qry(c[k][]);
}
}
}kd; int main(){
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++) scanf("%d%d",&a[i].p[],&a[i].p[]);
kd.build(rt,,n,);
while(m--){
int k;
scanf("%d%d%d",&k,&kd.q.p[],&kd.q.p[]);
if(k==) kd.ins(rt,);
else{
ans=inf; kd.qry(rt); printf("%d\n",ans);
}
}
return ;
}
2、将2维普遍化。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1000005
#define inf (1<<30)
using namespace std;
int n,m,dim,rt,ans,k;
struct node{int p[],minn[],maxx[];}a[N];
bool cmp(node x,node y){ return x.p[dim]<y.p[dim]; }
struct kd_tree{
int c[N][];
node s[N],q;
void update(int o){//管辖范围
int l=c[o][],r=c[o][];
for(int i=;i<k;i++){
if(l){ s[o].minn[i]=min(s[o].minn[i],s[l].minn[i]); s[o].maxx[i]=max(s[o].maxx[i],s[l].maxx[i]); }
if(r){ s[o].minn[i]=min(s[o].minn[i],s[r].minn[i]); s[o].maxx[i]=max(s[o].maxx[i],s[r].maxx[i]); }
}
}
void add(int o,node t){ for(int i=;i<k;i++)s[o].minn[i]=s[o].maxx[i]=s[o].p[i]=t.p[i]; }
int dist(node t,int o){
int tmp=;
for(int i=;i<k;i++) tmp+=max(,s[o].minn[i]-t.p[i]);
for(int i=;i<k;i++) tmp+=max(,t.p[i]-s[o].maxx[i]);
return tmp;
}//?
void build(int &o,int l,int r,int now){
o=(l+r)>>; dim=now%k;
nth_element(a+l,a+o,a+r+,cmp);
add(o,a[o]);
if(l<o) build(c[o][],l,o-,now+);
if(o<r) build(c[o][],o+,r,now+);
update(o);
} void ins(int o,int now){
now%=k;
if(q.p[now]<s[o].p[now]){
if(c[o][]) ins(c[o][],now+);
else c[o][]=++n,add(n,q);
}
else{
if(c[o][]) ins(c[o][],now+);
else c[o][]=++n,add(n,q);
}
update(o);
}
void qry(int o){//曼哈顿距离,且只求最短,dis是最短距离
int tmp=;
for(int i=;i<k;i++) tmp+=abs(s[o].p[i]-q.p[i]);
ans=min(ans,tmp);
int dl=c[o][]?dist(q,c[o][]):inf,dr=c[o][]?dist(q,c[o][]):inf;
if(dl<dr){
if(dl<ans) qry(c[o][]);
if(dr<ans) qry(c[o][]);
}else{
if(dr<ans) qry(c[o][]);
if(dl<ans) qry(c[o][]);
}
}
}kd; int main(){
k=;
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++) scanf("%d%d",&a[i].p[],&a[i].p[]);
kd.build(rt,,n,);
while(m--){
int k;
scanf("%d%d%d",&k,&kd.q.p[],&kd.q.p[]);
if(k==) kd.ins(rt,);
else{
ans=inf; kd.qry(rt); printf("%d\n",ans);
}
}
return ;
}
[bzoj2648]SJY摆棋子(带插入kd-tree)的更多相关文章
- 【BZOJ-2648&2716】SJY摆棋子&天使玩偶 KD Tree
2648: SJY摆棋子 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 2459 Solved: 834[Submit][Status][Discu ...
- 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 ...
- luogu4169 [Violet]天使玩偶/SJY摆棋子 / bzoj2648 SJY摆棋子 k-d tree
k-d tree + 重构的思想,就能卡过luogu和bzoj啦orz #include <algorithm> #include <iostream> #include &l ...
- Bzoj2648 SJY摆棋子
Time Limit: 20 Sec Memory Limit: 128 MB Submit: 3128 Solved: 1067 Description 这天,SJY显得无聊.在家自己玩.在一个 ...
- 2019.01.14 bzoj2648: SJY摆棋子(kd-tree)
传送门 kd−treekd-treekd−tree模板题. 题意简述:支持在平面上插入一个点,求对于一个点的最近点对. 思路:cdqcdqcdq是一种很不错的分治方法 只是好像码量有点窒息 所以我用了 ...
- BZOJ2648 SJY摆棋子(KD-Tree)
板子题. #include<iostream> #include<cstdio> #include<cmath> #include<cstdlib> # ...
- 【kd-tree】bzoj2648 SJY摆棋子
#include<cstdio> #include<cmath> #include<algorithm> using namespace std; #define ...
- KDTree(Bzoj2648: SJY摆棋子)
题面 传送门 KDTree 大概就是一个分割\(k\)维空间的数据结构,二叉树 建立:每层选取一维为关键字,把中间的点拿出来,递归左右,有个\(STL\)函数nth_element可以用一下 维护:维 ...
随机推荐
- Linux环境下安装zookeeper
1. 下载安装文件zookeeper-3.4.6.tar.gz 镜像地址1: http://apache.fayea.com/zookeeper/ 镜像地址2: http://mirrors.hust ...
- Android 从上层到底层-----app层
CPU:RK3288 系统:Android 5.1 功能:上层 app 控制 led 亮灭 开发板:Firefly RK3288 MainActivity.java package com.aaron ...
- xunsearch之php索引维护(四)
1.添加文档 $xs = new XS('njw'); $index = $xs->index; $data = array( 'pid' => 234, // 此字段为主键,必须指定 ' ...
- 【Leetcode 371】Sum of Two Integers
问题描述:不使用+是或-操作符进行整数的加法运算 int getSum(int a, int b); 我的思路:把整数化成二进制进行运算,注意类型是int,也就是要考虑负数.关于负数的二进制表示可见之 ...
- Android UI一些技巧
(1)去掉EditText的背景颜色 android:background="@null" (2)ListView 修改某行的值,因为一些功能逻辑,需要修改ListView某行的 ...
- 第六章 hbase shell 命令
hbase shell命令 描述 alter 修改列族(Column Family)模式 count 统计表中行的数量 create 创建表 ...
- select函数源码阅读
fd_set结构体 #undef __NFDBITS /* It's easier to assume 8-bit bytes than to get CHAR_BIT. */ #define __N ...
- 使用wifi网卡笔记3---工具wpa_supplicant(STA模式)
1. wpa_supplicant介绍 supplicant是恳求者的意思,是wpa的发起者,是发送认证请求的设备(手机),手机--AP--认证服务器,可用于上述4种"认证/加密" ...
- Spring表单标签
虽然我们可以使用HTML原生的form表单标签来轻松的写出一个表单,其实我一直都是这样做的,但是使用Spring表单标签可以更方便我们完成例如:验证失败后表单数据的回填功能(虽然你可以使用EL+JST ...
- PoJ 1595 PrimeCuts
Prime Cuts Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 9339 Accepted: 3562 Descri ...