bzoj2648

http://www.lydsy.com/JudgeOnline/problem.php?id=2648

kdtree裸题。。。

抄板子一边抄对了。。。

挺好理解的，就是说我们先找出中间的元素，然后小的放左边，大的放右边，这样我们就确定了这个点，然后维度换一下，递归建树，很像splay，其实就是把平面每次一分为二，先横着再竖着。。。

然后那个get挺神奇的，大概就是估计一下现在这个点距离下一个区域的距离？其实就是优化枚举，减少了很多重复状态。。。

#include<bits/stdc++.h>
using namespace std;
const int N = , inf =  << ;
struct data {
    int l, r, d[], min[], max[];
} t[N];
int n, m, root, ans, d;
bool cp(data x, data y)
{
    if(x.d[d] != y.d[d]) return x.d[d] < y.d[d];
    return x.d[d ^ ] < y.d[d ^ ];
}
void update(int k)
{
    for(int i = ; i < ; ++i)
    {
        if(t[k].l) t[k].min[i] = min(t[k].min[i], t[t[k].l].min[i]), t[k].max[i] = max(t[k].max[i], t[t[k].l].max[i]);
        if(t[k].r) t[k].min[i] = min(t[k].min[i], t[t[k].r].min[i]), t[k].max[i] = max(t[k].max[i], t[t[k].r].max[i]);
    }
}
int get(int k, data x)
{
    int ret = ;
    for(int i = ; i < ; ++i) ret += max(t[k].min[i] - x.d[i], ) + max(x.d[i] - t[k].max[i], );
    return ret;
}
int build(int l, int r, int now)
{
    d = now;
    int mid = (l + r) >> ;
    nth_element(t + l + , t + mid + , t + r + , cp);
    for(int i = ; i < ; ++i) t[mid].min[i] = t[mid].max[i] = t[mid].d[i];
    if(l < mid) t[mid].l = build(l, mid - , now ^ );
    if(r > mid) t[mid].r = build(mid + , r, now ^ );
    update(mid); return mid;
}
void insert(int k, int now, data x)
{
    if(x.d[now] < t[k].d[now])
    {
        if(t[k].l) insert(t[k].l, now ^ , x);
        else
        {
            t[k].l = ++n;
            for(int i = ; i < ; ++i) t[n].max[i] = t[n].min[i] = t[n].d[i] = x.d[i];
        }
    }
    else
    {
        if(t[k].r) insert(t[k].r, now ^ , x);
        else
        {
            t[k].r = ++n;
            for(int i = ; i < ; ++i) t[n].max[i] = t[n].min[i] = t[n].d[i] = x.d[i];
        }
    }
    update(k);
}
void query(int k, int now, data x)
{
    int d = abs(t[k].d[] - x.d[]) + abs(t[k].d[] - x.d[]), dl = inf, dr = inf;
    ans = min(ans, d);
    if(t[k].l) dl = get(t[k].l, x);
    if(t[k].r) dr = get(t[k].r, x);
    if(dl < dr)
    {
        if(dl < ans) query(t[k].l, now ^ , x);
        if(dr < ans) query(t[k].r, now ^ , x);
    }
    else
    {
        if(dr < ans) query(t[k].r, now ^ , x);
        if(dl < ans) query(t[k].l, now ^ , x);
    }
}
int main()
{
    scanf("%d%d", &n, &m);
    for(int i = ; i <= n; ++i) scanf("%d%d", &t[i].d[], &t[i].d[]);
    root = build(, n, );
    while(m--)
    {
        int opt; data x; scanf("%d%d%d", &opt, &x.d[], &x.d[]);
        if(opt == ) insert(root, , x);
        if(opt == )
        {
            ans = inf;
            query(root, , x);
            printf("%d\n", ans);
        }
    }
    return ;
}