题目链接:https://cometoj.com/contest/59/problem/D?problem_id=2713
题目描述
茶颜悦色也太好喝了!鸡尾酒在长沙的各种茶颜悦色的店铺中流连忘返。他发现长沙有炒鸡多的茶颜悦色店,走两步就能遇到一家。
“方圆一公里能有十家茶颜悦色!”鸡尾酒感叹了起来。
于是他想到了一个问题:最密集的地方,能有多少家茶颜悦色的店?
鸡尾酒将长沙地图用一个二维平面表示,他统计出了每个茶颜悦色店铺的坐标。
他想知道,在一个边长为 k且底边平行于 x轴的正方形中,最多有多少家茶颜悦色。
若茶颜悦色恰好在正方形的边上,也算在正方形之中。
输入描述
输入第一行包含两个正整数 n, k() 代表茶颜悦色店的数量和正方形的边长。
接下来 n 行每行有两个整数,描述一家茶颜悦色店的坐标 保证不会出现重复的坐标。
输出描述
输出一行一个正整数表示答案。
样例输入 1
4 2
1 1
3 1
3 4
2 2
样例输出 1
3
提示
以(1,0)(3,0)(1,2)(3,2)四个点为顶点组成的边长为2的正方形,包括了3家茶颜悦色店。
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <queue>
#include <map>
#include <iostream>
#include <algorithm>
using namespace std;
#define RD(x) scanf("%d",&x)
#define RD2(x,y) scanf("%d%d",&x,&y)
#define RD3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define clr0(x) memset(x,0,sizeof(x))
#define clr1(x) memset(x,-1,sizeof(x))
#define eps 1e-9
const double pi = acos(-1.0);
typedef long long LL;
#define lson i*2,l,m
#define rson i*2+1,m+1,r
const int MAXN=200000+5;//因为点有1W个,所以扫描线2W个,不同的Y坐标最多有2W个
int cnt[MAXN*4],sum[MAXN*4];
double Y[MAXN];
struct seg
{
double l,r,h;
int d;
seg(){}
seg(double a,double b,double c,int d):l(a),r(b),h(c),d(d){}
bool operator <(const seg&b)const
{
if(h == b.h) return d>b.d;
return h<b.h;
}
}ss[MAXN];
void PushUp(int i)
{
sum[i]=max(sum[i*2],sum[i*2+1]) + cnt[i];
}
void update(int ql,int qr,int v,int i,int l,int r)
{
if(ql<=l && r<=qr)
{
cnt[i]+=v;
sum[i]+=v;
return ;
}
int m=(l+r)>>1;
if(ql<=m) update(ql,qr,v,lson);
if(m<qr) update(ql,qr,v,rson);
PushUp(i);
}
int main()
{
int n;
double w,h;
while(~RD(n))
{
if(n == -1)
break;
scanf("%lf",&w);
h=w;
w+=0.5,h+=0.5;
double x,y;
int val;
int cnt_y=0,cnt_ss=0;//记录有多少个Y值和扫描线
for(int i=1;i<=n;i++)
{
scanf("%lf%lf",&x,&y);
//x+=20000,y+=20000;
ss[cnt_ss++] = seg(y-h/2,y+h/2,x-w/2,1);
ss[cnt_ss++] = seg(y-h/2,y+h/2,x+w/2,-1);
Y[cnt_y++] = y-h/2;
Y[cnt_y++] = y+h/2;
}
sort(ss,ss+cnt_ss);
sort(Y,Y+cnt_y);
cnt_y = unique(Y,Y+cnt_y)-Y;
int ans=0;
clr0(cnt),clr0(sum);
for(int i=0;i<cnt_ss-1;i++)
{
int ql=lower_bound(Y,Y+cnt_y,ss[i].l)-Y;
int qr=lower_bound(Y,Y+cnt_y,ss[i].r)-Y-1;
if(ql<=qr) update(ql,qr,ss[i].d,1,0,cnt_y-1);
ans=max(ans,sum[1]);
}
printf("%d\n",ans);
}
}