首页 > 代码库 > hdu5091 线段树

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\n",ans);    }    return 0;}
View Code

 

hdu5091 线段树