codves 3044 矩形面积求并

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,opl,opr,flag;
int col[801];
double hash[201],sum[1601],ans;
struct node
{
    double x1,x2,h;
    int f;
}e[202];
bool cmp(node a,node b) {return a.h<b.h;}
inline void up(int k,int l,int r)
{
    if(col[k]) sum[k]=hash[r+1]-hash[l];
    else sum[k]=sum[k<<1]+sum[(k<<1)+1];
}
inline void change(int k,int l,int r)
{
    if(l>=opl&&r<=opr)
    {
        col[k]+=flag; 
        up(k,l,r);
        return;
    }
    int m=l+r>>1;
    if(opl<=m) change(k<<1,l,m);
    if(opr>m) change((k<<1)+1,m+1,r);
    up(k,l,r);
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        if(!n) return 0;
        double x1,y1,x2,y2;
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            e[i*2-1].x1=e[i*2].x1=x1;
            e[i*2-1].x2=e[i*2].x2=x2;
            e[i*2-1].h=y1;e[i*2].h=y2;
            hash[i*2-1]=x1;hash[i*2]=x2;
            e[i*2-1].f=1;e[i*2].f=-1;
        }
        sort(hash+1,hash+2*n+1);
        sort(e+1,e+2*n+1,cmp);
        ans=0;
        for(int i=1;i<=2*n;i++)
        {
            opl=lower_bound(hash+1,hash+2*n+1,e[i].x1)-hash;
            opr=lower_bound(hash+1,hash+2*n+1,e[i].x2)-hash-1;
            flag=e[i].f;
            change(1,1,2*n);
            ans+=sum[1]*(e[i+1].h-e[i].h);
        }
        printf("%.2lf\n",ans);
    }
}