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

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5091

思路：扫描线的陈题，先线段树的功能，区间价值，维护的是区间上的最大值，然后将每个敌方单位分解成两个事件，不断更新就行。

code：

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int maxn=70050;
const int maxt=30500;
const int mo=20005;

struct sg_tree
{
    int l,r;
    int col,rmax;
} ss[maxn<<2];

struct Thing
{
    int l,r,y,kind;
} P[maxt];

int n;
int nn,w,h,ta;

void Push_up(int k)
{
    ss[k].rmax=max(ss[(k<<1)+1].rmax,ss[(k<<1)+2].rmax);
}
void Push_down(int k)
{
    ss[(k<<1)+1].col+=ss[k].col;
    ss[(k<<1)+2].col+=ss[k].col;
    ss[(k<<1)+1].rmax+=ss[k].col;
    ss[(k<<1)+2].rmax+=ss[k].col;
    ss[k].col=0;
}
void Build(int k,int l,int r)
{
    ss[k].col=0;
    ss[k].rmax=0;
    ss[k].l=l;
    ss[k].r=r;
    if(r==l) return ;
    Build((k<<1)+1,l,(l+r)/2);
    Build((k<<1)+2,(l+r)/2+1,r);
}
void init(int nn)
{
    n=1;
    while(n<nn) n*=2;
    Build(0,0,n-1);
}
void update(int a,int b,int k,int cc)
{
    if(b<ss[k].l||a>ss[k].r) return ;
    if(a<=ss[k].l&&ss[k].r<=b){
        ss[k].col+=cc;
        ss[k].rmax+=cc;
    }
    else{
        int mid=(ss[k].l+ss[k].r)/2;
        Push_down(k);
        update(a,b,(k<<1)+1,cc);
        update(a,b,(k<<1)+2,cc);
        Push_up(k);
    }
}
int query(int a,int b,int k)
{
    if(a<=ss[k].l&&ss[k].r<=b) return ss[k].rmax;
    else{
        int v1,v2;
        v1=v2=0;
        int mid=(ss[k].l+ss[k].r)/2;
        if(a<mid) v1=query(a,b,(k<<1)+1);
        if(b>mid) v2=query(a,b,(k<<1)+2);
        return max(v1,v2);
    }
}

bool cmp(Thing A,Thing B)
{
    if(A.y==B.y) return A.kind<B.kind;
    return A.y<B.y;
}

//最佳的放置状态下 矩阵的一边肯定是和每个节点重合
//我们假设左边是和敌人节点重叠的
//那么对每个敌方单位他能被包含的区域就是x+w
//对于每个敌方单位我们构造两种事件 l,r都为x和x+w y为y和y+h+1 kind分别为1和-1
//然后把所有事件按y的坐标排序 一次进行执行即可
//每次执行后都用根节点的值来更新res
void creat(int x,int y)
{
    x+=mo; y+=mo;
    P[ta].l=x; P[ta].y=y;
    if(x+w<=50000) P[ta].r=x+w;
    else P[ta].r=50000;
    P[ta++].kind=1;

    P[ta].l=P[ta-1].l;
    P[ta].r=P[ta-1].r;
    if(y+h+1<=50000) P[ta].y=y+h+1;
    else P[ta].y=50000;
    P[ta++].kind=-1;
}

int main()
{
    int x,y,res;
    while(scanf("%d",&nn),nn>=0){
        scanf("%d%d",&w,&h);
        res=ta=0;
        for(int i=0;i<nn;i++){
            scanf("%d%d",&x,&y);
            creat(x,y);
        }
        init(50005);
        sort(P,P+ta,cmp);
        //for(int i=0;i<ta;i++){
           // printf("l=%d r=%d y=%d kind=%d\n",P[i].l,P[i].r,P[i].y,P[i].kind);
       // }
        for(int i=0;i<ta;i++){
            update(P[i].l,P[i].r,0,P[i].kind);
            res=max(res,ss[0].rmax);
        }
        printf("%d\n",res);
    }
    return 0;
}