Codeforces Beta Round #19D. Points

链接：http://www.codeforces.com/contest/19/problem/D

题意：给定n个操作，add x y：在平面坐标系中[x,y]位置插入一个点，remove x y：将平面中[x,y]位置的点删掉，find x y：输出所有[x',y']，x'>x,y'>y中x'最小的那个点，如果有多个点输出y'最小的那个点。

分析：我们先离线将所有插入点和询问点都处理出来，将它们以x为第一关键字y为第二关键字从小到大排序。如果插入点的话，我们就找到这个点对应的位置标记有该点，查询的话我们只需要在大于x的范围中找y'>y的最小位置即可。比如我们现在查询[x,y]，那么我们只需要在排好序的数字中找到第一个大于x的位置开始，找到第一个y'>y的位置即可。这里找最左边的y'直接用线段树维护一下即可。O(nlogn)。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=200010;
const int MAX=1000000100;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=998244353;
const int INF=1000000010;
typedef double db;
typedef unsigned long long ull;
struct node {
    int op,x,y;
    node() {}
    node(int op,int x,int y):op(op),x(x),y(y) {}
    bool operator < (const node a) const {
        if (x!=a.x) return x<a.x;
        return y<a.y;
    }
}ope[N],add[N];
char s[5];
int k,sum[4*N],mx[4*N];
void build(int x,int l,int r) {
    sum[x]=mx[x]=0;
    if (l==r) return ;
    int mid=(l+r)>>1;
    build(2*x,l,mid);
    build(2*x+1,mid+1,r);
}
int find_point(int x,int y) {
    int l=0,r=k,mid=k>>1;;
    while (l+1<r)
    if (add[mid].x>x||(add[mid].x==x&&add[mid].y>=y)) { r=mid;mid=(l+r)>>1; }
    else { l=mid;mid=(l+r)>>1; }
    return r;
}
void add_point(int x,int l,int r,int w,int z) {
    if (l==r) { mx[x]=z;return ; }
    int mid=(l+r)>>1;
    if (w<=mid) add_point(2*x,l,mid,w,z);
    else add_point(2*x+1,mid+1,r,w,z);
    mx[x]=max(mx[2*x],mx[2*x+1]);
}
void del_point(int x,int l,int r,int w) {
    if (l==r) { mx[x]=0;return ; }
    int mid=(l+r)>>1;
    if (w<=mid) del_point(2*x,l,mid,w);
    else del_point(2*x+1,mid+1,r,w);
    mx[x]=max(mx[2*x],mx[2*x+1]);
}
int get(int x) {
    int l=1,r=k+1,mid=(l+r)>>1;
    while (l+1<r)
    if (add[mid].x<=x) { l=mid;mid=(l+r)>>1; }
    else { r=mid;mid=(l+r)>>1; }
    return r;
}
int query(int x,int l,int r,int L,int R,int w) {
    if (L>R||mx[x]<=w) return 0;
    if (l==r) return l;
    int ret,mid=(l+r)>>1;
    if (R<=mid) return query(2*x,l,mid,L,R,w);
    else if (L>mid) return query(2*x+1,mid+1,r,L,R,w);
        else {
            ret=query(2*x,l,mid,L,mid,w);
            if (!ret) ret=query(2*x+1,mid+1,r,mid+1,R,w);
            return ret;
        }
}
int main()
{
    int i,n,x,y,z;
    scanf("%d", &n);
    for (i=1;i<=n;i++) {
        scanf("%s%d%d", s, &x, &y);
        if (s[0]=='a') {
            ope[i]=node(1,x,y);add[++k]=node(1,x,y);
        } else if (s[0]=='r') ope[i]=node(2,x,y);
            else {
                ope[i]=node(3,x,y);add[++k]=node(1,x,y);
            }
    }
    sort(add+1,add+k+1);
    build(1,1,k);
    for (i=1;i<=n;i++)
    if (ope[i].op==1) add_point(1,1,k,find_point(ope[i].x,ope[i].y),ope[i].y);
    else if (ope[i].op==2) del_point(1,1,k,find_point(ope[i].x,ope[i].y));
        else {
            z=query(1,1,k,get(ope[i].x),k,ope[i].y);
            if (!z) printf("-1\n");
            else printf("%d %d\n", add[z].x, add[z].y);
        }
    return 0;
}