hdu 4007 Dave 求矩形圈点最大值

#include<algorithm>
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<vector>
#include<queue>
#include<cmath>
#include<math.h>
using namespace std;
#define inf 0x3f3f3f3f
const int dir[4][2]={0,1,0,-1,1,0,-1,0};
#define maxn 2005

struct Tree
{
int l,r;//左右端点
int c;//覆盖情况
int cnt,lf,rf;//cnt为层数，lf为实际左坐标，rf为实际右坐标
}tree[maxn*4];

struct Line
{
int x,y1,y2;
int f;//直线的标记
}line[maxn];

bool cmp(Line a,Line b)//sort排序的函数
{
if(a.x!=b.x)
return a.x < b.x;
else
return a.f < b.f;
}

int y[maxn];//记录y坐标

void build(int l,int r,int root)
{
tree[root].l=l;tree[root].r=r;
tree[root].c=tree[root].cnt=0;
tree[root].lf=y[l];
tree[root].rf=y[r];
if(l+1==r)return;
int mid=(l+r)>>1;
build(l,mid,root<<1);
build(mid,r,root<<1|1);
}

void calen(int root)
{
if(tree[root].l+1==tree[root].r) tree[root].cnt=tree[root].c;
else tree[root].cnt=max(tree[root<<1].cnt,tree[root<<1|1].cnt)+tree[root].c;
}

void update(int root,struct Line e)
{
if(e.y1==tree[root].lf&&e.y2==tree[root].rf)
{
tree[root].c+=e.f;
calen(root);
return;
}
if(e.y2<=tree[root<<1].rf)update(root<<1,e);
else if(e.y1>=tree[root<<1|1].lf)update(root<<1|1,e);
else
{
Line tmp=e;
tmp.y2=tree[root<<1].rf;
update(root<<1,tmp);
tmp=e;
tmp.y1=tree[root<<1|1].lf;
update(root<<1|1,tmp);
}
calen(root);
}

int main()
{
int i,j;
int n,r,k;
int x1,y1,x2,y2;
while(scanf("%d%d",&n,&r)!=EOF)
{
k=1;
for(i=1;i<=n;i++)
{
scanf("%d%d",&x1,&y1);
x2=x1+r+1;
y2=y1+r+1;
line[k].x=x1;
line[k].y1=y1;
line[k].y2=y2;
line[k].f=1;
y[k]=y1;
k++;
line[k].x=x2;
line[k].y1=y1;
line[k].y2=y2;
line[k].f=-1;
y[k]=y2;
k++;
}
k--;
sort(line+1,line+k+1,cmp);
sort(y+1,y+k+1);
int mm;
mm=unique(y+1,y+1+k)-y-1;
build(1,mm,1);
int ans=0;
for(i=1;i<=k;i++)
{
update(1,line[i]);
if(tree[1].cnt>ans)
ans=tree[1].cnt;
}
printf("%d\n",ans);
}
return 0;
}