CF 19D Points 【线段树+平衡树】

在平面上进行三种操作：

1、add x y：在平面上添加一个点（x，y）

2、remove x y：将平面上的点（x，y）删除

3、find x y：在平面上寻找一个点，使这个点的横坐标大于x，纵坐标大于y，而且要求他的横坐标尽量小，如果有多个点满足，则选取横坐标尽量小的前提下，纵坐标最小的点。

方法：

将横坐标x离散化，每一个坐标x对应的y用一颗平衡树维护（C++中的set），则这颗平衡树支持增加和删除以及查找比y大的最小值的操作。

在此基础上，对于每一个询问，只需要遍历大于x的set，并且找到最小的y即可。但是这样子依次向后遍历的复杂度为O(N)的，所以要用一颗线段树来维护横坐标区间段里面y的最大值，这样就可以在O(logN)的复杂度内找到最接近x的横坐标（而且满足当前坐标的最大纵坐标大于y），然后再set里面查找最小的纵坐标即可。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <set>
#define lowbit(x) ((x)&(-(x)))
#define lson (step<<1)
#define rson (step<<1|1)
using namespace std;
const int N = 200005;
struct Query{
    int k,x,y;
    void init(){
        char str[10];
        scanf("%s%d%d",str,&x,&y);
        if(str[0]=='a') k=0;
        else if(str[0]=='f') k=2;
        else k=1;
    }
}ask[N];
struct Seg_tree{
    int left,right,mx;
}L[N<<2];
int m,q,X[N];
set<int>s[N];
void bulid(int step,int l,int r){
    L[step].left=l;
    L[step].right=r;
    L[step].mx=-1;
    if(l==r) return ;
    int m=(l+r)>>1;
    bulid(lson,l,m);
    bulid(rson,m+1,r);
}
inline void push_up(int step){
    L[step].mx=max(L[lson].mx,L[rson].mx);
}
void update(int step,int pos){
    if(L[step].left==L[step].right){
        L[step].mx=((int)s[pos].size()==0)?-1:(*(--s[pos].end()));
        return ;
    }
    int m=(L[step].left+L[step].right)>>1;
    if(pos<=m) update(lson,pos);
    else update(rson,pos);
    push_up(step);
}
int query(int step,int l,int r,int y){
    if(l>r) return -1;
    if(L[step].mx<=y) return -1;
    if(L[step].left==L[step].right) return L[step].left;
    int m=(L[step].left+L[step].right)>>1;
    if(r<=m) return query(lson,l,r,y);
    else if(l>m) return query(rson,l,r,y);
    else{
        int t=query(lson,l,m,y);
        if(t>=0) return t;
        return query(rson,m+1,r,y); 
    }
}
int main(){
    scanf("%d",&q);
    for(int i=0;i<q;i++){
        ask[i].init();    
        X[i]=ask[i].x;
    }
    sort(X,X+q);
    m=unique(X,X+q)-X;
    bulid(1,1,m);
    for(int i=0;i<q;i++){
        int k=ask[i].k,x=lower_bound(X,X+m,ask[i].x)-X+1,y=ask[i].y;
        if(k==0){
            s[x].insert(y);
            update(1,x);
        }
        else if(k==1){
            s[x].erase(s[x].find(y));
            update(1,x);
        }
        else{
            int pos=query(1,x+1,m,y);
            if(pos==-1) printf("-1\n");
            else{
                printf("%d %d\n",X[pos-1],*s[pos].upper_bound(y));
            }
        }
    }
    return 0;
}