BZOJ3958: [WF2011]Mummy Madness

可以发现，你被抓住的时间满足二分性，因为若在时间x你被抓住了，那么对于任意时间y>x，木乃伊都一定能抓住你

那么二分这个时间，你能走到的地方是一个矩形，同时每个木乃伊能走到的地方也是一个矩形，扫描线扫一下，若整个矩形都被覆盖了就会被抓住

（只处理每个木乃伊的矩形和你走到的矩形的交的部分，不然可能会T）

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

inline void read(int &x)
{
char c; int f=1;
while(!((c=getchar())>='0'&&c<='9')) if(c=='-') f=-f;
x=c-'0';
while((c=getchar())>='0'&&c<='9') x=(x<<3)+(x<<1)+c-'0';
x*=f;
}
inline void up(int &x,const int y){if(x<y)x=y;}
inline int _max(const int &x,const int &y){return x>y?x:y;}
inline int _min(const int &x,const int &y){return x<y?x:y;}
const int maxn = 210000;
const int maxc = 1010000;

int n,C;
struct P{int x,y;}p[maxn];

struct edge
{
int l,r,c,nex;
edge(){}
edge(const int _l,const int _r,const int _c,const int _nex){l=_l;r=_r;c=_c;nex=_nex;}
}a[maxn<<1]; int len,fir[maxc<<1];
inline void ins(const int x,const int l,const int r,const int c)
{
a[++len]=edge(l,r,c,fir[x]);fir[x]=len;
}

struct segment
{
int flag,c;
}seg[maxc<<4];
int lx,rx,c;
void upd(const int x,const int l,const int r)
{
if(lx<=l&&r<=rx) { seg[x].flag+=c; seg[x].c+=c; return; }
if(rx<l||r<lx) return;
upd(x<<1,l,l+r>>1),upd((x<<1)|1,(l+r>>1)+1,r);
seg[x].c=seg[x].flag+_min(seg[x<<1].c,seg[(x<<1)|1].c);
}

bool Judge(const int m)
{
len=0; for(int i=0;i<=2*C+1;i++) fir[i]=0;

for(int i=1;i<=n;i++)
{
int d=m+max(p[i].y-m,-m),u=m+min(p[i].y+m,m);
int l=m+max(p[i].x-m,-m),r=m+min(p[i].x+m,m);
if(u>=0&&d<=u) ins(d,l,r,1),ins(u+1,l,r,-1);
}

bool flag=false;
for(int i=0;i<=2*m;i++)
{
for(int k=fir[i];k;k=a[k].nex)
{
lx=a[k].l,rx=a[k].r,c=a[k].c;
if(lx<=rx) upd(1,0,m<<1);
}
if(!seg[1].c) flag=true;
}
for(int k=fir[2*m+1];k;k=a[k].nex)
{
lx=a[k].l,rx=a[k].r,c=a[k].c;
if(lx<=rx) upd(1,0,m<<1);
}
return flag;
}
int solve()
{
int l=1,r=C;
while(l<=r)
{
int mid=l+r>>1;
if(Judge(mid)) l=mid+1;
else r=mid-1;
}
return l-1;
}

int main()
{
for(int t=1;scanf("%d",&n)&&n!=-1;t++)
{
C=0;
for(int i=1;i<=n;i++)
{
int x,y; read(x); read(y);
p[i].x=x; p[i].y=y;
up(C,abs(x)); up(C,abs(y));
}
int re=solve()+1;
printf("Case %d: ",t);
if(re==C+1) printf("never\n");
else printf("%d\n",re);
}

return 0;
}