[BZOJ2716]天使玩偶

[BZOJ2716]天使玩偶

题目大意：

一个平面直角坐标系，坐标\(1\le x,y\le10^6\)。\(n(n\le10^6)\)次操作，操作包含以下两种：

1. 新增一个点\((x,y)\)；
2. 询问离\((x,y)\)最近的点的距离。

思路：

分别统计左下、左上、右上、右下的最近的点，每次使用CDQ分治。树状数组维护最小值。

有一个一样的题是SJY摆棋子，不过CDQ会被卡，需要用KD树才能过。

源代码:

#include<cstdio>
#include<cctype>
#include<climits>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=1e6+1,M=5e5;
struct Node {
    int t,x,y,type;
};
Node a[N];
int ans[M],lim;
inline bool cmp1(const Node &p1,const Node &p2) {
    if(p1.t==p2.t) {
        if(p1.type==p2.type) {
            if(p1.x==p2.x) {
                return p1.y<p2.y;
            }
            return p1.x<p2.x;
        }
        return p1.type>p2.type;
    }
    return p1.t<p2.t;
}
inline bool cmp2(const Node &p1,const Node &p2) {
    if(p1.x==p2.x) {
        return p1.y<p2.y;
    }
    return p1.x<p2.x;
}
class FenwickTree {
    private:
        int val[N];
        int lowbit(const int &x) const {
            return x&-x;
        }
    public:
        void reset() {
            std::fill(&val[1],&val[N],INT_MAX);
        }
        void modify(int p,const int &x) {
            for(;p<N;p+=lowbit(p)) val[p]=std::min(val[p],x);
        }
        int query(int p) const {
            int ret=INT_MAX;
            for(;p;p-=lowbit(p)) ret=std::min(ret,val[p]);
            return ret;
        }
        void clear(int p) {
            for(;p<N;p+=lowbit(p)) val[p]=INT_MAX;
        }
};
FenwickTree t;
void cdq(const int &b,const int &e) {
    if(b==e) return;
    const int mid=(b+e)>>1;
    cdq(b,mid);
    cdq(mid+1,e);
    int p=b,q=mid+1;
    for(;q<=e;q++) {
        if(a[q].type==1) continue;
        for(;p<=mid&&a[p].x<=a[q].x;p++) {
            if(a[p].type==1) t.modify(a[p].y,N-a[p].y-a[p].x);
        }
        const int tmp=t.query(a[q].y);
        if(tmp==INT_MAX) continue;
        ans[a[q].t]=std::min(ans[a[q].t],tmp-(N-a[q].y-a[q].x));
    }
    while(--p>=b) {
        if(a[p].type==1) t.clear(a[p].y);
    }
    std::inplace_merge(&a[b],&a[mid]+1,&a[e]+1,cmp2);
}
int main() {
    const int n=getint(),m=getint();
    int cnt=-1;
    for(register int i=1;i<=n;i++) {
        a[i].type=1;
        a[i].t=cnt;
        a[i].x=getint();
        a[i].y=getint();
    }
    for(register int i=1;i<=m;i++) {
        a[n+i].type=getint();
        if(a[n+i].type==2) cnt++;
        a[n+i].t=cnt;
        a[n+i].x=getint();
        a[n+i].y=getint();
        ans[cnt]=INT_MAX;
    }
    t.reset();
    std::sort(&a[1],&a[n+m]+1,cmp1);
    cdq(1,n+m);
    for(register int i=1;i<=n+m;i++) {
        a[i].x=N-a[i].x;
    }
    std::sort(&a[1],&a[n+m]+1,cmp1);
    cdq(1,n+m);
    for(register int i=1;i<=n+m;i++) {
        a[i].y=N-a[i].y;
    }
    std::sort(&a[1],&a[n+m]+1,cmp1);
    cdq(1,n+m);
    for(register int i=1;i<=n+m;i++) {
        a[i].x=N-a[i].x;
    }
    std::sort(&a[1],&a[n+m]+1,cmp1);
    cdq(1,n+m);
    for(register int i=0;i<=cnt;i++) {
        printf("%d\n",ans[i]);
    }
    return 0;
}