BZOJ1941: [Sdoi2010]Hide and Seek

题解:KDtree模板题 估价函数决策最右两边即可

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <stack>
#include <queue>
#include <cmath>
#include <set>
#include <map>
#define mp make_pair
#define pb push_back
#define pii pair<int,int>
#define link(x) for(edge *j=h[x];j;j=j->next)
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
const int MAXN=5e5+10;
const double eps=1e-8;
const int inf=1e9;
#define ll long long
using namespace std;
struct edge{int t,v;edge*next;}e[MAXN<<1],*h[MAXN],*o=e;
void add(int x,int y,int vul){o->t=y;o->v=vul;o->next=h[x];h[x]=o++;}
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int d,rt,n;
typedef struct node{
	int p[2],maxx[2],minn[2],c[2];
	friend bool operator<(node aa,node bb){
		if(aa.p[d]!=bb.p[d])return aa.p[d]<bb.p[d];
		return aa.p[d^1]<bb.p[d^1];
	}
}node;
node a[MAXN];
void up(int x,int y){
	for(int i=0;i<=1;i++)a[x].maxx[i]=max(a[x].maxx[i],a[y].maxx[i]),a[x].minn[i]=min(a[x].minn[i],a[y].minn[i]);
}
int built(int l,int r,int now){
	int mid=(l+r)>>1;
	d=now;nth_element(a+l,a+mid,a+r+1);
	for(int i=0;i<=1;i++)a[mid].maxx[i]=a[mid].minn[i]=a[mid].p[i];
	a[mid].c[0]=a[mid].c[1]=0;
	if(l<mid)a[mid].c[0]=built(l,mid-1,now^1),up(mid,a[mid].c[0]);
	if(r>mid)a[mid].c[1]=built(mid+1,r,now^1),up(mid,a[mid].c[1]);
	return mid;
}
int ans;node t;
int dist(node x,node y){
	int res=0;
	for(int i=0;i<=1;i++)res+=abs(x.p[i]-y.p[i]);
	return res;
}
int get_min(node x,node y){
	int res=0;
	for(int i=0;i<=1;i++){
		if(y.p[i]<x.minn[i])res+=x.minn[i]-y.p[i];
		if(y.p[i]>x.maxx[i])res+=y.p[i]-x.maxx[i];
	}
	return res;
}
int get_max(node x,node y){
	int res=0;
	for(int i=0;i<=1;i++){
		res+=max(abs(y.p[i]-x.minn[i]),abs(y.p[i]-x.maxx[i]));
	}
	return res;
}
void Min(int x){
	if(!x)return ;
	if(a[x].p[0]!=t.p[0]||a[x].p[1]!=t.p[1])ans=min(ans,dist(a[x],t));
	int dl=(a[x].c[0])?get_min(a[a[x].c[0]],t):inf;
	int dr=(a[x].c[1])?get_min(a[a[x].c[1]],t):inf;
	if(dl<dr){
		if(dl<ans)Min(a[x].c[0]);
		if(dr<ans)Min(a[x].c[1]);
	}
	else{
		if(dr<ans)Min(a[x].c[1]);
		if(dl<ans)Min(a[x].c[0]);
	}
}
void Max(int x){
	if(!x)return ;
	if(a[x].p[0]!=t.p[0]||a[x].p[1]!=t.p[1])ans=max(ans,dist(a[x],t));
	int dl=(a[x].c[0])?get_max(a[a[x].c[0]],t):-inf;
	int dr=(a[x].c[1])?get_max(a[a[x].c[1]],t):-inf;
	if(dl>dr){
		if(dl>ans)Max(a[x].c[0]);
		if(dr>ans)Max(a[x].c[1]);
	}
	else{
		if(dr>ans)Max(a[x].c[1]);
		if(dl>ans)Max(a[x].c[0]);
	}
}
int main(){
	n=read();
	inc(i,1,n)inc(j,0,1)a[i].p[j]=read();
	rt=built(1,n,0);
	int res=inf,res1,res2;
	inc(i,1,n){
		ans=inf;t=a[i];Min(rt);res1=ans;
		ans=-inf;t=a[i];Max(rt);res2=ans;
		res=min(res,res2-res1);
	}
	printf("%d\n",res);
}

　　

1941: [Sdoi2010]Hide and Seek

Time Limit: 16 Sec  Memory Limit: 162 MB
Submit: 1846  Solved: 1004
[Submit][Status][Discuss]

Description

小猪iPig在PKU刚上完了无聊的猪性代数课，天资聪慧的iPig被这门对他来说无比简单的课弄得非常寂寞，为了消除寂寞感，他决定和他的好朋友giPi（鸡皮）玩一个更加寂寞的游戏---捉迷藏。 但是，他们觉得，玩普通的捉迷藏没什么意思，还是不够寂寞，于是，他们决定玩寂寞无比的螃蟹版捉迷藏，顾名思义，就是说他们在玩游戏的时候只能沿水平或垂直方向走。一番寂寞的剪刀石头布后，他们决定iPig去捉giPi。由于他们都很熟悉PKU的地形了，所以giPi只会躲在PKU内n个隐秘地点，显然iPig也只会在那n个地点内找giPi。游戏一开始，他们选定一个地点，iPig保持不动，然后giPi用30秒的时间逃离现场（显然，giPi不会呆在原地）。然后iPig会随机地去找giPi，直到找到为止。由于iPig很懒，所以他到总是走最短的路径，而且，他选择起始点不是随便选的，他想找一个地点，使得该地点到最远的地点和最近的地点的距离差最小。iPig现在想知道这个距离差最小是多少。 由于iPig现在手上没有电脑，所以不能编程解决这个如此简单的问题，所以他马上打了个电话，要求你帮他解决这个问题。iPig告诉了你PKU的n个隐秘地点的坐标，请你编程求出iPig的问题。

Input

第一行输入一个整数N 第2~N+1行，每行两个整数X，Y，表示第i个地点的坐标

Output

一个整数，为距离差的最小值。

Sample Input

4
0 0
1 0
0 1
1 1

Sample Output

1

HINT

 

对于30%的数据,N<=1000 对于100%的数据，N<=500000,0<=X,Y<=10^8 保证数据没有重点保证N>=2