HDOJ 1542 Atlantis

线段树扫描线的应用，每一个矩形将其抽象为上下两条边，上下边区分为f=-1，f=1，离散化横坐标，按横坐标建立线段树，然后将每条边按纵坐标高度排序。

接下来的操作可以想象为一条线由下往上扫描，每次扫过两纵坐标之间的区域。用cover[rt]记录出现下边比上边多的次数，当cover[rt]不为0时表示此段纵坐标区间内rt表示横坐标之间的区域矩形面积可以计算。

每次pushup，sum表示当前区间可计算的区间长度。因每次区间可计算值将更新，所以每次操作后累加当前区域面积值。最后输出累加值即可。

注意节点并不表示端点，表示当前端点到下一段点间的区间。注意r+1,r-1的问题。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 220
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
double sum[N<<2];
int cover[N<<2];
double X
;
struct BIAN
{
double l;
double r;
double y;
int f;
}s
;
bool cmp(BIAN a,BIAN b)
{
return a.y<b.y;
}
void pushup(int l,int r,int rt)
{
if(cover[rt])sum[rt]=X[r+1]-X[l];
else if(l==r)sum[rt]=0;
else sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void update(int L,int R,int c,int l,int r,int rt)
{
if(L<=l&&r<=R)
{
cover[rt]+=c;
pushup(l,r,rt);
return ;
}
int m=(l+r)>>1;
if(L<=m)update(L,R,c,lson);
if(R>m)update(L,R,c,rson);
pushup(l,r,rt);
}
int main()
{
int n,i,cas=1;
while(scanf("%d",&n),n)
{
i=0;
double a,b,c,d;
while(n--)
{
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
X[i]=a;
s[i].y=b;s[i].l=a;s[i].r=c;s[i].f=1;
i++;
X[i]=c;
s[i].y=d;s[i].l=a;s[i].r=c;s[i].f=-1;
i++;
}
int m=i;
sort(X,X+i);
sort(s,s+i,cmp);
int lx=unique(X,X+i)-X;;
memset(sum,0,sizeof(sum));
memset(cover,0,sizeof(cover));
double ret=0;
for(i=0;i<=m-1;i++)
{
int l=lower_bound(X,X+lx,s[i].l)-X;
int r=lower_bound(X,X+lx,s[i].r)-X-1;
update(l,r,s[i].f,0,lx-1,1);
ret+=sum[1]*(s[i+1].y-s[i].y);
}
printf("Test case #%d\nTotal explored area: %.2lf\n\n",cas++,ret);

}
return 0;
}