HDU 5091 Beam Cannon(扫描线 + 线段树)

题目链接：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;
}