hdu 1542 Atlantis(线段树+扫描线)

先离散y，然后以y为区间建线段树，从左到右，枚举分割线x，一段一段面积累加。

Cover记录当前的y区间被几个矩形包含。len记录被覆盖的区间长度。

比较机智的地方是Cover并不需要在更新线段树的时候进行更新。

一旦这个区间的cover更改了，他可以由之前的小区间的值反更新上来。

#include<iostream>
#include<algorithm>

using namespace std;
#define lson 2*o,l,mid
#define rson 2*o+1,mid+1,r
#define INF 1000000000
const int maxn=5010;
typedef long long ll;
struct TREE
{
	double len;
	int cover;
}tree[maxn];
double y[maxn];
struct LINE
{
	double x,y1,y2;
	int flag;
}line[maxn];
int find(int l,int r,double yy)
{
	int ans=-1;
	while(l<=r){
		int mid=(l+r)>>1;
		if(y[mid]==yy)return mid;
		if(y[mid]<yy)
			l=mid+1;
		else
			r=mid-1;
	}
	return ans;
}//找到离散之后的位置
bool cmp(LINE a,LINE b){
	return a.x<b.x;
}
void update(int o,int l,int r)
{
	if(tree[o].cover){//之前有过标记
		tree[o].len=y[r+1]-y[l];//为区间长度
	}
	else{//没有标记
		if(l==r)tree[o].len=0;
		else tree[o].len=tree[2*o].len+tree[2*o+1].len;
	}
}
void insert(int o,int l,int r,int a,int b,int cover){
	if(b<l||a>r)return;
	if(a<=l&&r<=b){
		tree[o].cover+=cover;
		update(o,l,r);//
		return;
	}
	int mid=(l+r)>>1;
	insert(lson,a,b,cover);
	insert(rson,a,b,cover);
	update(o,l,r);
}
void init()
{
	memset(line,0,sizeof(line));
	memset(y,0,sizeof(y));
	memset(tree,0,sizeof(tree));
}
int main()
{
	int n;int ica=1;
	while(scanf("%d",&n),n){
		init();
		double x1,x2,y1,y2;
		for(int i=1;i<=2*n;i+=2){
			scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
			line[i].x=x1;line[i+1].x=x2;
			line[i+1].y1=line[i].y1=y1;
			line[i+1].y2=line[i].y2=y2;
			line[i].flag=1;line[i+1].flag=-1;
			y[i]=y1;y[i+1]=y2;
		}
		sort(y+1,y+2*n+1);
		int num=unique(y+1,y+2*n+1)-y-1;//纵向区间个数
		sort(line+1,line+2*n+1,cmp);
		double area=0;
		for(int i=1;i<2*n;i++){
			int l=find(1,num,line[i].y1);
			int r=find(1,num,line[i].y2);
			insert(1,1,num,l,r-1,line[i].flag);
			area+=tree[1].len*(line[i+1].x-line[i].x);
		}
		printf("Test case #%d\n",ica++);
		printf("Total explored area: %.2lf\n\n",area);
	}
	return 0;
}