codevs 3044 矩形面积求并

3044 矩形面积求并

 

题目描述 Description

输入n个矩形，求他们总共占地面积（也就是求一下面积的并）

输入描述 Input Description

可能有多组数据，读到n=0为止(不超过15组)

每组数据第一行一个数n，表示矩形个数(n<=100)

接下来n行每行4个实数x1,y1,x2,y1(0 <= x1 < x2 <= 100000;0 <= y1 < y2 <= 100000)，表示矩形的左下角坐标和右上角坐标

输出描述 Output Description

每组数据输出一行表示答案

样例输入 Sample Input

2
10 10 20 20
15 15 25 25.5
0

样例输出 Sample Output

180.00

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 210
using namespace std;
struct node{
    double l,r,h;int d;
    node(){}
    node(double l,double r,double h,int d):l(l),r(r),h(h),d(d){}
    bool operator < (const node &b)const{
        return h<b.h;
    }
}a[maxn];
int cnt[maxn*],n;
double sum[maxn*],all[maxn];
void pushup(int l,int r,int k){
    if(cnt[k])sum[k]=all[r+]-all[l];
    else if(l==r)sum[k]=;
    else sum[k]=sum[k<<]+sum[k<<|];
}
void update(int opl,int opr,int opv,int l,int r,int k){
    if(l>=opl&&r<=opr){
        cnt[k]+=opv;
        pushup(l,r,k);
        return;
    }
    int mid=(l+r)>>;
    if(opl<=mid)update(opl,opr,opv,l,mid,k<<);
    if(opr>mid)update(opl,opr,opv,mid+,r,k<<|);
    pushup(l,r,k);
}
int main(){
    while(){
        scanf("%d",&n);
        if(n==)return ;
        double x1,x2,y1,y2;
        for(int i=;i<=n;i++){
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            a[i]=node(x1,x2,y1,);
            a[i+n]=node(x1,x2,y2,-);
            all[i]=x1;all[i+n]=x2;
        }
        n<<=;
        sort(a+,a+n+);
        sort(all+,all+n+);
        int m=unique(all+,all+n+)-all-;
        memset(cnt,,sizeof(cnt));
        memset(sum,,sizeof(sum));
        double ans=;
        for(int i=;i<n;i++){
            int l=lower_bound(all+,all+m+,a[i].l)-all;
            int r=lower_bound(all+,all+m+,a[i].r)-all;
            if(l<r)update(l,r-,a[i].d,,m,);
            ans+=sum[]*(a[i+].h-a[i].h);
        }
        printf("%.2lf\n",ans);
    }
    return ;
}