【HDU 1255】【覆盖的面积】

题目：

给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积. 

 

Input

输入数据的第一行是一个正整数T(1<=T<=100),代表测试数据的数量.每个测试数据的第一行是一个正整数N(1<=N<=1000),代表矩形的数量,然后是N行数据,每一行包含四个浮点数,代表平面上的一个矩形的左上角坐标和右下角坐标,矩形的上下边和X轴平行,左右边和Y轴平行.坐标的范围从0到100000. 

注意:本题的输入数据较多,推荐使用scanf读入数据. 

Output

对于每组测试数据,请计算出被这些矩形覆盖过至少两次的区域的面积.结果保留两位小数. 

Sample Input

2
5
1 1 4 2
1 3 3 7
2 1.5 5 4.5
3.5 1.25 7.5 4
6 3 10 7
3
0 0 1 1
1 0 2 1
2 0 3 1

Sample Output

7.63
0.00

之前做了扫描线的入门题，二者没有很大差别，只是在求和的时候需要注意咱们找的是重复出现次数超过两次及以上的

入门题博客链接：https://mp.csdn.net/postedit/81516859

ac代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define maxn 1005
using namespace  std;
struct edge{
	double l,r;
	double h;
	int f;
}e[maxn<<1];
bool cmp(edge a,edge b)
{
	return a.h<b.h;
}
struct node{
	int l,r;
	int s;
	double len1;
	double len2;
}q[maxn*8];
double x[maxn*2];
void pushup(int cur)
{
	if(q[cur].s) q[cur].len1=x[q[cur].r+1]-x[q[cur].l];
	else if(q[cur].l==q[cur].r) q[cur].len1=0;
	else
	{
		q[cur].len1=q[cur<<1].len1+q[cur<<1|1].len1;
	}
	if(q[cur].s>1)	q[cur].len2=x[q[cur].r+1]-x[q[cur].l];
	else if(q[cur].l==q[cur].r) q[cur].len2=0;
	else if(q[cur].s==1)  q[cur].len2=q[cur<<1].len1+q[cur<<1|1].len1;
	else  q[cur].len2=q[cur<<1].len2+q[cur<<1|1].len2;
}

void build(int l,int r,int cur)
{
	q[cur].l=l;
	q[cur].r=r;
	q[cur].len1=q[cur].len2=0;
	q[cur].s=0;
	if(l==r)
		return;
	int mid=(q[cur].l+q[cur].r)>>1;
	build(l,mid,cur<<1);
	build(mid+1,r,cur<<1|1);
}
void updata(int l,int r,int cur,int xx)
{
	if(q[cur].l==l&&q[cur].r==r)
	{
		q[cur].s+=xx;
		pushup(cur);
		return ; 
	}
	int mid=(q[cur].l+q[cur].r)>>1;
	if(r<=mid) updata(l,r,cur<<1,xx);
	else if(l>mid) updata(l,r,cur<<1|1,xx);
	else
	{
		updata(l,mid,cur<<1,xx);
		updata(mid+1,r,cur<<1|1,xx);
	}
	pushup(cur);
}
int main()
{
	int t;
	scanf("%d",&t);
	while(t--)
	{
		int n;
		scanf("%d",&n);
		int tot=0;
		double x1,x2,y1,y2;
		for(int i=1;i<=n;i++)
		{
			scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
			edge t1,t2;
			t1.l=t2.l=x1;
			t1.r=t2.r=x2;
			t1.h=y1;
			t1.f=1;
			t2.h=y2;
			t2.f=-1;
			e[tot]=t1;
			e[tot+1]=t2;
			x[tot]=x1;
			x[tot+1]=x2;
			tot+=2;
		}
		sort(e,e+tot,cmp);
		sort(x,x+tot);
		int k=1;
		for(int i=1;i<tot;i++)
		{
			if(x[i]!=x[i-1])
				x[k++]=x[i];
		}
		build(0,k-1,1);
		double ans=0.0;
		for(int i=0;i<tot;i++)
		{
			int l=lower_bound(x,x+k,e[i].l)-x;
			int r=lower_bound(x,x+k,e[i].r)-x-1;
			updata(l,r,1,e[i].f);
			ans+=(e[i+1].h-e[i].h)*q[1].len2;	
		}
		printf("%.2f\n",ans);
	}
	return 0;
}