hdu 1007 最近点对问题（Splay解法）

为什么要写这个题、、经典啊，当然，别以为我用分治做的，不过主要思想还是那神奇的六个点共存（一个h*2h的矩形中最多能放下多少个点使得两两距离不超过h）

其实我是在这里看到的

http://community.topcoder.com/tc?module=Static&d1=tutorials&d2=lineSweep

排个序，然后扫描过去，每次确定Y的范围，暴力找每个点（其实这是O(1)的）

蛮不错的哦。

写完后发现比分治写的快了好多啊，估计是我不会写分治吧T_T，总之，现在跑到了rank 4了。。。

数据结构就是这样的直接啊！

 

/* **********************************************
Author      : wuyiqi
Created Time: 2013-8-24 12:54:35
File Name   : ruocai.cpp
*********************************************** */
#include <cstdio>
#include <cmath>
#include <stack>
#include <algorithm>
using std::stack;
    using std::pair;
    using std::make_pair;
#define L x->c[0]
#define R x->c[1]
#define KT root->c[1]->c[0]
const int maxn = 200010;
struct node{
    node *c[2] , *fa;
    int id,sz;
    double val,x;
    inline bool d() {
        return fa->c[0] == this;
    }
    inline void setc(int d,node *s) {
        c[d] = s;
        s->fa = this;
    }
    inline void up() {
        sz = c[0]->sz + c[1]->sz + 1;
    }
}NODE[maxn],*null=&NODE[0];
node *ID[maxn];
int top;
struct Splay_tree{
    node *root;
    void Rotate(node *x,int f){    
        node *y = x->fa;    
        y->setc(!f,x->c[f]);    
        x->fa = y->fa;    
        if(y->fa != null) y->fa->setc(!y->d(),x);    
        x->setc(f,y);    
        y->up();    
    }    
    void Splay(node *x,node *goal) {    
        while(x->fa!=goal) {    
            if(x->fa->fa == goal) Rotate(x,x->d());    
            else {    
                int f = x->fa->d();    
                x->d() == f ? Rotate(x->fa,f) : Rotate(x,!f);    
                Rotate(x,f);    
            }    
        }    
        x->up();    
        if(goal == null)  {    
            root = x;    
        }    
    }    
    void RTO(int k,node *goal) {    
        node *x = root;    
        while(L->sz + 1 != k) {    
            if(k < L->sz + 1) x = L;    
            else {    
                k -= L->sz + 1;    
                x = R;    
            }    
        }    
        Splay(x,goal);    
    }    
    void insert(node* &x,node *y) {    
        if(x == null) {    
            x = y;    
            return ;    
        }    
        if(y->val <= x->val) {    
            insert(x->c[0],y);    
            x->c[0]->fa = x;    
        } else {    
            insert(x->c[1],y);    
            x->c[1]->fa = x;    
        }    
        x->up();    
    }    
    node* new_node(int id,double x,double y) {
        node* tmp = &NODE[++top];
        tmp->fa = null;
        tmp->val = y;
        tmp->x = x;
        tmp->sz = 1;
        tmp->id = top;
        tmp->c[0] = tmp->c[1] = null;
        ID[id] = tmp;
        return tmp;
    }
    void insert(int id,double x,double y) {
        node *tmp = new_node(id,x,y);
        insert(root,tmp);
    }
    void init() {
        root = null;
    }
    void Del_root() {    
        node *t = root;    
        if(t->c[1] != null) {    
            root = t->c[1];    
            RTO(1,null);    
            root->c[0] = t->c[0];    
            if(root->c[0] != null)     
                root->c[0]->fa = root;    
        } else  {    
            root = root->c[0];    
        }    
        root->fa = null;    
        if(root != null) root->up();    
    }    

    void Del(node *tmp) {
        Splay(tmp,null);
        Del_root();
    }
    node* find_pre(node *x,double v) {
        if(x == null) return null;
        if(x->val < v) {
            node *tmp = find_pre(x->c[1],v);
            return tmp == null ? x : tmp;
        }
        else {
           return  find_pre(x->c[0],v);
        }
    }
    node* find_succ(node *x,double v) {
        if(x == null) return null;
        if(x->val > v ) {
            node *tmp = find_succ(x->c[0],v);
            return tmp == null ? x : tmp;
        } else {
            return find_succ(x->c[1],v);
        }
    }
    void search(double a,double b,double &ans,node* x){
        if(x == null) return ;
        double  c = x->x , d = x->val;
        double cost = sqrt((a-c)*(a-c) + (b-d)*(b-d));
        if(cost  < ans) ans = cost;
        if(x->c[0] != null) search(a,b,ans,x->c[0]);
        if(x->c[1] != null) search(a,b,ans,x->c[1]);
    }
    void gao(double x,double y,double &ans) {
        node *pre = find_pre(root,y-ans) ;
        node *succ = find_succ(root,y+ans);
        if(pre == null || succ == null) while(1);
        Splay(pre,null);
        Splay(succ,root);
        search(x,y,ans,KT);
    }
    void vist(node *x) {
        if(x!=null) {
            printf("%d x = %lf y = %lf lson=%d rson=%d\n",x->id,x->x,x->val,x->c[0]->id,x->c[1]->id);
            if(x->c[0]!=null) vist(x->c[0]);
            if(x->c[1]!=null) vist(x->c[1]);
        }
    }
    void debug(){
        puts("*******");
        vist(root);
        puts("*****");
    }
}spt;
void prepare() {
    null->id = 0;
    null->c[0] = null->c[1] = null->fa = NULL;
    null->sz = null->val = 0;
    top = 0;
}
double x[maxn] , y[maxn];
bool by_x(int a,int b) {
    return x[a] > x[b];
}
int id[maxn];
double sqr(double a) {  return a * a; }
int main()
{
    int n;
    while(scanf("%d",&n),n) {
        prepare();
        spt.init();
        for(int i = 0; i < n; i++) {
            id[i] = i;
            scanf("%lf%lf",&x[i],&y[i]);
        }
        if(n == 1){
            puts("0.00");
            continue;
        }
        std::sort(id,id+n,by_x);
        spt.insert(0,-1e100,-1e100);
        spt.insert(0,-1e100,1e100);
        double ans = 1e50 ;
        int pt = 0;
        for(int i = 0; i < n; i++) {
            if(ans == 0) break;
            while(pt < i && x[id[pt]] > x[id[i]] + ans ) {
                if(spt.root->sz == 2) break;
                spt.Del(ID[id[pt]]);
                pt++;
            }
            if(spt.root->sz == 2) {
                spt.insert(id[i],x[id[i]],y[id[i]]);
            } else {
                spt.gao(x[id[i]],y[id[i]],ans);
                spt.insert(id[i],x[id[i]],y[id[i]]);
            }
        }
        printf("%.2f\n",ans*0.5);
    }
    return 0;
}