矩形覆盖求并 二维离散化 扫描线算法

Problem:
给了n个矩形的左下角坐标和右上角坐标，求这n个矩形的面积和，覆盖部分只计算一次。
2017广西南宁acm区域赛:https://nanti.jisuanke.com/t/17313
Solution:
利用扫描线算法进行离散化，然后再求和。

#include<cstdio>
#include<iostream>
#include<sstream>
#include<cstdlib>
#include<cmath>
#include<cctype>
#include<string>
#include<cstring>
#include<algorithm>
#include<stack>
#include<queue>
#include<set>
#include<map>
#include<ctime>
#include<vector>
#include<fstream>
#include<list>
#include <iomanip>

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(){
    //    freopen("/Users/really/Documents/code/input","r",stdin);

    int n,i,k,j,h1,h2,v1,v2;
    while(scanf("%d",&n) != EOF){
        if(n==0) {
            printf("*\n");
            break;
        }
        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));
    }
}