hdu 1542 Atlantis

题目大意：

求矩形的并面积；

解题思路：

从下往上扫描，实际上就是把整个图案分成了一横块，一横块的，矩形的左边的边，标记为1，右边的边标记为-1，用于对线段长度的处理，每次查询的实际上上是横块的横边长度。因为范围略大，对横坐标要进行离散，需要注意的时候线段树中【l,r】并不是表示表示pos【l】，pos【r】的长度，是pos【l】，pos【r+1】的长度，因为离散后的建立的线段树的一个点是一个单位长度。网上一个大牛画的的图1----------2----------3----------4离散的点

                                            1               2             3 建立线段树的点

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

#define maxn 100010
#define mn 1100
#define L(a) (a)<<1
#define R(a) (a)<<1|1

struct Line{
	double l,r,h;
	int f;
	Line(){}
	Line(double a,double b,double c,double d):l(a),r(b),h(c),f(d){}
	bool operator <(const Line &a)const{return h<a.h;}
}line[2*mn];

struct Node{
	int lc,rc;
	int cnt;
	double len;
	int mid()
	{
		return (lc+rc)>>1;
	}
}node[2010];

double pos[2*mn];
int num;

int find(double val,int l,int r)
{
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(pos[mid]==val) return mid;
		else if(val<pos[mid]) r=mid-1;
		else l=mid+1;
	}
	return -1;
}
void pushup(int k)
{
	if(node[k].cnt)
	node[k].len=pos[node[k].rc+1]-pos[node[k].lc];
	else if(node[k].lc==node[k].rc)
	node[k].len=0;
	else
	node[k].len=node[L(k)].len+node[R(k)].len;
}
void create(int k,int l,int r)
{
	node[k].lc=l;
	node[k].rc=r;
	node[k].len=node[k].cnt=0;
	if(l==r)
	return ;
	int mid=node[k].mid();
	create(L(k),l,mid);
	create(R(k),mid+1,r);
}

void updata(int k,int l,int r,int val)
{
	if(node[k].lc==l&&node[k].rc==r)
	{
		node[k].cnt+=val;
		pushup(k);
		return ;
	}
	int mid=node[k].mid();
	if(r<=mid)
	updata(L(k),l,r,val);
	else if(l>mid)
	updata(R(k),l,r,val);
	else
	{
		updata(L(k),l,mid,val);
		updata(R(k),mid+1,r,val); 
	}
	pushup(k);
}

int main()
{
	int count=0;
	int n;
	while(~scanf("%d",&n),n)
	{
		num=0;
		for(int i=0;i<n;i++)
		{
			double x1,x2,y1,y2;
			scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
			line[num]=Line(x1,x2,y1,1);
			pos[num++]=x1;
			line[num]=Line(x1,x2,y2,-1);
			pos[num++]=x2;
		}
		sort(line,line+num);
		sort(pos,pos+num);
		int m=1;
		for(int i=1;i<num;i++)
		if(pos[i]!=pos[i-1])
		pos[m++]=pos[i];
		create(1,0,m-1);
		double ans=0;
		for(int i=0;i<num-1;i++)
		{
			int l=find(line[i].l,0,m-1);
			int r=find(line[i].r,0,m-1)-1;  
			updata(1,l,r,line[i].f);
			ans+=(line[i+1].h-line[i].h)*node[1].len;
		}
		printf("Test case #%d\n",++count);
		printf("Total explored area: %.2lf\n",ans);
	}
}