http://www.elijahqi.win/2018/01/13/bzoj4533-beijing2014-wintercamp-%e6%95%b0%e6%8d%ae/
Description
为了写论文,Alex经常要整理大量的数据。 这一次,Alex面临一个严峻的考验:他需要实现一个数据结构来维护一个点集。
现在,二维平面上有N个点。Alex 需要实现以下三种操作:
1. 在点集里添加一个点;
2. 给出一个点,查询它到点集里所有点的曼哈顿距离的最小值;
3. 给出一个点,查询它到点集里所有点的曼哈顿距离的最大值。
两个点的曼哈顿距离定义为它们的横坐标差的绝对值与纵坐标差的绝对值的和。这么困难的问题,Alex当然不会做,只好再次请你帮忙了。
Input
第一行包含一个整数N,表示点集最初的点数。
接下来N行,每行两个整数,依次表示每个点的横坐标和纵坐标。
第N+2行包含一个整数Q,表示询问的数目。
接下来Q行,每行三个整数,依次表示询问的类型,点的横坐标和纵坐标。0类型表示添加一个点,1类型表示查询到该点的曼哈顿距离的最小值,2类型表示查询最大值。
1 ≤ N, Q ≤ 100,000,点的坐标是不超过10^9的非负整数
Output
输出若干行,依次表示每个查询操作的答案。
Sample Input
3
7 5
6 2
3 1
5
1 6 1
1 5 5
2 7 1
0 3 2
1 1 0
Sample Output
1
2
4
3
复习熟练一下kd-tree的写法qwq
求曼哈顿距离最小值:
估价函数为该点到边框矩形的最小值 然后优先去做估价小的子树 程序中的calc1函数
求曼哈顿距离最大值:
估价函数为该点到边框矩形的最大值 然后优先去做估价大的子树 程序中的calc2函数
为了保证是时间复杂度 每10000个暴力重构 注意在建树以及insert1的时候要将左右子树指针清空
#include<cstdio>
#include<algorithm>
#define N 220000
#define inf 0x7f7f7f7f
using namespace std;
inline char gc(){
static char now[1<<16],*S,*T;
if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0;char ch=gc();
while(ch<'0'||ch>'9') ch=gc();
while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
return x;
}
int root,D,n,ans,m;
struct node1{
int d[2];
int& operator[](int x){return d[x];}
friend bool operator<(node1 a,node1 b){return a[D]<b[D];}
}point[N],P;
struct node{
node1 x;int left,right,min[2],max[2];
}tree[N];
inline void update(int x){
int l=tree[x].left,r=tree[x].right;
for (int i=0;i<2;++i) tree[x].min[i]=min(tree[x].min[i],min(tree[l].min[i],tree[r].min[i]));
for (int i=0;i<2;++i) tree[x].max[i]=max(tree[x].max[i],max(tree[l].max[i],tree[r].max[i]));
}
inline void build(int &x,int l,int r,int dim){
D=dim;int mid=l+r>>1;x=mid;nth_element(point+l,point+mid,point+r+1);tree[x].left=tree[x].right=0;
tree[x].x=point[mid];for (int i=0;i<2;++i) tree[x].min[i]=tree[x].max[i]=point[mid][i];
if (l<mid) build(tree[x].left,l,mid-1,dim^1);
if (r>mid) build(tree[x].right,mid+1,r,dim^1);update(x);
}
inline void insert1(int &x,int dim){
if (!x) {x=n;tree[x].x=point[n];tree[x].left=tree[x].right=0;
for (int i=0;i<2;++i) tree[x].min[i]=tree[x].max[i]=point[n][i];return;}
if (point[n][dim]<=tree[x].x[dim]) insert1(tree[x].left,dim^1);else insert1(tree[x].right,dim^1);update(x);
}
inline int calc(node1 a){
return abs(a[0]-P[0])+abs(a[1]-P[1]);
}
inline int calc1(int x){
int tmp=0;
for (int i=0;i<2;++i) tmp+=max(0,tree[x].min[i]-P[i]),tmp+=max(0,P[i]-tree[x].max[i]);return tmp;
}
inline int calc2(int x){
int tmp=0;
for (int i=0;i<2;++i) tmp+=max(abs(tree[x].max[i]-P[i]),abs(P[i]-tree[x].min[i]));return tmp;
}
inline void query(int x){
ans=min(ans,calc(tree[x].x));int disl=inf,disr=inf;
if (tree[x].left) disl=calc1(tree[x].left);if(tree[x].right) disr=calc1(tree[x].right);
if (disl<disr){
if (disl<ans) query(tree[x].left);if (disr<ans) query(tree[x].right);
}else{
if (disr<ans) query(tree[x].right);if (disl<ans) query(tree[x].left);
}
}
inline void query1(int x){
ans=max(ans,calc(tree[x].x));int disl=0,disr=0;
if (tree[x].left) disl=calc2(tree[x].left);if (tree[x].right) disr=calc2(tree[x].right);
if(disl>disr){
if(disl>ans) query1(tree[x].left);if (disr>ans) query1(tree[x].right);
}else{
if (disr>ans) query1(tree[x].right);if (disl>ans) query1(tree[x].left);
}
}
int main(){
freopen("bzoj4533.in","r",stdin);
n=read();for (int i=0;i<2;++i) tree[0].min[i]=inf;
for (int i=1;i<=n;++i) for (int j=0;j<2;++j) point[i][j]=read();
build(root,1,n,0);m=read();int delta=0;
for (int i=1;i<=m;++i){
int op=read();
if (!op){
++delta;++n;for (int j=0;j<2;++j) point[n][j]=read();
if (delta%10000==0) build(root,1,n,0);else insert1(root,0);
}
if (op==1){
ans=inf;for (int j=0;j<2;++j) P[j]=read();
query(root);printf("%d\n",ans);
}
if (op==2){
ans=0;for (int j=0;j<2;++j) P[j]=read();
query1(root);printf("%d\n",ans);
}
}
return 0;
}