poj1151平面矩形面积并，双离散化

线段树维护扫描线，表示这里的线段树维护多个矩形覆盖同一段区间时不好搞

这里的线段树需要灵活一点

扫描线记录当前扫描所在位置，扫描了多少区域（线段树维护），乘上距离，相加得答案即可

线段树维护，加一条边，在区间内加

删除一条边，区间内删除，tree【i】.s记录区间内被多少线段覆盖

如果没有被线段覆盖，则通过左右子树被线段覆盖的长度相加得到

struct seg//扫描线线段树的一般框架，关键是有一个被覆盖次数，以及特有计算长度的函数
{
         int a,b;
         int len;
         int s; //被覆盖的次数
};

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdlib>
using namespace std;


int n,kase=0;
struct aa
{
	double x,l,r;int t;
	bool operator <(const aa &b)const
	{
		return x<b.x;
	}
}a[2009];
struct aaa
{
	int l,r,s;//s表示覆盖情况，覆盖了几层
	double len;//表示真实长度
}tree[1009];//线段树
double hs[209];//y坐标的hash
void build(int i,int l,int r)
{
	tree[i].l=l;tree[i].r=r;
	tree[i].len=tree[i].s=0;
	if (l+1==r) return ;
	int mid=(l+r)>>1;
	build(i<<1,l,mid);
	build(i<<1|1,mid,r);
}
void pan(int i)//计算长度，如果被覆盖，则为总长，否则为左右儿子长度之和
{
	if (tree[i].s) tree[i].len=hs[tree[i].r]-hs[tree[i].l];//在线段树中，有的区间是被多个矩形覆盖 
	else 												   //这里的s，表示整个区间被覆盖了几次 
	if (tree[i].l+1==tree[i].r) tree[i].len=0;			   //因为保证删除后和插入之前一样，所有没有 
	else tree[i].len=tree[i<<1].len+tree[i<<1|1].len;	   //标记下传，QAQ感觉思路有点乱 
}
void updata(int i,int l,int r,int t)//插入或删除一条线段
{
	if (tree[i].l==l&&tree[i].r==r)
	{
		tree[i].s+=t;
		pan(i);
		return ;
	}
	int mid=(tree[i].l+tree[i].r)>>1;
	if (mid>=r) updata(i<<1,l,r,t);
	else if (mid<=l) updata(i<<1|1,l,r,t);
	else updata(i<<1,l,mid,t),updata(i<<1|1,mid,r,t);
	pan(i);
}
void work()
{
	double aa,bb,cc,dd;
	int tot=0;
	for (int i=1;i<=n;i++)
	{
		scanf("%lf%lf%lf%lf",&aa,&bb,&cc,&dd);
		a[++tot].x=aa;a[tot].l=bb;a[tot].r=dd;a[tot].t=1;hs[tot]=bb;
		a[++tot].x=cc;a[tot].l=bb;a[tot].r=dd;a[tot].t=-1;hs[tot]=dd;
	}
	sort(hs+1,hs+tot+1);
	sort(a+1,a+tot+1);
	int nn=unique(hs+1,hs+tot+1)-hs-1;
	build(1,1,nn);
	int l,r;
	double ans=0;
	for (int i=1;i<=tot;i++)
	{
		l=lower_bound(hs+1,hs+nn+1,a[i].l)-hs;
		r=lower_bound(hs+1,hs+nn+1,a[i].r)-hs;
		ans+=(a[i].x-a[i-1].x)*tree[1].len;
		updata(1,l,r,a[i].t);
	}
	printf("Test case #%d\n",++kase);
	printf("Total explored area: %.2lf\n\n",ans);
}
int main()
{
	while (scanf("%d",&n)!=EOF)
	{
		if (n==0) break;
		sizeof(a,0,sizeof(a));
		sizeof(hs,0,sizeof(hs));
		sizeof(tree,0,sizeof(tree));
		work();
	}
	return 0;
}

暴力：

离散化后，数据范围小，直接枚举，标记所有可行点就行

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;

double x[109],y[109],xx[109],yy[109];

int n,kase=0;
double hsy[209],hsx[209];
bool f[209][209];
void work()
{
	int tot=0;
	for (int i=1;i<=n;i++) 
	{
		scanf("%lf%lf%lf%lf",&x[i],&y[i],&xx[i],&yy[i]);
		hsx[++tot]=x[i];hsy[tot]=y[i];
		hsx[++tot]=xx[i];hsy[tot]=yy[i];
	}
	sort(hsx+1,hsx+tot+1);
	sort(hsy+1,hsy+tot+1);
	int nx=unique(hsx+1,hsx+tot+1)-hsx-1;
	int ny=unique(hsy+1,hsy+tot+1)-hsy-1;
	double ans=0;
	int a,b,c,d;
	for (int i=1;i<=n;i++)
	{
		a=lower_bound(hsx+1,hsx+nx+1,x[i])-hsx;
		b=lower_bound(hsx+1,hsx+nx+1,xx[i])-hsx;
		c=lower_bound(hsy+1,hsy+ny+1,y[i])-hsy;
		d=lower_bound(hsy+1,hsy+ny+1,yy[i])-hsy;
		for (int i=a;i<b;i++)
			for (int j=c;j<d;j++) f[i][j]=true;
	}
	for (int i=1;i<nx;i++)
		for (int j=1;j<ny;j++)
		if (f[i][j]) ans+=(hsy[j+1]-hsy[j])*(hsx[i+1]-hsx[i]);
	printf("Test case #%d\n",++kase);
	printf("Total explored area: %.2lf\n\n",ans);
	return ;
}
int main()
{
	while (scanf("%d",&n)!=EOF)
	{
		if (n==0) break;
		memset(hsx,0,sizeof(hsx));
		memset(hsy,0,sizeof(hsy));
		memset(f,false,sizeof(f));
		work();
	}
	return 0;
}