BZOJ 1941: [Sdoi2010]Hide and Seek KDtree + 估价函数

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

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

题解： 求对于每个点最远距离-最近距离的最小值， 依次枚举距离每个点最远距离与最近距离即可.

#include<bits/stdc++.h>
#define maxn 200000
#define inf 1000000000 
using namespace std;
void getmin(int &a,int b) { if(b<a) a=b;  }
void getmax(int &a,int b) { if(b>a) a=b;  }
void setIO(string s)
{
    string in=s+".in"; 
    string out=s+".out"; 
    freopen(in.c_str(),"r",stdin); 
    freopen(out.c_str(),"w",stdout);   
}
namespace KDtree
{
    #define lson t[x].ch[0]
    #define rson t[x].ch[1] 
    #define mid ((l+r)>>1)
    int d,ansmin,ansmax,answer;      
    struct Node
    {
        int ch[2],minv[2],maxv[2],p[2]; 
    }t[maxn],T;   
    bool cmp(Node a,Node b)
    {
        if(a.p[d] == b.p[d]) return a.p[d^1] < b.p[d^1]; 
        return a.p[d] < b.p[d]; 
    }        
    void pushup(int o,int x)
    {
        getmin(t[o].minv[0], t[x].minv[0]); 
        getmax(t[o].maxv[0], t[x].maxv[0]); 
        getmin(t[o].minv[1], t[x].minv[1]); 
        getmax(t[o].maxv[1], t[x].maxv[1]);                           
    }
    int build(int l,int r,int o)
    {
        d=o; 
        nth_element(t+l,t+mid,t+1+r,cmp); 
        t[mid].minv[0]=t[mid].maxv[0]=t[mid].p[0]; 
        t[mid].minv[1]=t[mid].maxv[1]=t[mid].p[1];  
        t[mid].ch[0]=t[mid].ch[1]=0; 
        if(l<mid) 
        {
            t[mid].ch[0]=build(l,mid-1,o^1); 
            pushup(mid,t[mid].ch[0]); 
        }
        if(r>mid) 
        {
            t[mid].ch[1]=build(mid+1,r,o^1); 
            pushup(mid,t[mid].ch[1]);         
        }
        return mid;   
    }
    int _min(Node a)
    {
        int ans=0;                      
        for(int i=0;i<2;++i)
        {
            ans += max(T.p[i]-a.maxv[i], 0); 
            ans += max(a.minv[i]-T.p[i], 0);            
        }
        return ans; 
    }
    int _max(Node a)
    {
        int ans=0; 
        for(int i=0;i<2;++i)
        {
            ans += max(abs(T.p[i]-a.minv[i]), abs(T.p[i]-a.maxv[i])); 
        }
        return ans; 
    }
    void qmin(int x,int x1,int y1)
    {
        int dn=abs(t[x].p[0]-x1) + abs(t[x].p[1]-y1),dl,dr; 
        if(dn>0) getmin(ansmin,dn); 
        dl=lson?_min(t[lson]):inf; 
        dr=rson?_min(t[rson]):inf;   
        if(dl<dr)
        {
            if(dl<ansmin) qmin(lson,x1,y1);        
            if(dr<ansmin) qmin(rson,x1,y1); 
        }
        else
        {
            if(dr<ansmin) qmin(rson,x1,y1); 
            if(dl<ansmin) qmin(lson,x1,y1);        
        }
    }
    void qmax(int x,int x1,int y1)
    {
        int dn=abs(t[x].p[0]-x1) + abs(t[x].p[1]-y1),dl,dr; 
        getmax(ansmax, dn); 
        dl=lson?_max(t[lson]):-inf; 
        dr=rson?_max(t[rson]):-inf; 
        if(dl>dr) 
        {
            if(dl>ansmax) qmax(lson,x1,y1); 
            if(dr>ansmax) qmax(rson,x1,y1); 
        }
        else 
        {
            if(dr>ansmax) qmax(rson,x1,y1); 
            if(dl>ansmax) qmax(lson,x1,y1); 
        }
    }
    void solve(int i,int rt)
    {
        T=t[i], ansmin=inf,ansmax=-inf;   
        qmin(rt, T.p[0], T.p[1]); 
        qmax(rt, T.p[0], T.p[1]); 
        getmin(answer, ansmax-ansmin);              
    }
}; 
int main()
{
    //setIO("A");       
    int n,root,i;
    KDtree::answer=inf; 
    scanf("%d",&n);          
    for(i=1;i<=n;++i) 
        scanf("%d%d",&KDtree::t[i].p[0],&KDtree::t[i].p[1]);  
    root=KDtree::build(1,n,0); 
    for(i=1;i<=n;++i)      
    { 
        KDtree::solve(i,root); 
    }      
    printf("%d\n",KDtree::answer); 
    return 0; 
}