poj2482Stars in Your Window（线段树+扫描线+离散化）

题目链接：http://poj.org/problem?id=2482

链接题目大意：给出n个星星的坐标，每个星星有一个亮度，给出一个矩形的长和宽，问矩形能包括的星星的最大亮度和（不包括边框）。

思路：假设每个星星都是矩形的最左下点，那么每一个星星都可以得到一个矩形，(x,y)->(x,y,x+w,y+h)

这里的亮度相当于下边和上边的值不是1和-1而是c和-c

左右两边即使出现在边界处，也不会将其记录，所以我们只要考虑上下边！！即扫描时应先扫描值为负的高线，意思就是结构体排序时，高度相同时，上边在前（负的），下边在后（正的），这与求面积交时略有区别！！

这题相当于在求底边最多能重叠几层！！！

解释下pushup函数

cnt数组意思相当于当前区间长度的边的条数

sum数组是当前区间最多叠了几层

那么pushup函数的意思 就是对于左右两个子区间中最对叠了几层（这个层数不包括当前区间长的边的条数），所以还要加上cnt【rt】当前区间的边的条数！！

#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
const int maxn = 50000;
long long  cnt[maxn << 2];
long long  sum[maxn << 2];
long long   X[maxn];
struct Seg
{
   long long  h , l , r;
    long long  s;
    Seg(){}
    Seg(long long   a,long long   b,long long  c,long long  d) : l(a) , r(b) , h(c) , s(d) {}
    bool operator < (const Seg &cmp) const
    {
         if( h==cmp.h)return s<cmp.s;
        return h < cmp.h;
    }
}ss[maxn];
void PushUp(int rt)
 {
    sum[rt]=max(sum[rt<<1],sum[rt<<1|1])+cnt[rt];
}
void update(long long  L,long long  R,long long  c,long long  l,long long  r,long long  rt)
{
    if (L <= l && r <= R)
    {
        cnt[rt] += c;
        sum[rt]+=c;
        return ;
    }

    long long  m = (l + r) >> 1;
    if (L <= m)
        update(L , R , c , lson);
    if (m < R)
        update(L , R , c , rson);
    PushUp(rt);
}
int main()
{
   long long  n,w,h;
    while (~scanf("%lld%lld%lld",&n,&w,&h))
    {
       long long  m = 0;
        while (n --)
        {
            long long  a , b , c;
            scanf("%lld%lld%lld",&a,&b,&c);
            X[m] = a;
            ss[m++] = Seg(a , a+w , b , c);
            X[m] = a+w;
            ss[m++] = Seg(a ,a+w , b+h , -c);
        }
        sort(X , X + m);
        sort(ss , ss + m);
         int k=unique(X,X+m)-X;

        memset(cnt , 0 , sizeof(cnt));
        memset(sum , 0 , sizeof(sum));
        long long   ret = 0;
        for (int i = 0 ; i < m - 1 ; i ++)
        {
            long long  l=lower_bound(X,X+k,ss[i].l)-X;
            long long  r =lower_bound(X,X+k,ss[i].r)-X-1;
            if (l <= r) update(l , r , ss[i].s , 0 , k-1 , 1);//扫描线段时更新底边长度和底边相差个数
                ret= max(ret,sum[1]) ;//底乘以高
        }
        printf("%lld\n", ret);
    }
    return 0;
}