bzoj 1218: [HNOI2003]激光炸弹 （扫描线+线段树）

题目描述

传送门

题目大意：平面内有n个点，求用长度为R的矩形可以框住的点的最大价值，边界上的不算。

题解

如果直接算边长为R的矩形框住的点，不容易计算。
我们把问题做变形，对于每个点我们把其扩大成一个矩形，表示只要矩形的左上角在该区域中就可以得到这个点的价值，那么我们其实就是要找到被覆盖的价值最大的位置。
根据矩形的左右边界把每个矩形拆成两个操作，在左边界+val,在右边界+1 -val。
然后按照纵坐标从小到大的顺序扫描，每扫到一列就把这一列之前的操作都加入，然后用线段树中的最大值更新答案。
因为边界上的不算，所以我们在每行每列中间再加入新的行列。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 20003
using namespace std;
struct data{
    int l,r,x,val;
}a[N];
int n,R,cnt,tr[N*4],delta[N*4];
void update(int now)
{
    tr[now]=max(tr[now<<1],tr[now<<1|1]);
}
void change(int now,int val)
{
    tr[now]+=val; delta[now]+=val;
}
void pushdown(int now)
{
    if (delta[now]) {
        change(now<<1,delta[now]);
        change(now<<1|1,delta[now]);
        delta[now]=0;
    }
}
void query(int now,int l,int r,int ll,int rr,int val)
{
    if (ll<=l&&r<=rr) {
        change(now,val);
        return;
    }
    int mid=(l+r)/2;
    pushdown(now);
    if (ll<=mid) query(now<<1,l,mid,ll,rr,val);
    if (rr>mid) query(now<<1|1,mid+1,r,ll,rr,val);
    update(now);
}
int cmp(data a,data b)
{
    return a.x<b.x;
}
int main()
{
   freopen("a.in","r",stdin);
   freopen("my.out","w",stdout);
   scanf("%d%d",&n,&R); int cnt=0; int mx=0;
   for (int i=1;i<=n;i++) {
     int x,y,v; 
     scanf("%d%d%d",&x,&y,&v); x++; y++; mx=max(mx,y+1); mx=max(mx,x+1);
     ++cnt; a[cnt].l=max((x-R)*2+1,1); a[cnt].r=max(x*2-1,1); a[cnt].x=max((y-R)*2+1,1); a[cnt].val=v;
     ++cnt; a[cnt].l=a[cnt-1].l; a[cnt].r=a[cnt-1].r; a[cnt].x=y*2; a[cnt].val=-v;
   }    
   sort(a+1,a+cnt+1,cmp);
   int l=1; n=mx*2; int ans=0;
   for (int i=1;i<=n;i++){
        while (a[l].x<=i&&l<=cnt) {
             query(1,1,n,a[l].l,a[l].r,a[l].val);
             l++;
        }
        ans=max(ans,tr[1]);
   }
   printf("%d\n",ans);
}