【BZOJ 2716】天使玩偶

最新推荐文章于 2019-07-04 10:58:00 发布

forever_dreams

最新推荐文章于 2019-07-04 10:58:00 发布

阅读量3.9k

点赞数

分类专栏： --------数据结构-------- # CDQ分治文章标签： CDQ分治

本文链接：https://blog.csdn.net/forever_dreams/article/details/88593861

版权

--------数据结构-------- 同时被 2 个专栏收录

149 篇文章 1 订阅

订阅专栏

CDQ分治

5 篇文章 0 订阅

订阅专栏

【题目】

洛谷传送门

题目描述：
在这里插入图片描述
输入格式：

输出格式：

样例数据：

输入
2 3
1 1
2 3
2 1 2
1 3 3
2 4 2

输出
1
2

备注：
在这里插入图片描述

【分析】

我们把四个象限分开来考虑。

对于一个点 $(x, y)$ ，也就是求 $m i n (∣ x^{'} - x ∣ + ∣ y^{'} - y ∣)$ 。

假设 $(x^{'}, y^{'})$ 在 $(x, y)$ 的左下角。

那我们要求的就是 $m i n (x - x^{'} + y - y^{'}) = m i n (x + y - (x^{'} + y^{'}))$ 。

相当于求最大的 $x^{'} + y^{'}$ ，这个可以用 $\mathrm{CDQ}$ 分治搞定。

对于 $(x^{'}, y^{'})$ 在 $(x, y)$ 的其它位置的情况，直接翻转象限就行。

然后就可以 $\mathrm A$ 掉这道题啦。

【代码】

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1000005
#define lowbit(x) (x&-x)
using namespace std;
int n,m,Max,cnt,bit[N],ans[N];
struct node{int x,y,op,id,pos;}a[N],A[N];
bool comp(const node &p,const node &q)  {return p.x<q.x;}
void modify(int i,int x)  {for(;i<=Max;i+=lowbit(i))bit[i]=max(bit[i],x);}
int query(int i,int ans=0)  {for(;i;i-=lowbit(i))ans=max(ans,bit[i]);return ans;}
void Clear(int i)  {for(;i<=Max;i+=lowbit(i))bit[i]=0;}
void solve(int l,int r)
{
	if(l==r)  return;
	int mid=(l+r)>>1;
	int i,j,s=l,t=mid+1;
	for(i=l;i<=r;++i)
	{
		if(a[i].id<=mid)  A[s++]=a[i];
		else  A[t++]=a[i];
	}
	for(i=l;i<=r;++i)  a[i]=A[i];
	i=l;
	for(j=mid+1;j<=r;++j)
	{
		if(a[j].op!=2)  continue;
		for(;i<=mid&&a[i].x<=a[j].x;++i)
		  if(a[i].op==1)  modify(a[i].y,a[i].x+a[i].y);
		int now=query(a[j].y);
		if(now)  ans[a[j].pos]=min(ans[a[j].pos],a[j].x+a[j].y-now);
	}
	for(j=l;j<i;++j)  if(a[j].op==1)  Clear(a[j].y);
	solve(l,mid),solve(mid+1,r);
}
int main()
{
	int i;
	scanf("%d%d",&n,&m);
	memset(ans,127,sizeof(ans));
	for(i=1;i<=n;++i)
	{
		a[i].id=i,a[i].op=1;
		scanf("%d%d",&a[i].x,&a[i].y);
		a[i].x++,a[i].y++;
		Max=max(Max,max(a[i].x,a[i].y));
	}
	for(i=1;i<=m;++i)
	{
		a[n+i].id=n+i;
		scanf("%d%d%d",&a[n+i].op,&a[n+i].x,&a[n+i].y);
		a[n+i].x++,a[n+i].y++;
		if(a[n+i].op==2)  a[n+i].pos=++cnt;
		Max=max(Max,max(a[n+i].x,a[n+i].y));
	}
	Max++;
	sort(a+1,a+n+m+1,comp),solve(1,n+m);
	for(i=1;i<=n+m;++i)  a[i].x=-a[i].x+Max;
	sort(a+1,a+n+m+1,comp),solve(1,n+m);
	for(i=1;i<=n+m;++i)  a[i].y=-a[i].y+Max;
	sort(a+1,a+n+m+1,comp),solve(1,n+m);
	for(i=1;i<=n+m;++i)  a[i].x=-a[i].x+Max;
	sort(a+1,a+n+m+1,comp),solve(1,n+m);
	for(i=1;i<=cnt;++i)  printf("%d\n",ans[i]);
	return 0;
}