模板：扫描线

矩形有重叠，求总面积。

横纵两条扫描线：

#include<bits/stdc++.h>  
using namespace std;  
  
double x[2002],y[2002];//最多100个矩形，所以最多只有200条横线或纵线  
double a[1002][4];//矩形实际坐标，分别是左下角横、纵，右下角横、纵
bool cover[2002][2002];  
const double EPS = 1e-7;  
  
int cmp(const void * b,const  void * a){//b>a 返回1  
    if(fabs(*(double* )b-*(double *)a)<EPS)return 0;  
    else if(*(double* )b-*(double *)a>0)return 1;  
    else return -1;  
}  
  
int main(){  
    int n,i,k,j,h1,h2,v1,v2;  
    while(scanf("%d",&n)&&n){  
        for(i=0,k=0;i<n;i++,k++){  
            scanf("%lf%lf%lf%lf",&a[i][0],&a[i][1],&a[i][2],&a[i][3]);  
            x[k]=a[i][0],y[k]=a[i][1],x[++k]=a[i][2],y[k]=a[i][3];  
        }  
        qsort(x,2*n,sizeof(x[0]),cmp);  
        qsort(y,2*n,sizeof(y[0]),cmp);//将横线，竖线排好序  
        memset(cover,0,sizeof(cover));//初始时无覆盖区  
        //对每个矩形，进行一次横纵线的扫描，覆盖范围  
        //扫描线条数作为下标，离散化  
        for(i=0;i<n;i++){  
            k=0;  
            while(fabs(x[k]-a[i][0])>EPS)k++;h1=k;  
            k=0;  
            while(fabs(y[k]-a[i][1])>EPS)k++;v1=k;  
            k=0;  
            while(fabs(x[k]-a[i][2])>EPS)k++;h2=k;  
            k=0;  
            while(fabs(y[k]-a[i][3])>EPS)k++;v2=k;  
            //标记该矩形覆盖的区域  
            for(j=h1;j<h2;j++)      //不考虑最右  
                for(k=v1;k<v2;k++)  //不考虑最下  
                    cover[j][k]=true;  
        }  
        double sum=0;  
        //这时再用扫描线真实的位置计算面积，由矩形左下角的点(判断是否覆盖)及四条扫描线确定。  
        for(i=0;i<2*n-1;i++)  
            for(j=0;j<2*n-1;j++)  
                sum+=(cover[i][j]*(x[i+1]-x[i])*(y[j+1]-y[j]));  
        printf("%d\n", (int)round(sum));  
    }  
}