poj2482 Stars in Your Window(线段树+扫描线+离散化)

这题目好浪漫啊~~~~~~~~~~~~~
说正事：对于每个点，我们把它变成一个范围，这个范围的意义是，只要矩形的右上角在这个范围内，这个点就可以被这个矩形覆盖。则这个范围显然就是(x,y,x+w-1,y+h-1)的一个矩形。现在我们的题目就变为：平面上有n个矩形，每个矩形内的点都有一个权值，可以互相叠加，问你权值最大的点是多少？我们可以用扫描线解决这个问题，按x排序，每次入边就是对区间[y1,y2]做区间加法，出边就是对区间[y1,y2]做区间减法，我们要找的就是整个区间最大的数，用线段树来维护即可。

#include <cstdio>#include <cstring>#include <algorithm>using namespace std;#define ll long long#define inf 0x3f3f3f3f#define N 20010inline int read(){    int x=0,f=1;char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();    return x*f;}int n,w,h,aa[N],num,m,ans;struct Lines{    unsigned int x,y1,y2;int val;}a[N];struct node{    int mx,add;}tree[N<<2];inline bool cmp(Lines x,Lines y){    return x.x==y.x?x.val>y.val:x.x<y.x;}inline void pushup(int p){    tree[p].mx=max(tree[p<<1].mx,tree[p<<1|1].mx);}inline void pushdown(int p){    if(!tree[p].add) return;    tree[p<<1].add+=tree[p].add;tree[p<<1|1].add+=tree[p].add;    tree[p<<1].mx+=tree[p].add;tree[p<<1|1].mx+=tree[p].add;    tree[p].add=0;}void build(int p,int l,int r){    tree[p].add=tree[p].mx=0;    if(l==r) return;    int mid=l+r>>1;    build(p<<1,l,mid);build(p<<1|1,mid+1,r);}void change(int p,int l,int r,int x,int y,int val){    if(x<=l&&r<=y){tree[p].add+=val;tree[p].mx+=val;return;}    int mid=l+r>>1;pushdown(p);    if(x<=mid) change(p<<1,l,mid,x,y,val);    if(y>mid) change(p<<1|1,mid+1,r,x,y,val);    pushup(p);}int main(){//  freopen("a.in","r",stdin);    while(~scanf("%d",&n)){        w=read();h=read();num=0;m=0;ans=0;        while(n--){            unsigned int x=read(),y=read(),val=read();            a[++num].x=x;a[num].y1=y;a[num].y2=y+h-1;a[num].val=val;            a[++num].x=x+w-1;a[num].y1=y;a[num].y2=y+h-1;a[num].val=-val;            aa[++m]=y;aa[++m]=y+h-1;        }sort(aa+1,aa+m+1);m=unique(aa+1,aa+m+1)-aa-1;        sort(a+1,a+num+1,cmp);build(1,1,m);        for(int i=1;i<=num;++i){            change(1,1,m,lower_bound(aa+1,aa+m+1,a[i].y1)-aa,lower_bound(aa+1,aa+m+1,a[i].y2)-aa,a[i].val);            ans=max(ans,tree[1].mx);        }        printf("%d\n",ans);    }    return 0;}