线段树+扫描线【HDU1542】Atlantis

Description

给定一些二维空间上的矩形，求它们的面积并.

一道线段树+扫描线的板子题

然而即使我会打了,也不能灵活运用这种算法。QAQ

遇到题还是不太会.

但是这种板子题还是随随便便切的.

时间特别特别多的话,就练线段树！(flag)

这里放下代码好了,网上讲解有很多.就不写了 qwq

代码

#include<cstdio>
#include<algorithm>
#define ls o<<1
#define rs o<<1|1
#define R register
using namespace std;
struct code
{
	double l,r,h;
	int f;
	bool operator < (const code&a)const
	{
		return h<a.h;
	}
}edge[2008];
struct cod
{
	int l,r,s;
	double len;
}tr[8888];
double x[2333];
void build(int o,int l,int r)
{
	tr[o].l=l,tr[o].r=r;
	tr[o].s=0;tr[o].len=0;
	if(l==r)return;
	int mid=(l+r)>>1;
	build(ls,l,mid);
	build(rs,mid+1,r);
}
inline void up(int o)
{
	if(tr[o].s)
		tr[o].len=x[tr[o].r+1]-x[tr[o].l];
	else if(tr[o].l==tr[o].r)
		tr[o].len=0;
	else tr[o].len=tr[ls].len+tr[rs].len;
}
void change(int o,int l,int r,int del)
{
	if(tr[o].l==l && tr[o].r==r)
	{
		tr[o].s+=del;
		up(o);
		return;
	}
	int mid=(tr[o].l+tr[o].r)>>1;
	if(r<=mid)change(ls,l,r,del);
	else if(l>mid) change(rs,l,r,del);
	else change(ls,l,mid,del),change(rs,mid+1,r,del);
	up(o);
}
int main()
{
	int n,cas=0;
	for(;;)
	{
		scanf("%d",&n);
		if(n==0)break;
		int tot=0;
		for(R int i=1;i<=n;i++)
		{
			R double x1,x2,y1,y2;
			scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
			edge[++tot].l=x1;edge[tot].r=x2;edge[tot].h=y1;
			edge[tot].f=1;x[tot]=x1;
			edge[++tot].l=x1;edge[tot].r=x2;edge[tot].h=y2;
			edge[tot].f=-1;x[tot]=x2;
		}
		sort(x+1,x+tot+1);
		sort(edge+1,edge+tot+1);
		int new_n=1;
		for(R int i=2;i<=tot;i++)
			if(x[new_n]!=x[i])x[++new_n]=x[i];
		build(1,1,new_n);
		double ans=0;
		for(R int i=1;i<=tot;i++)
		{
			int l=lower_bound(x+1,x+new_n+1,edge[i].l)-x;
			int r=lower_bound(x+1,x+new_n+1,edge[i].r)-x-1;
			change(1,l,r,edge[i].f);
			ans+=(edge[i+1].h-edge[i].h)*tr[1].len;
		}
		printf("Test case #%d
",++cas);
		printf("Total explored area: %.2f

",ans);
	}
}