题意：

求所给矩形的覆盖面积

题解：

利用扫描线的思想，先将坐标离散化，之后以y轴分成多个矩形求解，可以让下边界+1上边界-1

问题就转化为了：求区间中有多少个非0数，要求支持区间+1 -1操作

我们可以通过维护区间最小值以及最小值的个数来完成这件事情

代码：

#include
      
       #include
       
        #include
        
         #include
         
          #include 
          
           #define maxn 5000#define INF 59999999#define eps 1e-6#define mid (h+t)/2using namespace std;struct re {    double a;int b;}a[maxn],b[maxn];struct ree{    int h,t,x,lazy;    double sum,tot;}p[maxn*4];double c[maxn];int a1[maxn];bool cmp(re x,re y){  if (x.a
           
            t|| p[x].t
            
             t||p[x].t
             
              >n&&n!=0){ o++; memset(p,0,sizeof(p)); for (int i=1;i<=2*n;i++) a[i].b=i,b[i].b=i; for (int i=1;i<=n;i++) { cin>>a[2*i-1].a>>b[2*i-1].a>>a[2*i].a>>b[2*i].a; } sort(a+1,a+1+2*n,cmp); sort(b+1,b+1+2*n,cmp); int ll=0; a[0].a=INF; for (int i=1;i<=2*n;i++) { if (abs(a[i].a-a[i-1].a)>eps) ll++; a1[a[i].b]=ll; c[ll]=a[i].a; } double ans=0; build(1,1,ll-1); for (int i=1;i<2*n;i++) { int pp,tmp=b[i].b; if (tmp%2==1) pp=1 ;else pp=-1; insert(1,a1[(tmp+1)/2*2-1],a1[(tmp+1)/2*2]-1,pp); ans+=(b[i+1].a-b[i].a)*query(1,1,ll-1); } cout<<"Test case #"<
              
               <
               
                <<"Total explored area: " ; printf("%.2f\n\n",ans); }}