【VIOLET3】【BZOJ2716】天使玩偶（cdq分治）

传送门

把四个象限分开讨论（不然没法做对吧）

对于一个点$(x,y)$,也就是求min(∣x′−x∣+∣y′−y∣)min(|x'-x|+|y'-y|)min(∣x′−x∣+∣y′−y∣)
考虑左下角，其他的类似
则求min(x−x′+y−y′)=min(x+y−x′−y′)min(x-x'+y-y')=min(x+y-x'-y')min(x−x′+y−y′)=min(x+y−x′−y′)
那就相当于求最大的x′+y′x'+y'x′+y′
$cdq$就行了

注意坐标翻转会出现负数，加一个极大值

#include
using namespace std;
inline int read(){
	char ch=getchar();
	int res=0,f=1;
	while(!isdigit(ch)){if(ch=='-')f=-f;ch=getchar();}
	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=getchar();
	return res*f;
}
const int N=1000005;
struct ask{
	int x,y,pos,idx,op;
}p[N],q[N];
inline bool comp(const ask&a,const ask&b){
	return a.x<b.x||(a.x==b.x&&a.op<b.op);
}
int ans[N],n,m,mx,cnt,tot,tr[N];
inline int lowbit(int x){
	return x&(-x);
}
inline void update(int p,int k){
	for(;p<=mx;p+=lowbit(p))tr[p]=max(tr[p],k);
}
inline int query(int p,int res=0){
	for(;p;p-=lowbit(p))res=max(res,tr[p]);
	return res;
}
inline void clear(int p){
	for(;p<=mx;p+=lowbit(p))tr[p]=0;
}
void cdq(int l,int r){
	if(l==r)return;
	int mid=(l+r)>>1;
	int st=l,des=mid+1;
	for(int i=l;i<=r;i++){
		if(p[i].idx<=mid)q[st++]=p[i];
		else q[des++]=p[i];
	}
	for(int i=l;i<=r;i++)p[i]=q[i];
	int i=l;
	for(int j=mid+1;j<=r;j++){
		if(p[j].op!=2)continue;
		while(i<=mid&&p[i].x<=p[j].x){
			if(p[i].op==1)update(p[i].y,p[i].x+p[i].y);
			i++;
		}
		int t=query(p[j].y);
		if(t){
			t=p[j].x+p[j].y-t;
			ans[p[j].pos]=min(ans[p[j].pos],t);
		}
	}
	for(int j=l;j<=i;j++){
		if(p[j].op==1)clear(p[j].y);
	}
	cdq(l,mid),cdq(mid+1,r);
}
int main(){
	n=read(),m=read();
	memset(ans,127,sizeof(ans));
	for(int i=1;i<=n;i++){
		p[++tot].x=read()+1,p[tot].y=read()+1;
		p[tot].idx=tot,p[tot].op=1;
		mx=max(mx,p[tot].x),mx=max(mx,p[tot].y);
	}
	for(int i=1;i<=m;i++){
		p[++tot].op=read();
		p[tot].x=read()+1,p[tot].y=read()+1;
		p[tot].idx=tot;if(p[tot].op==2)p[tot].pos=++cnt;
		mx=max(mx,p[tot].x),mx=max(mx,p[tot].y);
	}mx++;
	sort(p+1,p+tot+1,comp);
	cdq(1,tot);
	for(int i=1;i<=tot;i++)p[i].x=-p[i].x+mx;
	sort(p+1,p+tot+1,comp);
	cdq(1,tot);
	for(int i=1;i<=tot;i++)p[i].y=-p[i].y+mx;
	sort(p+1,p+tot+1,comp);
	cdq(1,tot);
	for(int i=1;i<=tot;i++)p[i].x=-p[i].x+mx;
	sort(p+1,p+tot+1,comp);
	cdq(1,tot);
	for(int i=1;i<=cnt;i++)cout<<ans[i]<<'\n';
}

转载于:https://www.cnblogs.com/stargazer-cyk/p/11145623.html