题目：
题目链接：[POJ Atlantis]
题解：
（这个题作为学习线段树的尾巴，，，结果还是学了好长时间）
先安利一下sdau_blue大佬的优秀博客，太有用了。
这个题就是有许多的矩形进行随意放置在二维平面内，重叠的面积只算一次，求总面积。
在这里就是离散化处理横坐标，对于纵坐标进行结构体处理一下，然后再以横坐标作为线段(区间)，对横坐标线段进行扫描，扫描的作用是每次更新下底边总长度和下底边个数,增加新面积即可。
代码：
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define D double
using namespace std;
inline int read()
{
int s=0,w=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
while(ch<='9'&&ch>='0')s=s*10+ch-'0',ch=getchar();
return s*w;
}
const int sea=1e5+7;
int T,n,m,t,cnt[sea];
D ans,x[sea],sum[sea];
struct hit
{
D l,r,h;int d;
hit(){} hit(D l,D r,D h,int d):l(l),r(r),h(h),d(d) {}
}tr[sea*4];//结构体函数，，，，lrd的经典
bool cmp(hit x,hit y){return x.h<y.h;}
void up(int l,int r,int i)
{
if(cnt[i]) sum[i]=x[r+1]-x[l];
else sum[i]=sum[i<<1]+sum[i<<1|1];
}
void build(int x,int y,int l,int r,int i,int tt)
{
if(x<=l&&y>=r){cnt[i]+=tt; up(l,r,i);return ;}
int mid=(l+r)/2;
if(x<=mid) build(x,y,l,mid,i<<1,tt);
if(y>mid) build(x,y,mid+1,r,i<<1|1,tt);
up(l,r,i);
}
int main()
{
while(1)
{
T=read(); if(!T) return 0;
n=m=0; ans=0;
memset(sum,0,sizeof(sum));
memset(cnt,0,sizeof(cnt));
for(int i=1;i<=T;i++)
{
D x1,y1,x2,y2;
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
x[++n]=x1; x[++n]=x2;
tr[++m]=hit(x1,x2,y1,1);
tr[++m]=hit(x1,x2,y2,-1);
}
sort(x+1,x+n+1); sort(tr+1,tr+m+1,cmp);
int k=unique(x+1,x+n+1)-x-1;
for(int i=1;i<m;i++)
{
int l=lower_bound(x+1,x+k+1,tr[i].l)-x;
int r=lower_bound(x+1,x+k+1,tr[i].r)-x; r--;//这里的r-1，，还不是很懂，应该是为了防止少算mid———mid-1之间少出现的小数所以上面要加1，下面减1应该是与其抵消的，，还是不很懂
if(l<=r) build(l,r,1,k-1,1,tr[i].d);
ans+=sum[1]*(tr[i+1].h-tr[i].h);
}
printf("Test case #%d\nTotal explored area: %.2f\n\n",++t,ans);
}
return 0;
}