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 = ;
typedef int ll;
struct point{
  ll x,y;
  point(ll a=, ll b=){
       x=a; y=b;
  }
}P[maxn];
ll cL,cR,value,N,H,W,ynum,xnum;
struct Itree{
    ll ma[maxn*],se[maxn*];
    void build(ll o, ll L, ll R){
         se[o]=ma[o]=;
         if(L==R){
             return ;
         }
         ll mid=(L+R)/;
         build(o*,L,mid);
         build(o*+,mid+,R);
    }
    void matain(ll o){
        ma[o]=max(ma[o*]+se[o*],ma[o*+]+se[o*+]);
    }
    void push(ll o){
        se[o*]+=se[o];
        se[o*+]+=se[o];
        se[o]=;
    }
    void update( ll o, ll L, ll R){
         if(cL<=L&&R<=cR){
             se[o]+=value;
             return ;
         }
         ll mid =(L+R)/;
         if(se[o]!=){
             push(o);
         }
         if(cL<=mid) update(o*,L,mid);
         if(cR>mid) update(o*+,mid+,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=; i<siz; ++i){
        ll y = F[loc][i];
        cL = lower_bound(Y,Y+ynum,y-H)-Y+;
        cR = lower_bound(Y,Y+ynum,y)-Y+;
        Q.update(,,ynum);
    }
}
int main()
{

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