[HDU1542]Atlantis

Description:

给你n个矩形,求矩形面积并

Hint:

\(n\le 100\)

Solution:

扫描线模板

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#define ls p<<1 
#define rs p<<1|1
using namespace std;
typedef long long ll;
const int mxn=1e5+5;
int n,m,tot,tot,cnt,hd[mxn];
double tr[mxn<<2];
int cov[mxn<<2];
double x[mxn];

inline int read() {
	char c=getchar(); int x=0,f=1;
	while(c>'9'||c<'0') {if(c=='-') f=-1;c=getchar();}
	while(c<='9'&&c>='0') {x=(x<<3)+(x<<1)+(c&15);c=getchar();}
	return x*f;
}
inline void chkmax(int &x,int y) {if(x<y) x=y;}
inline void chkmin(int &x,int y) {if(x>y) x=y;}

struct T {
	double l,r,h; int val;
	friend bool operator < (T x,T y) {
		return x.h<y.h;
	}
}t[mxn];

void push_up(int p,int l,int r) {
	if(cov[p]) tr[p]=x[r+1]-x[l];
	else if(l==r) tr[p]=0;
	else tr[p]=tr[ls]+tr[rs];
}

void update(int l,int r,int ql,int qr,int val,int p) {
	if(ql<=l&&r<=qr) {
		cov[p]+=val;
		push_up(p,l,r);
		return ;
	} 
	int mid=(l+r)>>1;
	if(ql<=mid) update(l,mid,ql,qr,val,ls);
	if(qr>mid) update(mid+1,r,ql,qr,val,rs);
	push_up(p,l,r);
}

int main()
{
	int cas=1;
	while(scanf("%d",&n)!=EOF) {
		if(n==0) break;
		cnt=tot=0; 
		for(int i=0;i<n;++i) {
			double a,b,c,d;
			scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
			t[++cnt]=(T){a,c,b,1}; x[cnt]=a;
			t[++cnt]=(T){a,c,d,-1}; x[cnt]=c;
		}
		sort(x+1,x+cnt+1);
		sort(t+1,t+cnt+1);
		++tot;
		for(int i=2;i<=cnt;++i) 
			if(x[i]!=x[i-1]) x[++tot]=x[i];
		memset(tr,0,sizeof(tr));
		memset(cov,0,sizeof(cov));
		double ans=0;
		for(int i=1;i<cnt;++i) {
			int l=lower_bound(x+1,x+tot+1,t[i].l)-x;
			int r=lower_bound(x+1,x+tot+1,t[i].r)-x-1;
			update(1,tot,l,r,t[i].val,1);
			ans+=tr[1]*(t[i+1].h-t[i].h);
		}
		printf("Test case #%d\nTotal explored area: %.2lf\n\n",cas++,ans);
	}
    return 0;
}