题目链接:Beam Cannon
解题思路:之前只是听说过扫描线没有实际用过,想了好久,只能看别人博客找找解答,不知道是他们看别人代码后写的博客自己说不清楚,还是懒得说,半天看不懂。后来自己想了想,总结出了比较好理解的
扫描线从左向右扫,先将即x小的先被扫到,x大的后被扫到。扫描线垂直于x轴,平行于y轴。扫描线对于矩形来说很好用,这个题就是矩形。
我们以一个点为矩形左上角的点,扫出与题目矩形大小相同的矩形出来。
矩形如下
(x, y) (x + w, y)
(x, y - h) (x + w, y - h)
这就是我们扫出的矩形,那么这个面积又有什么含义呢?这个面积就代表着当以这个区域中任意一个点为矩形右下角的点,画出宽为w高为h的矩形,都能包含(x, y)这个点。
当两个不同的点画出矩形有重合的话,那就说明这个区域中任意一个点为右下角的点,画出题示矩形,都能包含这两个点。那么这就成了求区域覆盖数的最大值,就是很简单的扫描线了,最后再用线段树统计。
#include<cstdio>
#include<vector>
#include<string>
#include<algorithm>
#include<iostream>
#define MAX 60010
#define C 20001
#define lson l, m, rt << 1
#define rson m + 1, r, rt << 1 | 1
using namespace std;
struct point{
int x, y, v;
point(){}
point(int xx, int yy, int vv){
x = xx, y = yy, v = vv;
}
bool operator < (const point &a)const{
return x < a.x || x == a.x && v > a.v;
}
};
vector<point> p;
int v[MAX << 2], add[MAX << 2];
int n, w, h;
void pushUp(int rt){
v[rt] = max(v[rt << 1], v[rt << 1 | 1]);
}
void pushDown(int rt){
if(add[rt]){
add[rt << 1] += add[rt];
add[rt << 1 | 1] += add[rt];
v[rt << 1 | 1] += add[rt];
v[rt << 1] += add[rt];
add[rt] = 0;
}
}
void build(int l, int r, int rt){
if(l == r){
add[rt] = v[rt] = 0;
return;
}
int m = (l + r) >> 1;
build(lson);
build(rson);
pushUp(rt);
}
void modify(int a, int b, int ad, int l, int r, int rt){
if(a <= l && r <= b){
add[rt] += ad;
v[rt] += ad;
return;
}
int m = (l + r) >> 1;
pushDown(rt);
if(a <= m){
modify(a, b, ad, lson);
}
if(m < b){
modify(a, b, ad, rson);
}
pushUp(rt);
}
int querry(int a, int b, int l, int r, int rt){
if(a <= l && r <= b){
return v[rt];
}
pushDown(rt);
int m = (l + r) >> 1;
int ret = 0;
if(a <= m){
ret = max(ret, querry(a, b, lson));
}
if(m < b){
ret = max(ret, querry(a, b, rson));
}
return ret;
}
int main(){
int x, y;
int i, j, k, ans;
while(~scanf("%d", &n)){
if(n == -1) break;
scanf("%d%d", &w, &h);
ans = 0;
p.clear();
for(i = 0; i < n; i++){
scanf("%d%d", &x, &y);
x += C;
y += C;
p.push_back(point(x, y, 1));
p.push_back(point(x + w, y, -1));
}
sort(p.begin(), p.end());
memset(v, 0, sizeof(v));
memset(add, 0, sizeof(add));
for(i = 0; i < p.size(); i++){
modify(p[i].y, p[i].y + h, p[i].v, 1, 60001, 1);
int tem = querry(1, 60001, 1, 60001, 1);
ans = max(ans, tem);
}
printf("%d\n", ans);
}
return 0;
}