线段树专题：poj 2464 Brownie Points II

非常金典的线段树题目~~~

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ls rt<<1
#define rs rt<<1|1
#define lson l,m,ls
#define rson m+1,r,rs
#define sf scanf 
#define pf printf
using namespace std;
//非常经典的线段树题目，题意：平面上有n个点，选中一个点将平面分成四个象限，统计
//一三象限与二四象限点的个数（坐标轴上的不算），求选中一个点划分，使得一三象限点
//的个数最多，及一三点数最多的情况下，二四点数的可能值。 

const int N=200000+10;
int aby[4*N][2];
int a[N],b[N],y[N],q[N];
struct node{
       int x,y;
}p[N];

bool cmp(node a,node b){
     return a.x<b.x;
}
 
void pushdown(int rt,int o){ 
     if(aby[rt][o]!=0){
         aby[ls][o]+=aby[rt][o];
         aby[rs][o]+=aby[rt][o];
         aby[rt][o]=0;
     }
}

void build(int l, int r ,int rt){
     aby[rt][0]=aby[rt][1]=0;
     if(l==r){
         aby[rt][0]=a[l];
         aby[rt][1]=b[l];
         return;
     }
     int m=(l+r)>>1;
     build(lson);
     build(rson);
}

void updata(int ql,int qr,int val,int l,int r,int rt,int o){
     if(ql<=y[l] && qr>=y[r]){
         aby[rt][o]+=val;
         return;
     }
     pushdown(rt,o);
     int m=(l+r)>>1;
     if(ql<=y[m]) updata(ql,qr,val,lson,o);
     if(qr>=y[m+1]) updata(ql,qr,val,rson,o);
}

int query(int yy,int l,int r,int rt){
     if(l==r) return rt;
     pushdown(rt,0);
     pushdown(rt,1);
     int m=(l+r)>>1;
     if(yy<=y[m]) return query(yy,lson);
     return query(yy,rson);
}
     
int main()
{
    int i,j,k,n;
    while(sf("%d",&n)!=EOF,n){
        for(i=0;i<n;i++){
            sf("%d%d",&p[i].x,&p[i].y);
            y[i]=p[i].y;//记录点的 y 值 
        }
        sort(p,p+n,cmp); //所有点按 x 大小排序 ； 
        sort(y,y+n); //注意：得先排序，再进行下面的记录； 
        int m=0;
        b[0]=0;
        for(i=0;i<n;i++){
            if(y[m]<y[i]){
                a[m]=n-i;
                y[++m]=y[i];
                b[m]=i;
            }
        }
        a[m]=0;
        //记录不同 y 值的点个数 m ，写在 y[]中，记录比 y[i] 点大的个数 a[i] ，比
        // y[i] 小的点的个数 b[i]； 
        
        build(0,m,1);//建立两颗线段树，分别统计一、三象限与二、四象限点的个数，
        //初始aby[0]=a[],aby[1]=b[]； 
        int ans=0,cnt=0;
        for(i=0;i<n;){ //从左往右遍历，对于 x 坐标相等的点要集体操作； 
            for(j=i;j<n;j++){ //删除直线上的点，同时对应的 [y[0],y[j]) 的 a值-1，
            // (y[j],y[m]] 的 b值-1；   
                if(p[j].y>y[0]) updata(y[0],p[j].y-1,-1,0,m,1,0);
                if(p[j].y<y[m]) updata(p[j].y+1,y[m],-1,0,m,1,1);
                if(p[j].x!=p[j+1].x) break;
            }
            int maxx=0,minn=n;
            for(k=i;k<=j;k++){
                int id=query(p[k].y,0,m,1);
                maxx=max(aby[id][1],maxx);
                minn=min(aby[id][0],minn);
            }
            if(minn==ans) q[cnt++]=maxx;
            if(minn>ans){
                ans=minn;
                q[0]=maxx;
                cnt=1;
            }
            for(;i<=j;i++){
                if(p[i].y>y[0]) updata(y[0],p[i].y-1,1,0,m,1,1);
                if(p[i].y<y[m]) updata(p[i].y+1,y[m],1,0,m,1,0);
            }//注意：加点与上面的删点是相反的，因为点要么从一象限变到二象限，要么
            //从四象限飙到三象限，所以等于点从一棵树移到了另一颗树上； 
        }
        sort(q,q+cnt);
        for(i=k=0;i<cnt;i++)
            if(q[k]<q[i]) q[++k]=q[i];
        pf("Stan: %d; Ollie:",ans);  
        for(i=0;i<=k;i++)pf(" %d",q[i]);  
        puts(";");  
    }
    return 0;
}