POJ 1389 Area of Simple Polygons 线段树求矩形面积

题意：给一些矩形，让求这些矩形的面积。矩形是给出的左下角和右下角的坐标，都是正整数。

思路：简单的线段树题目，就是求矩形面积。因为都是整数，而且数据范围不大，因此不需要离散化。

代码：

 

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string.h>
using namespace std;

const int N = 50010,M = 1010;
struct tree{
    int lp,rp,len,num;
    int getmid(){
       return (lp + rp) / 2;
    }
}tt[N * 4];
struct line{
    int lp,rp,cnt,value;
}LL[M * 2];
void built_tree(int lpos,int rpos,int pos){
    tt[pos].lp = lpos;
    tt[pos].rp = rpos;
    tt[pos].len = tt[pos].num = 0;
    if(tt[pos].lp + 1 == tt[pos].rp)
        return;
    int mid = tt[pos].getmid();
    built_tree(lpos,mid,pos*2);
    built_tree(mid,rpos,pos*2+1);
}
bool cmp(line a,line b){
    return a.value > b.value;
}
void getlen(int pos){
    if(tt[pos].num > 0)
        tt[pos].len = tt[pos].rp - tt[pos].lp;
    else if(tt[pos].lp + 1 == tt[pos].rp)
        tt[pos].len = 0;
    else
        tt[pos].len = tt[pos * 2].len + tt[pos * 2 + 1].len;
}
void update(int pos,int id){
    if(tt[pos].lp >= LL[id].lp && tt[pos].rp <= LL[id].rp){
        tt[pos].num += LL[id].cnt;
        getlen(pos);
        return;
    }
    if(tt[pos].lp + 1 == tt[pos].rp)
        return;
    int mid = tt[pos].getmid();
    if(LL[id].lp >= mid){
       update(pos * 2 + 1,id);
    }
    else if(LL[id].rp <= mid){
       update(pos * 2,id);
    }
    else{
       update(pos * 2,id);
       update(pos * 2 + 1,id);
    }
    getlen(pos);
}
int main(){
    //freopen("1.txt","r",stdin);
    int a,b,c,d,tot = 0;
    while(scanf("%d%d%d%d",&a,&b,&c,&d) && (a + b + c + d) != -4){
        tot = 0;
        LL[tot].lp = a; LL[tot].rp = c; LL[tot].value = d; LL[tot++].cnt = 1;
        LL[tot].lp = a; LL[tot].rp = c; LL[tot].value = b; LL[tot++].cnt = -1;
        int x1,y1,x2,y2;
        while(1){
           scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
           if(x1 + y1 + x2 + y2 == -4)
               break;
           LL[tot].lp = x1; LL[tot].rp = x2; LL[tot].value = y2; LL[tot++].cnt = 1;
           LL[tot].lp = x1; LL[tot].rp = x2; LL[tot].value = y1; LL[tot++].cnt = -1;
        }
        built_tree(0,N,1);
        sort(LL,LL+tot,cmp);
        update(1,0);
        long long ans = 0;
        for(int i = 1; i < tot; ++i){
           ans += tt[1].len * (LL[i-1].value - LL[i].value);
           update(1,i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}