标签：style blog http io color ar os for sp
题意: 给了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; }
标签：style blog http io color ar os for sp
原文地址：http://www.cnblogs.com/Opaser/p/4069823.html
