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 MBSubmit: 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
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
