矩形面积交

例题：hdu1255
题意：裸题，不过是要交两次以上的面积
基本上都和面积并差不多

#include<bits/stdc++.h>
using namespace std;
#define lson rt<<1
#define rson rt<<1|1

const int maxn=2005;
struct node//和面积并一样，存放横向的线（把一个矩形拆成上下两根线）
{
    double x1,x2,h;
    int sta;
    node(){}
    node(double x1,double x2,double h,int sta):x1(x1),x2(x2),h(h),sta(sta){}
    bool operator < (const node &a)const//这些横向线按y升序排序
    {
        if(h==a.h)return x1<a.x1;
        return h<a.h;
    }
}p[maxn<<1];

double x[maxn];//存放所有出现过的x值，用来离散化
double tree1[maxn<<2],tree2[maxn<<2];//分别维护只有一次和两次的x区间长度
int lazy[maxn<<2];//对x区间的标记

void pushup(int rt,int l,int r)
{
    if(lazy[rt]>=2)tree1[rt]=tree2[rt]=x[r+1]-x[l];//如果当前这个区间已经有两次了，一次和二次直接加上这段长度
    else if(lazy[rt]==1)//如果当前这个区间只有一次，一次的要更新
    {
        tree1[rt]=x[r+1]-x[l];
        if(l==r)tree2[rt]=0;
        else tree2[rt]=tree1[lson]+tree1[rson];
    }
    else
    {
        if(l==r)tree1[rt]=tree2[rt]=0;
        else tree1[rt]=tree1[lson]+tree1[rson],tree2[rt]=tree2[lson]+tree2[rson];
    }
}

void updata(int rt,int l,int r,int x,int y,int v)
{
    if(l==r)
    {
        lazy[rt]+=v;
        pushup(rt,l,r);
        return;
    }
    int mid=(l+r)>>1;
    if(x<=mid)updata(lson,l,mid,x,y,v);
    if(y>=mid+1)updata(rson,mid+1,r,x,y,v);
    pushup(rt,l,r);
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);
        double x1,x2,y1,y2;
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            p[i]=node(x1,x2,y1,1);
            p[n+i]=node(x1,x2,y2,-1);
            x[i]=x1,x[i+n]=x2;
        }
        n<<=1;
        sort(p+1,p+1+n);
        sort(x+1,x+1+n);
        memset(tree1,0,sizeof(tree1));
        memset(tree2,0,sizeof(tree2));
        memset(lazy,0,sizeof(lazy));
        double ans=0;
        int sz=unique(x+1,x+1+n)-x-1;
        for(int i=1;i<n;i++)
        {
            int l=lower_bound(x+1,x+1+sz,p[i].x1)-x;
            int r=lower_bound(x+1,x+1+sz,p[i].x2)-x;
            updata(1,1,sz,l,r-1,p[i].sta);
            ans+=tree2[1]*(p[i+1].h-p[i].h);
        }
        printf("%.2f\n",ans);
    }
    return 0;
}