luogu2280&bzoj1218 hnoi2003激光炸弹 扫描线

http://www.elijahqi.win/archives/1097
Description

一种新型的激光炸弹，可以摧毁一个边长为R的正方形内的所有的目标。现在地图上有n(N<=10000)个目标，用整数Xi,Yi(其值在[0,5000])表示目标在地图上的位置，每个目标都有一个价值。激光炸弹的投放是通过卫星定位的，但其有一个缺点，就是其爆破范围，即那个边长为R的正方形的边必须和x，y轴平行。若目标位于爆破正方形的边上，该目标将不会被摧毁。 0

Input

输入文件的第一行为正整数n和正整数R，接下来的n行每行有3个正整数，分别表示

Output

输出文件仅有一个正整数，表示一颗炸弹最多能炸掉地图上总价值为多少的目标（结果不会超过32767）。

Sample Input

2 1
0 0 1
1 1 1
Sample Output

1
这题和poj那题real 像

我们以每个炸弹的坐标为矩形的左下建立矩形

然后我们按照他们纵坐标的高度排序并且保证删除的那个操作在添加操作之前

这样统计线段树上的最大值就可以了

画个图就知道为什么了，因为那个框里是都能框住的 所以假如两个框重合那么他们中间的部分就是我可以安防矩形的地方了

然后离散化一下的 别忘了，边界上是不可以统计的

另外没有下放的线段树真是快

#include<cstdio>
#include<algorithm>
#define N 11000
using namespace std;
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x;
} 
struct node{
    int l,r,h,c;
}line[N<<1];
struct node1{
    int l,r,left,right,sum,cnt;
}tree[N<<3];
int n,r,a,b,c,X[N<<1],num,cnt,root;
inline bool cmp(node a,node b){return a.h==b.h?a.c<b.c:a.h<b.h;}
void build(int &x,int l,int r){
    x=++cnt;tree[x].l=l;tree[x].r=r;
    if(l==r) return;
    int mid=tree[x].l+tree[x].r>>1;
    build(tree[x].left,l,mid);build(tree[x].right,mid+1,r);
}
inline void update(int x){
    int l=tree[x].left,r=tree[x].right;
    tree[x].sum=max(tree[l].sum,tree[r].sum)+tree[x].cnt;
}
void change(int x,int l,int r,int v){
    if(l<=tree[x].l&&r>=tree[x].r){tree[x].sum+=v;tree[x].cnt+=v;return;}
    int mid=tree[x].l+tree[x].r>>1;
    if (l<=mid) change(tree[x].left,l,r,v);
    if (r>mid) change(tree[x].right,l,r,v);
    update(x);
}
int main(){
    freopen("2280.in","r",stdin);
    n=read();r=read();
    for (int i=1;i<=n;++i){
        int a=read(),b=read(),c=read();
        X[++num]=a;line[num].l=a;line[num].r=a+r;line[num].h=b;line[num].c=c;
        X[++num]=a+r;line[num].l=a;line[num].r=a+r;line[num].h=b+r;line[num].c=-c;
    }
    sort(X+1,X+num+1);sort(line+1,line+num+1,cmp);
    int nn=unique(X+1,X+num+1)-X-1;build(root,1,nn); int ans=0;
    for (int i=1;i<=num;++i){
        int l=lower_bound(X+1,X+nn+1,line[i].l)-X,r=lower_bound(X+1,X+nn+1,line[i].r)-X-1;
        if (l<=r) change(root,l,r,line[i].c);
        ans=max(ans,tree[root].sum);
    }
    printf("%d",ans);
    return 0;
}