矩形有重叠,求总面积。
横纵两条扫描线:
#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));
}
}