AcWing 247. 亚特兰蒂斯

有几个古希腊书籍中包含了对传说中的亚特兰蒂斯岛的描述。

其中一些甚至包括岛屿部分地图。

但不幸的是，这些地图描述了亚特兰蒂斯的不同区域。

您的朋友Bill必须知道地图的总面积。

你自告奋勇写了一个计算这个总面积的程序。

输入格式
输入包含多组测试用例。

对于每组测试用例，第一行包含整数n，表示总的地图数量。

接下来n行，描绘了每张地图，每行包含四个数字x1,y1,x2,y2（不一定是整数）,(x1,y1)和(x2,y2)分别是地图的左下角位置和右上角位置。

当输入用例n=0时，表示输入终止，该用例无需处理。

输出格式
每组测试用例输出两行。

第一行输出”Test case #k”，其中k是测试用例的编号，从1开始。

第二行输出“Total explored area：a”，其中a是总地图面积（即此测试用例中所有矩形的面积并，注意如果一片区域被多个地图包含，则在计算总面积时只计算一次），精确到小数点后两位数。

在每个测试用例后输出一个空行。

数据范围
1≤n≤100,//数据应该给错了，n会再大一点。怪不得我一直Segmentation Fault。。。
0≤x1<x2≤100000,
0≤y1<y2≤100000
输入样例：
2
10 10 20 20
15 15 25 25.5
0
输出样例：
Test case #1
Total explored area: 180.00

注意
①double不能++
②同一数据输出时应该在同一语句输出，不然printf和cout混用输出顺序会出错
③ios::sync_with_stdio(false);很好用

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=400;
int n,tot,cc;
double d[N],ans;
struct node
{
	double x,yu,yd;
	int k;
}b[N];
bool cmp(node a,node b)
{
	return a.x<b.x;
}

struct segmentree
{
	int l,r,cnt;
	double len;
	#define l(p) p*2
	#define r(p) p*2+1
}t[N*4];

void decrete()
{
	sort(d+1,d+1+2*n);
	tot=unique(d+1,d+1+2*n)-d-1;
}
int val(double y)
{
	return lower_bound(d+1,d+1+tot,y)-d;
}
double raw(int y)
{
	return d[y];
}
void build(int p,int l,int r)
{
	t[p].l=l;t[p].r=r;
	if(l==r)return ;
	int mid=(l+r)>>1;
	build(l(p),l,mid);build(r(p),mid+1,r);
}
void get_len(int p)
{
	if(t[p].cnt)t[p].len=raw(t[p].r+1)-raw(t[p].l);
	else t[p].len=t[l(p)].len+t[r(p)].len;
}
void change(int p,int l,int r,int v)
{
	if(t[p].l==l&&t[p].r==r)
	{
		t[p].cnt+=v;
        get_len(p);
		return ;
	}
	int mid=(t[p].l+t[p].r)>>1;
	if(r<=mid)change(l(p),l,r,v);
	else if(l>mid)change(r(p),l,r,v);
	else 
	{
		change(l(p),l,mid,v);
		change(r(p),mid+1,r,v);
	}
	get_len(p);
}

int main()
{
	ios::sync_with_stdio(false);
	while(cin>>n&&n){
		
		double x,y,xx,yy;ans=0;tot=0;
		memset(d,0,sizeof(d));
		memset(b,0,sizeof(b));
		memset(t,0,sizeof(t));
		for(int i=1; i<=n; i++)
		{
			cin>>x>>yy>>xx>>y;
			d[i]=y;d[i+n]=yy;
			b[i].x=x,b[i].yu=y,b[i].yd=yy,b[i].k=1;
			b[i+n].x=xx,b[i+n].yu=y,b[i+n].yd=yy,b[i+n].k=-1;
		}

		decrete();
		n*=2;
		build(1,1,tot-1);
		sort(b+1,b+1+n,cmp);
		for(int i=1; i<n; i++)
		{
			change(1,val(b[i].yd),val(b[i].yu)-1,b[i].k);
			ans+=t[1].len*(b[i+1].x-b[i].x);
		}
		printf("Test case #%d\nTotal explored area: %.2lf\n\n",++cc,ans);
	}
	
	return 0;
}