poj2464扫描线好题，回头再做

扫描线+区间更新

题解

/*
st[i],ol[i]表示y坐标大于y[i]和小于y[i]的点
两颗线段树建立在y轴上，区间[l,r]ol线选在[l,r]时st的分数
每次查询完成后再更新一次
遍历每条st线上的ol线，这是单点查询
*/
#include<iostream>
#include<cstring>
#include<cstdio>
#include<map>
#include<vector>
#include<algorithm>
using namespace std;
#define N 200005
#define LL(x) (x<<1)
#define RR(x) (x<<1|1)
#define MID(a,b) (a+((b-a)>>1))

struct Point{
    int x,y;
    void get(){scanf("%d%d",&x,&y);}
    bool operator<(const Point & b)const
    { return x<b.x; }
} p[N];
struct node{//线段树结点
    int lft,rht;
    int flag[2];
    int mid(){return MID(lft,rht);}
    void init(int a,int b){flag[0]=a;flag[1]=b;}
};
int n,m,mi;
vector<int> mx;
map<int,int> H;
int y[N],st[N],ol[N];//离散化数轴，大于等于y[i]的点数，小于等于y[i]的点
struct Segtree{
    node tree[N*4];
    void down(int ind){
        for(int i=0;i<2;i++){
            tree[LL(ind)].flag[i]+=tree[ind].flag[i];
            tree[RR(ind)].flag[i]+=tree[ind].flag[i];
            tree[ind].flag[i]=0;
        }
    }
    void build(int lft,int rht,int ind){
        tree[ind].lft=lft;tree[ind].rht=rht;
        tree[ind].init(0,0);
        if(lft==rht) tree[ind].init(st[lft],ol[lft]);
        else {
            int mid=tree[ind].mid();
            build(lft,mid,LL(ind));
            build(mid+1,rht,RR(ind));
        }
    }
    void update(int st,int ed,int ind,int type,int valu){
        int lft=tree[ind].lft,rht=tree[ind].rht;
        if(st<=lft && ed>=rht) 
            tree[ind].flag[type]+=valu;
        else {
            down(ind);
            int mid=tree[ind].mid();
            if(st<=mid) update(st,ed,LL(ind),type,valu);
            if(ed>mid) update(st,ed,LL(ind),type,valu);
        }
    }
    void query(int pos,int ind,int &mi,int &mx){
        if(tree[ind].lft==tree[ind].rht){
            mi=tree[ind].flag[0];
            mx=tree[ind].flag[1];
        }
        else {
            down(ind);
            int mid=tree[ind].mid();
            if(pos<=mid) return query(pos,LL(ind),mi,mx);
            if(pos>mid) return query(pos,RR(ind),mi,mx);
        }
    }
}seg;
int main(){
    while(scanf("%d",&n),n){
        H.clear();mx.clear();mi=m=0;

        int id1=0,id2=0;
        for(int i=0;i<n;i++){
            p[i].get();
            y[i]=p[i].y;
        }
        sort(y,y+n);
        sort(p,p+n);//把点按从左往右顺序排好
        H[y[0]]=0;
        for(int i=0;i<n;i++)//遍历一次纵坐标
            if(y[m]!=y[i]){
                st[m]=n-i;//高于y[i]的点数
                y[++m]=y[i];
                ol[m]=i;//低于y[i]的点数
                H[y[m]]=m;
            }
        st[m]=0;
        seg.build(0,m,1);
        while(id1<n){
            id2=id1;
            while(p[id1].x==p[id2].x){//这一步删掉被st线穿过的点
                //不是最低的点
                if(p[id2].y!=y[0]) //ol线低于p[id2].y情况的st分数就少了一份
                    seg.update(H[y[0]],H[p[id2].y]-1,1,0,-1);
                //不是最高的点
                if(p[id2].y!=y[m]) //ol线高于p[id2].y的时候
                    seg.update(H[p[id2].y]+1,H[y[m]],1,1,-1);
                if(++id2>=n) break;
            }
            int mii=n,mxx=0;
            for(int i=id1;i<id2;i++){//依次在st线上的点建立ol线并进行查询，找到使st得分最低的那个点
                int tmp1,tmp2;
                seg.query(H[p[i].y],1,mii,mxx);
                mii=min(mii,tmp1);
                mxx=max(mxx,tmp2);
            }

            if(mii==mi) mx.push_back(mxx);//找到了ol的新解
            else if(mii>mi){//找到了st的更优解
                mi=mii;
                mx.clear();mx.push_back(mxx);
            }
            //再更新一次，把这条线上删掉的反向加回去：因为扫描线往后扫描这条线上原来应该属于ol的现在属于st，原来属于st的现在应该属于ol
            for(int i=id1;i<id2;i++){
                if(p[i].y!=y[0])
                    seg.update(H[y[0]],H[p[i].y]-1,1,1,1);
                if(p[i].y!=y[m])
                    seg.update(H[p[i].y]+1,H[y[m]],1,0,1);
            }
            id1=id2;
        }
        sort(mx.begin(),mx.end());
        mx.erase(unique(mx.begin(),mx.end()),mx.end());
        printf("Stan: %d; Ollie:",mi);
        for(int i=0;i<mx.size();i++)
            printf(" %d",mx[i]);
        printf(";
");
    }
    return 0;
}