Acwing 254. 天使玩偶

Ayu在七年前曾经收到过一个天使玩偶，当时她把它当做时间囊埋在了地下。

而七年后的今天，Ayu却忘了她把天使玩偶埋在了哪里，所以她决定仅凭一点模糊的记忆来寻找它。

我们把Ayu生活的小镇看做一个二维平面直角坐标系，而Ayu会不定时的记起可能在某个点(x,y)埋下了天使玩偶。

或者Ayu会询问你，假如她在(x,y)，那么她离最近的天使玩偶可能埋下的地方有多远。

因为Ayu只会沿着平行坐标轴的方向来行动，所以在这个问题里我们定义两个点之间的距离为曼哈顿距离：

dist(A,B)=|Ax−Bx|+|Ay−By|

其中Ax

表示点A的横坐标，其余类似。
输入格式

第一行包含两个整数n和m，在刚开始时，Ayu已经知道有n个点可能埋着天使玩偶，接下来Ayu要进行m次操作。

接下来n行，每行两个非负整数xi,yi

，表示初始n个点的坐标。

再接下来m行，每行三个非负整数 t,x,y 。

如果t=1，表示Ayu又回忆起了一个可能埋着玩偶的点(x,y)。

如果t=2，表示Ayu询问如果她在坐标(x,y)，那么在已经回忆出的点里，离她最近的那个点有多远。
输出格式

对于每个t=2的询问，在单独的一行内输出该询问的结果。
数据范围

n,m≤5∗105
,坐标范围为 0~106

。
输入样例：

2 3
1 1
2 3
2 1 2
1 3 3
2 4 2

输出样例：

1
2

基于时间的分治算法。
注意排序时三个关键字（横坐标、纵坐标、操作先后）都必须判断。
PS：同一点注意修改在查询前。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
const int max_n=500005,max_x=1e6+6,p=1e6,mim=-2139062144;
int c[max_x],n,m,ans[max_x];
const int dx[4]={1,1,-1,-1},dy[4]={1,-1,-1,1};
struct node
{
	int x,y,k,i;
}nd[max_n<<1];
struct Node
{ 
    int x,y,k,i;
}pd[max_n<<1];
bool cmpu(Node a,Node b){
    if(a.x==b.x&&a.y==b.y)return a.k<b.k;
    if(a.x==b.x)return a.y<b.y;
	return a.x<b.x;
}
bool cmpd(Node a,Node b){
    if(a.x==b.x&&a.y==b.y)return a.k<b.k;
	if(a.x==b.x)return a.y>b.y;
	return a.x>b.x;
}
int lowbit(int x){return x&(-x);}
void mp(int x,int y)
{
	pd[x].x=nd[y].x;
	pd[x].y=nd[y].y;
    pd[x].k=nd[y].k;
    pd[x].i=nd[y].i;
}
int get(int x)
{
	int cnt=mim;
	for(;x;x-=lowbit(x))cnt=max(cnt,c[x]);
	return cnt;
}
void add(int x,int v)
{
	for(;x<=1e6;x+=lowbit(x))c[x]=max(c[x],v);
}
void remo(int x,int v)
{
	for(;x<=1e6;x+=lowbit(x))c[x]=mim;
}
void solve(int l,int r)
{
	if(l==r)return ;
	int mid=(l+r)>>1;
	int tot=0;
	for(int i=l;i<=mid;i++){
		if(nd[i].k==1){
			mp(++tot,i);
		}
	}
	for(int i=mid+1;i<=r;i++){
		if(nd[i].k==2){
			mp(++tot,i);
		}
	}
	for(int i=0;i<4;i++)
	{
		if(i==0)sort(pd+1,pd+1+tot,cmpu);
		if(i==2)sort(pd+1,pd+1+tot,cmpd);
		for(int j=1;j<=tot;j++)
		{
			if(pd[j].k==1){
				add((p+pd[j].y*dy[i])%p,pd[j].x*dx[i]+pd[j].y*dy[i]);
			}else {
				ans[pd[j].i]=min(ans[pd[j].i],pd[j].x*dx[i]+pd[j].y*dy[i]-get((p+pd[j].y*dy[i])%p));
			}
		}
		for(int j=1;j<=tot;j++)
		{
			if(pd[j].k==1)
			remo((p+pd[j].y*dy[i])%p,pd[j].x*dx[i]+pd[j].y*dy[i]);
		}
    }
    
    solve(l,mid);solve(mid+1,r);
}
int main()
{
	cin>>n>>m;
	for(int i=1; i<=n;i++){
		cin>>nd[i].x>>nd[i].y;
		nd[i].k=1;
	}
	for(int i=1;i<=m;i++){
		cin>>nd[i+n].k;nd[i+n].i=i;
		cin>>nd[i+n].x>>nd[i+n].y;
	}
	memset(c,-0x80,sizeof(c));
	memset(ans,0x3f,sizeof(ans));
	solve(1,n+m);
	for(int i=1;i<=m;i++)
	if(nd[i+n].k==2)cout<<ans[i]<<endl;
	return 0;
}
/*
2 3
1 1
2 3
2 1 2
1 3 3
2 4 2
*/