POJ_2482_Stars in Your Window_线段树、离散化处理

手机屏摔碎了哦凑。

题意

一个平面里有N(N<10000)个点，每个点 有一个值Vi，用一个长宽为W，H的矩形框，在平面内平移，不能旋转，最多能使框内（不包括边界）的点的V值和为多少 ？

IO

Input
There are several test cases in the input. The first line of each case contains 3 integers: n, W, H, indicating the number of stars, the horizontal length and the vertical height of the rectangle-shaped window. Then n lines follow, with 3 integers each: x, y, c, telling the location (x, y) and the brightness of each star. No two stars are on the same point.

There are at least 1 and at most 10000 stars in the sky. 1<=W，H<=1000000, 0<=x，y<2^31.

Output
For each test case, output the maximum brightness in a single line.

分析

首先要对问题进行一次转化，对每个点来说，矩形框中心的位置必须在以当前点为中心的一个w,h的矩形框内，那么我们就可以把点转化为矩形，而把矩形转化为点。
之所以要进行这样的转化，是因为我们要移动矩形，而点是不动的，显然矩形这一模型比点要复杂，因此我们把静态的点转化为更复杂的矩形，而把动态的矩形转化为简单的点。
用一条扫描线的从左往右扫，维护沿y轴方向矩形数量的分布信息，注意对于转化后的问题，扫描线进入一个矩形以后，分布信息中就应在对应位置加上它，知道出这个矩形，删除对应信息。用线段树来维护当前沿y轴矩形数量分布信息。
实现会遇到一些小问题：
从左到右的距离太大——排序矩形的左右边，扫描矩形左右边。
纵向距离太大——离散化纵向坐标。
按照题意点在框上的情况不算——由于输入都是整数，用一个小于0.5的小数修正即可。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define MXN 10010
int ll[MXN<<4],rr[MXN<<4],sum[MXN<<4],flag[MXN<<4];
int n,w,h;
int ans;
void build(int id,int l,int r){
    ll[id]=l,rr[id]=r,sum[id]=flag[id]=0;
    int m=(l+r)>>1,ls=id<<1,rs=ls|1;
    if(l!=r){
        build(ls,l,m);
        build(rs,m+1,r);
    }
}
struct line{
    double x,l,h,v;
    bool end;
    line(){}
    line(double a,double b,double c,long long e,bool d){x=a,l=b,h=c,end=d,v=e;}
    bool operator < (const line& rgt)const{
        return x<rgt.x||x==rgt.x&&end;
    }
}l[MXN<<2];
long long lcnt,tot;
double dy[MXN<<2];
void input(){
    ans=lcnt=tot=0;
    for(int i=0;i<n;++i){
        int tx,ty,tc;
        scanf("%d%d%d",&tx,&ty,&tc);
        double tl=ty-h/2.0+0.1;
        double th=ty+h/2.0-0.1;
        double ttx=tx-w/2.0+0.1;
        l[lcnt++]=line(ttx,tl,th,tc,false);
        ttx=tx+w/2.0-0.1;
        l[lcnt++]=line(ttx,tl,th,tc,true);
        dy[tot++]=tl;
        dy[tot++]=th;
    }
    sort(l,l+lcnt);
    sort(dy,dy+tot);
    tot=unique(dy,dy+tot)-dy;
}
void down(int id);
void add(int id,double l,double r,int v,bool end){
    if(flag[id])    down(id);
    if(l==dy[ll[id]]&&r==dy[rr[id]]){
        flag[id]+=end ? -v : v;
        sum[id]+=flag[id];
        return;
    }
    int m=(ll[id]+rr[id])>>1,ls=id<<1,rs=ls|1;
    if(r<=dy[m])    add(ls,l,r,v,end);
    else if(l>dy[m])    add(rs,l,r,v,end);
    else{
        add(ls,l,dy[m],v,end);
        add(rs,dy[m+1],r,v,end);
    }
    sum[id]=max(sum[ls],sum[rs]);
}
void down(int id){
    if(ll[id]==rr[id]){
        flag[id]=0;
        return;
    }
    int m=(ll[id]+rr[id])>>1,ls=id<<1,rs=ls|1;
    add(ls,dy[ll[id]],dy[m],abs(flag[id]),flag[id]<0);
    add(rs,dy[m+1],dy[rr[id]],abs(flag[id]),flag[id]<0);
    flag[id]=0;
}
int query(int id){
    if(flag[id]){
        down(id);
    }
    return sum[id];
}
void solve(){
    for(int i=0;i<lcnt;++i){
        add(1,l[i].l,l[i].h,l[i].v,l[i].end);
        ans=max(query(1),ans);
    }
}
int main(){
    //freopen("data.txt","r",stdin);
    //freopen("ans.txt","w",stdout);
    while(scanf("%d%d%d",&n,&w,&h)!=EOF){
        input();
        build(1,0,tot-1);
        solve();
        printf("%d\n",ans);
    }
    return 0;
}