hdu 1255 覆盖的面积 矩阵面积交

题意：

求n个矩阵面积交

分析：

其实和求矩阵面积并是一样的思想

我们只需增加一个nsum数组储存区间中覆盖2次及以上的长度就行

#include <bits/stdc++.h>#define maxn 200002#define tmp (st<<1)#define mid ((l+r)>>1)#define lson l,mid,tmp#define rson mid+1,r,tmp|1using namespace std;int cnt[maxn<<2];double sum[maxn<<2];double nsum[maxn<<2];double x[maxn];struct Seg{    double h,l,r;    int s;    Seg(){}    Seg(double a,double b,double c,int d):l(a),r(b),h(c),s(d){}    bool operator<(const Seg &cmp)const{        return h<cmp.h;    }}ss[maxn];void push_up(int st,int l,int r){    if(cnt[st]==1){        sum[st]=x[r+1]-x[l];        nsum[st]=sum[tmp]+sum[tmp|1];    }    else if(cnt[st]>=2)sum[st]=nsum[st]=x[r+1]-x[l];    else if(l==r)nsum[st]=sum[st]=0;    else {        sum[st]=sum[tmp]+sum[tmp|1];        nsum[st]=nsum[tmp]+nsum[tmp|1];    }}void update(int L,int R,int c,int l,int r,int st){    if(L<=l&&r<=R){        cnt[st]+=c;        push_up(st,l,r);        return ;    }    if(L<=mid)update(L,R,c,lson);    if(R>mid)update(L,R,c,rson);    push_up(st,l,r);}int main(){    int n,m,loop;    scanf("%d",&loop);    while(loop--){        scanf("%d",&n);        double a,b,c,d,ans=0;        m=0;        while(n--){            scanf("%lf%lf%lf%lf",&a,&b,&c,&d);            x[m]=a;            ss[m++]=Seg(a,c,b,1);            x[m]=c;            ss[m++]=Seg(a,c,d,-1);        }        sort(x,x+m);        sort(ss,ss+m);        for(int i=0;i<m;++i){            int l=lower_bound(x,x+m,ss[i].l)-x;            int r=lower_bound(x,x+m,ss[i].r)-x-1;            update(l,r,ss[i].s,0,m-1,1);            ans+=nsum[1]*(ss[i+1].h-ss[i].h);        }        printf("%.2lf\n",ans);    }    return 0;}