【codevs3044】矩形面积求并

这几天入了数据结构的坑，其实简单的线段树、树状数组、左偏树之类的以前就会的，但是不是很会应用，特别是例如优化DP、扩展到高维这些应用上就是渣啊，因为之前对离散化不了解，所以看了hwzer的代码好久都没懂，到最后发现原来这个线段树维护的并不是一个数轴，只是一个数组而已，一个数组有被标记的元素，然后用离散化过的扫描线法扫出来每条线之间被覆盖的矩形距离，附上原汁原味的~~我的~~hzwer的代码。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
struct data{
    double x1,x2,y;
    int flag;
}a[808];
int col[801];
double hash[201];
double sum[801];
inline bool cp(data a,data b){return a.y<b.y;}
void pushup(int size,int l,int r)
{
    if (col[size])sum[size]=hash[r+1]-hash[l];
    else if (l==r)sum[size]=0;
    else sum[size]=sum[size<<1]+sum[size<<1|1];
}
void updata(int L,int R,int flag,int l,int r,int size)
{
    if (L<=l&&R>=r)
    {
        col[size]+=flag;
        pushup(size,l,r);
        return ;
    }
    int m=(l+r)>>1;
    if (L<=m)updata(L,R,flag,l,m,size*2);
    if (R>m)updata(L,R,flag,m+1,r,size*2|1);
    pushup(size,l,r);
}
int main()
{
    int n;
    while(cin>>n)
    {
        if (n==0)break;
        double x1,y1,x2,y2;
        for (int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            a[i*2-1].x1=a[i*2].x1=x1;
            a[i*2-1].x2=a[i*2].x2=x2;
            a[i*2-1].y=y1,a[i*2].y=y2;
            a[i*2-1].flag=1;a[i*2].flag=-1;
            hash[i*2-1]=x1;hash[i*2]=x2;
        }
        sort(a+1,a+n*2+1,cp);sort(hash+1,hash+n*2+1);
        memset(col,0,sizeof(col));
        memset(sum,0,sizeof(sum));
        double ans=0;
        for (int i=1;i<=n*2;i++)
        {
            int l=lower_bound(hash+1,hash+2*n+1,a[i].x1)-hash;
            int r=lower_bound(hash+1,hash+2*n+1,a[i].x2)-hash-1;
            if (l<=r)updata(l,r,a[i].flag,1,n*2,1);
            ans+=sum[1]*(a[i+1].y-a[i].y);
        }
        printf("%.2f\n",ans);
    }
    return 0;
}