POJ 1151 Atlantis(线段树+离散化+扫描线)

题目很短，大意为有多个矩形求合计面积。（会有重叠，要剔除）

1.线段树，将y轴当作线段，然后在计算的时候从x轴扫描过去

2.数据离散化，而这里区别于前面的题目，只有有区间才可能有长度，所以这里的叶子节点不是l=r,而是l+1=r

3.扫描线，这是我最郁闷的一个地方。刚开始不理解。其实质上是从x轴扫描过去，（把一段段平行于y轴的线段表示成数组 ）遇到了看线段树中的标记位。若大于0则为正，表明覆盖了这段线，否则表明没有覆盖或者重复。

具体过程：

2
10 10 20 20
15 25 25 35
0
//分离
2
10 10 20 20
15 15 25 25
0
//相交
2
10 10 20 20
12 12 18 18
0
//包含

这几个模拟一下。都是算(xi-1,x)段内，y的长度

注意：

对于线段树明确是结点还是区间，这个有根本区别，尽量从1开始取下标。节点数一般为3*N，但保险还是4*N

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <algorithm>
#define maxn 210
#define mem(a) memset(a,0,sizeof(a))

struct node
{
	int l,r;
	double lf,rf;
	int cover;
	double len;
}nodes[maxn*4];
/*
离散化了区间l,r，而其实际长度lf,rf还需保存。因为计算长度用到
cover为标记位，若大于0则表明覆盖，Calen时计算
len为真实长度
*/
struct Line
{
	double x,y1,y2;
	int tag;
}lines[maxn];
/*
用于每次Update，将一个矩形的前后两条边分为两条线段，前面的tag=1,后面-1.tag>0就计算长度
这个原因可以画图得知。大于0的情况出现于多个矩形占据同一x点。
*/
 double y[maxn];
int cmp(const void *a, const void *b)
{
	if( (*(struct Line *)a).x > (*(struct Line *)b).x)
		return 1;
	else
		return -1;
}
void Build(int l,int r,int t)
{
	nodes[t].l=l;
	nodes[t].r=r;
	nodes[t].lf=y[l];
	nodes[t].rf=y[r];
	if(l+1==r)
		return;//这里格外注意，求面积只有区间才有意义，当到l==r的叶子无意义，同时与后面Calen相对应
	int mid=(l+r)>>1;
	Build(l,mid,t<<1);
	Build(mid,r,t<<1|1);
}
void Calen(int t)
{
	if(nodes[t].cover>0)
		nodes[t].len=nodes[t].rf-nodes[t].lf;
	else if(nodes[t].l+1==nodes[t].r) nodes[t].len=0;
	else
		nodes[t].len=nodes[t<<1].len+nodes[t<<1|1].len;
}
/*
对于线段树的离散化，y[l]=y[r],y[l]<y[r]两种情况。但不论哪种第一步判断cover都解决了问题。
如果不能再分(nodes[t].l+1==nodes[t].r)则直接为0，否则等于两部分之和
*/
 void Update(int t,struct Line L)
{
	if(L.y1==nodes[t].lf && L.y2==nodes[t].rf)
	{
		nodes[t].cover+=L.tag;
		Calen(t);return;
	}
	if(L.y2<=nodes[t<<1].rf) Update(t<<1,L);
	else if(L.y1>=nodes[t<<1|1].lf) Update(t<<1|1,L);
	else
	{
		struct Line temp;
		temp=L;
		temp.y2=nodes[t<<1].rf;
		Update(t<<1,temp);
		temp=L;
		temp.y1=nodes[t<<1|1].lf;
		Update(t<<1|1,temp);
	}
	Calen(t);
}

int main()
{
	int n,m,Case=1;
	int i;
	double a,b,c,d;
	while (scanf("%d",&n),n)
	{
		m=1;
		for(i=0;i<n;i++)
		{
			scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
			lines[m].tag=1;lines[m].x=a;
			lines[m].y1=b; lines[m].y2=d;
			y[m]=b;m++;
			lines[m].tag=-1;lines[m].x=c;
			lines[m].y1=b; lines[m].y2=d;
			y[m]=d;m++;
		}
		std::sort(y+1,y+m);//这里注意浮点数的判断问题，尽量使用不等于的符号，直接用减号会在返回时转int出错，也避免==,WA了好久
		qsort(&lines[1],m-1,sizeof(lines[0]),cmp);
		Build(1,m-1,1);
		Update(1,lines[1]);
		a=0;
		for(i=2;i<m;i++)
		{
			a+=nodes[1].len*(lines[i].x-lines[i-1].x);
			Update(1,lines[i]);
		}
		printf("Test case #%d\nTotal explored area: %.2f\n\n",Case++,a);
	}
	return 0;
}