[BZOJ]2716: [Violet 3]天使玩偶 CDQ分治+树状数组

Description

题解：

这类问题好像可以用kd-tree？然而蒟蒻并不会。不过CDQ分治的做法也挺巧妙的。考虑一个点(x,y)，如果只考虑它左下角的离他最近的点，那么我们要找的就是所有满足x1<=x,y1<=y的点(x1,y1)中，x+y−(x1+y1)的最大值，那么这个东西显然是可以用CDQ分治求的，那么我们只要做四次，就可以求出答案了，不过速度很慢……

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=500010;
const int Max=1000001;
const int inf=-1044266559;
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<<3)+(x<<1)+(ch^48),ch=getchar();
return x*f;
}
int s[Max+10],mxx=-1,mxy=-1;
void modify(int x,int y){if(!x)return;for(;x<=Max;x+=(x&-x))s[x]=((y==inf)?y:max(s[x],y));}
int getmax(int x){int re=inf;for(;x;x-=(x&-x))re=max(re,s[x]);return re;}
struct Opt{int x,y,type,id;}q[Maxn<<1],temp[Maxn<<1],p[Maxn<<1];
int n,m,ans[Maxn],qcnt=0;
void solve(int l,int r)
{
if(l==r)return;
int mid=l+r>>1;
solve(l,mid);solve(mid+1,r);
int i=l,j=mid+1,len=0;
while(i<=mid&&j<=r)
{
if(q[i].x<=q[j].x)
{
if(q[i].type==1)modify(q[i].y,q[i].x+q[i].y);
temp[++len]=q[i++];
}
else
{
if(q[j].type==2)ans[q[j].id]=min(ans[q[j].id],q[j].x+q[j].y-getmax(q[j].y));
temp[++len]=q[j++];
}
}
int t=i;
while(i<=mid)temp[++len]=q[i++];
while(j<=r)
{
if(q[j].type==2)ans[q[j].id]=min(ans[q[j].id],q[j].x+q[j].y-getmax(q[j].y));
temp[++len]=q[j++];
}
for(int p=l;p<t;p++)if(q[p].type==1)modify(q[p].y,inf);
for(int p=1;p<=len;p++)q[p+l-1]=temp[p];
}
int main()
{
memset(s,-63,sizeof(s));
memset(ans,127,sizeof(ans));
n=read(),m=read();
for(int i=1;i<=n;i++)p[i].type=1,p[i].x=read()+1,p[i].y=read()+1,mxx=max(mxx,p[i].x),mxy=max(mxy,p[i].y);
for(int i=1;i<=m;i++)
{
p[i+n].type=read(),p[i+n].x=read()+1,p[i+n].y=read()+1,mxx=max(mxx,p[i+n].x),mxy=max(mxy,p[i+n].y);
if(p[i+n].type==2)p[i+n].id=++qcnt;
}
mxx++;mxy++;
for(int i=1;i<=n+m;i++)q[i].x=p[i].x,q[i].y=p[i].y,q[i].type=p[i].type,q[i].id=p[i].id;solve(1,n+m);
for(int i=1;i<=n+m;i++)q[i].x=p[i].y,q[i].y=mxx-p[i].x,q[i].type=p[i].type,q[i].id=p[i].id;solve(1,n+m);
for(int i=1;i<=n+m;i++)q[i].x=mxy-p[i].y,q[i].y=mxx-p[i].x,q[i].type=p[i].type,q[i].id=p[i].id;solve(1,n+m);
for(int i=1;i<=n+m;i++)q[i].x=mxy-p[i].y,q[i].y=p[i].x,q[i].type=p[i].type,q[i].id=p[i].id;solve(1,n+m);
for(int i=1;i<=qcnt;i++)printf("%d\n",ans[i]);
}