HDU 4007 Dave（11年大连 线段树+离散化+扫描线）

转载请注明出处，谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents           by---cxlove 

题目：给出一些点，给出一个连长为R的正方形，问最多能覆盖多少个点 

http://acm.hdu.edu.cn/showproblem.php?pid=4007

记忆中这题在当时是比较坑的，题意不清。

正方形的边与坐标轴是平行的，这样就水多了。范围也不是很大

对于每一个点，以他为左下角，搞出一个正方形出来，然后对于所有的正方形求一次矩形交，找出覆盖次数最多的。

竟然对于扫描线不是很熟，囧。对于上下两条边，建立一个标号，很方便

#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<algorithm>
#include<set>
#define inf 110000
#define N 1005
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define pb(a) push_back(a)
#define mem(a,b) memset(a,b,sizeof(a))
#define eps 1e-9
#define zero(a) fabs(a)<eps
#define LL long long
#define lson (step<<1)
#define rson (step<<1|1)
#define MOD 1000000007
using namespace std;
struct Node{
	int left,right,mid,cover,val;
}L[N*2*4];
struct Seg{
	int y,x1,x2;
	int flag;
	Seg(){}
	Seg(int _y,int _x1,int _x2,int _flag):y(_y),x1(_x1),x2(_x2),flag(_flag){}
	bool operator<(const Seg s)const{
		return s.y==y?s.flag<flag:s.y>y;
	}
}seg[N*2];
int n,r,x[N*2],m;
int BinSearch(int num){
	int low=0,high=m-1,mid;
	while(low<=high){
		mid=(low+high)/2;
		if(x[mid]==num) return mid;
		if(x[mid]<num) low=mid+1;
		else high=mid-1;
	}
}
void Bulid(int step,int l,int r){
	L[step].left=l;
	L[step].right=r;
	L[step].mid=(l+r)/2;
	L[step].cover=L[step].val=0;
	if(l==r) return;
	Bulid(lson,l,L[step].mid);
	Bulid(rson,L[step].mid+1,r);
}
void Push_Down(int step){
	if(L[step].cover){
		L[lson].cover+=L[step].cover;
		L[rson].cover+=L[step].cover;
		L[lson].val+=L[step].cover;
		L[rson].val+=L[step].cover;
		L[step].cover=0;
	}
}
void Push_Up(int step){
	L[step].val=max(L[lson].val,L[rson].val);
}
void Update(int step,int l,int r,int c){
	if(L[step].left==l&&L[step].right==r){
		L[step].cover+=c;
		L[step].val+=c;
		return;
	}
	Push_Down(step);
	if(r<=L[step].mid) Update(lson,l,r,c);
	else if(l>L[step].mid) Update(rson,l,r,c);
	else{
		Update(lson,l,L[step].mid,c);
		Update(rson,L[step].mid+1,r,c);
	}
	Push_Up(step);
}
int main(){
	while(scanf("%d%d",&n,&r)!=EOF){
		for(int i=0;i<n;i++){
			int X,Y;
			scanf("%d%d",&X,&Y);
			x[i*2]=X;x[i*2+1]=X+r+1;
			seg[i*2]=Seg(Y,X,X+r+1,1);
			seg[i*2+1]=Seg(Y+r+1,X,X+r+1,-1);
		}
		sort(x,x+2*n);
		sort(seg,seg+2*n);
		m=1;
		for(int i=1;i<2*n;i++)
			if(x[i]!=x[m-1])
				x[m++]=x[i];
		Bulid(1,0,m-1);
		int ans=0;
		for(int i=0;i<2*n;i++){
			int left=BinSearch(seg[i].x1),right=BinSearch(seg[i].x2);
			Update(1,left,right,seg[i].flag);
			if(L[1].val>ans) ans=L[1].val;
		}
		printf("%d\n",ans);
	}
	return 0;
}