hdu5091 线段树

题意: 给了n个点在平面中 n<10000  然后 将这给了一个 宽为W 高为 H 的 矩形， 然后 使得这个矩形可以 涵盖最多的点有多少个，然后矩形的宽平行x 轴高平行y轴。可以将该矩形 水平或者上下移动，求他说能选中最多 多少个点，通过扫面线枚举每个x值的点 从小到大 ，选定区间后，将每个点的y值进行离散,然后以每个y为开始的点 分成 上下 的 区间 k个，然后建立一个 1到k 的 线段树， 对于每次选举的x 区间 操作这颗线段树， 因为我们知道 ， 对与 一个 y 他可能属于很多的区间， 这些区间是连续的，通过这个 我们可以用线段是的延迟更新解决这个问题。

#include <iostream>
#include <cstdio>
#include <string.h>
#include <algorithm>
#include <vector>
using namespace std;
const int maxn = 10005;
typedef int ll;
struct point{
  ll x,y;
  point(ll a=0, ll b=0){
       x=a; y=b;
  }
}P[maxn];
ll cL,cR,value,N,H,W,ynum,xnum;
struct Itree{
    ll ma[maxn*4],se[maxn*4];
    void build(ll o, ll L, ll R){
         se[o]=ma[o]=0;
         if(L==R){
             return ;
         }
         ll mid=(L+R)/2;
         build(o*2,L,mid);
         build(o*2+1,mid+1,R);
    }
    void matain(ll o){
        ma[o]=max(ma[o*2]+se[o*2],ma[o*2+1]+se[o*2+1]);
    }
    void push(ll o){
        se[o*2]+=se[o];
        se[o*2+1]+=se[o];
        se[o]=0;
    }
    void update( ll o, ll L, ll R){
         if(cL<=L&&R<=cR){
             se[o]+=value;
             return ;
         }
         ll mid =(L+R)/2;
         if(se[o]!=0){
             push(o);
         }
         if(cL<=mid) update(o*2,L,mid);
         if(cR>mid) update(o*2+1,mid+1,R);
         matain(o);
    }
}Q;
vector<ll> F[maxn];
ll X[maxn];
ll Y[maxn];
void solve(ll loc){
    ll siz=F[loc].size();
    for(ll i=0; i<siz; ++i){
        ll y = F[loc][i];
        cL = lower_bound(Y,Y+ynum,y-H)-Y+1;
        cR = lower_bound(Y,Y+ynum,y)-Y+1;
        Q.update(1,1,ynum);
    }
}
int main()
{

    while(scanf("%d",&N)==1){
         if(N<0) break;
         scanf("%d%d",&W,&H);
         for(ll i=0; i<N; ++i){
             scanf("%d%d",&X[i],&Y[i]);
             P[i]=point(X[i],Y[i]);
         }
         sort(X,X+N);
          xnum = unique(X,X+N)-X;
         for(ll i=0; i<=xnum; ++i) F[i].clear();
         sort(Y,Y+N);
         ynum = unique(Y,Y+N)-Y;
         Q.build(1,1,ynum);
         for(ll i=0; i<N; ++i){
             ll loc = lower_bound(X,X+xnum,P[i].x)-X;
             F[loc].push_back(P[i].y);
         }
         ll ans=0;
         for(ll i=0,j=0;j<xnum; ++j){
             while(X[j]-X[i]>W) {
                 value=-1;
                 solve(i);
                 i++;
             }
             value=1;
             solve(j);
             ans=max(ans,Q.ma[1]+Q.se[1]);
         }
         printf("%d
",ans);
    }
    return 0;
}