题目大意:给定一个二维平面,初始时有一些点了。然后有两种操作:
1.在(x,y)再加一个点
2.查询离(x,y)最近的点的距离。这里的距离是曼哈顿距离。
K-D tree,即K-Dimensional Tree,是一种高维索引树型数据结构。常用于大规模高维数据空间的最邻近或者K邻近查找,例如图像检索中高维图像特征向量的K邻近匹配,对KNN算法的优化等。
实现思路大体上就是:建出一棵BST.
怎么建呢?一维的我们都会建,多维的怎么办呢?简单来说,我们可以每次按不同的维度来划分。这样就把一个多维空间分成了若干小块,每一块内最多一个点。
calcmn函数:就是算出x点到p的所有子树的点的可能最小曼哈顿距离。用来判断p及其子树是否可能存在最优解。
至于复杂度分析,网上的讲解大都提到了要根据方差大小来均匀划分,保证复杂度,但是我还没有在oier的具体实现中看到这样的写法,均采用了比较简便的%k的方法来决定用哪一维来划分。不知这样的最坏复杂度是什么呢,望dalao指教。
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
inline int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x*f;
}
int n,m,D,ans,rt=0;
struct point{
int d[2];
point(int x=0,int y=0){d[0]=x;d[1]=y;}
int& operator[](int x){return d[x];}
friend bool operator<(point a,point b){return a[D]<b[D];}
friend int dis(point a,point b){return abs(a[0]-b[0])+abs(a[1]-b[1]);}
}P[N];
struct node{
point x;int mx[2],mn[2],lc,rc;
}tree[N<<1];
inline void update(int p){
int l=tree[p].lc,r=tree[p].rc;
for(int i=0;i<2;++i){
if(l){
tree[p].mn[i]=min(tree[p].mn[i],tree[l].mn[i]);
tree[p].mx[i]=max(tree[p].mx[i],tree[l].mx[i]);
}if(r){
tree[p].mn[i]=min(tree[p].mn[i],tree[r].mn[i]);
tree[p].mx[i]=max(tree[p].mx[i],tree[r].mx[i]);
}
}
}
inline void build(int &p,int l,int r,int op){
int mid=l+r>>1;p=mid;D=op;
nth_element(P+l,P+mid,P+r+1);tree[p].x=P[mid];
for(int i=0;i<2;++i) tree[p].mx[i]=tree[p].mn[i]=tree[p].x[i];
if(l<mid) build(tree[p].lc,l,mid-1,op^1);
if(r>mid) build(tree[p].rc,mid+1,r,op^1);
update(p);
}
inline void ins(int &p,point x,int op){
if(!p){
p=++n;tree[p].x=x;
for(int i=0;i<2;++i) tree[p].mn[i]=tree[p].mx[i]=tree[p].x[i];return;
}if(x[op]<=tree[p].x[op]) ins(tree[p].lc,x,op^1);
else ins(tree[p].rc,x,op^1);update(p);
}
inline int calcmn(int p,point x){
if(!p) return inf;int res=0;
for(int i=0;i<2;++i){
if(x[i]<tree[p].mn[i]) res+=tree[p].mn[i]-x[i];
if(x[i]>tree[p].mx[i]) res+=x[i]-tree[p].mx[i];
}return res;
}
inline void query(int p,point x){
ans=min(ans,dis(tree[p].x,x));
int dl=calcmn(tree[p].lc,x),dr=calcmn(tree[p].rc,x);
if(dl<dr){
if(dl<ans) query(tree[p].lc,x);
if(dr<ans) query(tree[p].rc,x);
}else{
if(dr<ans) query(tree[p].rc,x);
if(dl<ans) query(tree[p].lc,x);
}
}
int main(){
// freopen("a.in","r",stdin);
n=read();m=read();
for(int i=1;i<=n;++i)
for(int k=0;k<2;++k) P[i][k]=read();
build(rt,1,n,0);
while(m--){
int op=read(),x=read(),y=read();
if(op==1) ins(rt,point(x,y),0);
else{
ans=inf;query(rt,point(x,y));
printf("%d\n",ans);
}
}return 0;
}