BZOJ3958[WF2011]Mummy Madness 二分+扫描线+线段树

Posted CQzhangyu

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ3958[WF2011]Mummy Madness 二分+扫描线+线段树相关的知识,希望对你有一定的参考价值。

【BZOJ3958】[WF2011]Mummy Madness

Description

在2011年ACM-ICPC World Finals上的一次游览中,你碰到了一个埃及古墓。
不幸的是,你打开了坟墓之后,才发现这是一个坏主意:突然之间,原本空无一物的沙漠上已经爬满了暴躁的木乃伊。(如果你也沉睡几千年而突然被惊醒,你也会变得如此暴躁的。)(幸运的是,当你做完这道题的时候,你醒来了,发现你在弗罗里达的酒店里。那些木乃伊只是一场梦。)
面对这一大堆疯狂的木乃伊,你唯一的机会就是试图在他们抓到你之前逃跑。问题是:假如你与木乃伊永不疲倦,那么经过多长时间你会被木乃伊抓到?
我们把沙漠看成一个正方形的网格,你与木乃伊轮流移动(你走出第一步)。轮到你时,你可以移动到相邻的8个格子之一,或者站着不动。轮到木乃伊时,每个木乃伊会移动到其相邻的格子之一,使得他与你的欧几里得距离尽量小(假设你与木乃伊都站在格子的中心位置)。允许多个木乃伊同时占据同一个格子。
在每个单位时间内,你先做出移动,然后木乃伊做出移动。如果你与任何一个木乃伊站在同一位置,你会被抓住。当然,你试图尽量长时间避免被抓住。经过多少单位时间你会被抓住呢?
下图描述了你被4个木乃伊追逐的例子。H代表你的初始位置,而M代表木乃伊的初始位置。以你的初始位置为原点,则经过4个单位时间后,你被初始位置为(3,4)的木乃伊抓住。
技术分享图片

Input

输入文件包含若干组数据。每组数据的第一行为一个数n(0≤n≤10^5),表示沙漠中木乃伊的个数。接下来n行,每行两个整数x y,表示初始时在(x,y)有一个木乃伊。x,y的绝对值均不超过10^6。你的初始位置是(0,0),保证一开始这里没有木乃伊。
输入文件以一行-1结束。

Output

对于每组测试数据,输出一行,包括它的编号和被抓住经过的最长时间(即你做出决策的次数);或输出"never",如果你有办法永远不被抓住。
请以样例输出的格式输出数据。

Sample Input

4
-3 5
3 4
-6 -2
1 -5
1
0 -1
-1

Sample Output

Case 1: 4
Case 2: never

HINT

对于100%的数据,n≤10^5

题解:容易想到二分答案mid,因为如果在mid时刻能抓到以后也一直是抓到的。在mid时刻,人可能到达的位置是一个正方形,每个木乃伊能到达的位置也是正方形,我们可以发现如果在mid时刻被抓了,当且仅当木乃伊的移动位置将人的移动位置完全覆盖了。所以我们只需要判断若干个矩形的并是否能覆盖给定矩形即可,用扫描线+线段树可以实现。

但是本人从来没写过矩形并啊,这里学习一发:

考虑标记永久化,我们对于线段树上的每个节点维护cnt:区间中被覆盖的位置个数,sum:当前区间被覆盖了几层。其中cnt相当于值,可以pushup,而sum相当于标记,但不能pushdown。如果一个区间的sum>0,我们令cnt=r-l+1;否则cnt=cnt[lson]+cnt[rson]。

 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define lson x<<1
#define rson x<<1|1
using namespace std;
const int maxn=100010;
int n,tot;
struct node
{
	int l,r,x,k;
	node() {}
	node(int a,int b,int c,int d) {l=a,r=b,x=c,k=d;}
}p[maxn<<1];
int x[maxn],y[maxn],s1[maxn*80],s[maxn*80];
bool tc[maxn*80];
bool cmp(const node &a,const node &b)
{
	return a.x<b.x;
}
void updata(int l,int r,int x,int a,int b,int c)
{
	if(l!=r&&tc[x])	tc[lson]=tc[rson]=1,s1[lson]=s1[rson]=s[lson]=s[rson]=tc[x]=0;
	if(a<=l&&r<=b)
	{
		s[x]+=c;
		if(s[x])	s1[x]=r-l+1;
		else	if(l!=r)	s1[x]=s1[lson]+s1[rson];
		else	s1[x]=0;
		return ;
	}
	int mid=(l+r)>>1;
	if(a<=mid)	updata(l,mid,lson,a,b,c);
	if(b>mid)	updata(mid+1,r,rson,a,b,c);
	if(s[x])	s1[x]=r-l+1;
	else	if(l!=r)	s1[x]=s1[lson]+s1[rson];
	else	s1[x]=0;
}
inline bool check(int mid)
{
	int i;
	tot=0;
	for(i=1;i<=n;i++)
	{
		if(x[i]+mid<-mid||x[i]-mid>mid||y[i]+mid<-mid||y[i]-mid>mid)	continue;
		p[++tot]=node(max(-mid,y[i]-mid),min(mid,y[i]+mid),max(-mid,x[i]-mid),1);
		p[++tot]=node(max(-mid,y[i]-mid),min(mid,y[i]+mid),min(mid+1,x[i]+mid+1),-1);
	}
	sort(p+1,p+tot+1,cmp);
	tc[1]=1,s[1]=s1[1]=0;
	p[0].x=-mid,p[tot+1]=node(0,0,mid+1,0);
	for(i=1;i<=tot+1;i++)
	{
		if(p[i].x>p[i-1].x&&s1[1]<2*mid+1)	return 1;
		updata(-mid,mid,1,p[i].l,p[i].r,p[i].k);
	}
	return 0;
}
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<‘0‘||gc>‘9‘)	{if(gc==‘-‘)	f=-f;	gc=getchar();}
	while(gc>=‘0‘&&gc<=‘9‘)	ret=ret*10+(gc^‘0‘),gc=getchar();
	return ret*f;
}
int main()
{
	for(int cas=1;;cas++)
	{
		n=rd();
		if(n==-1)	return 0;
		int i,l=0,r=1000001,mid;
		for(i=1;i<=n;i++)	x[i]=rd(),y[i]=rd();
		while(l<r)
		{
			mid=(l+r)>>1;
			if(check(mid))	l=mid+1;
			else	r=mid;
		}
		printf("Case %d: ",cas);
		if(l==1000001)	printf("never\n");
		else	printf("%d\n",l);
	}
}//4 -3 5 3 4 -6 -2 1 -5 1 0 -1 -1

 

以上是关于BZOJ3958[WF2011]Mummy Madness 二分+扫描线+线段树的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ3963[WF2011]MachineWorks cdq分治+斜率优化

●BZOJ 3963 [WF2011]MachineWorks

BZOJ 2673[Wf2011]Chips Challenge

BZOJ 3963: [WF2011]MachineWorks [CDQ分治 斜率优化DP]

BZOJ 4078: [Wf2014]Metal Processing Plant

BZOJ-4950: [Wf2017]Mission Improbable (二分图最大匹配)