hdu 3642 Get The Treasury(扫描线、立方体交)

hdu 3642 Get The Treasury

和此题做对比

这个题相比之有不同的就是需要求覆盖两次以上部分的体积，而且不一定底面都是在xy平面上，但处理思路依然是一样的

主要问题是如果沿用上次求面积交的方法没办法处理处覆盖多次的情况，于是又回到了老办法上，没有多想先写下试试，结果没有超时，就这样了

不过还是要研究下更快的求面积交的算法，感觉挺重要

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<map>
using namespace std;
#define MAXN 2005
#define ll long long
#define lch p<<1
#define rch p<<1|1
#define mid (t[p].l+t[p].r)>>1
struct node
{
	int l,r;
	int cover;
	ll dy;
}t[MAXN<<2];
struct segment
{
	ll x,yd,yu;
	ll zd,zu;
	int lr;
	segment(){}
	segment(ll x,ll yd,ll yu,ll zd,ll zu,int lr):x(x),yd(yd),yu(yu),zd(zd),zu(zu),lr(lr){}
	bool operator < (const segment &a) const
	{
		return x<a.x;
	}
}seg[MAXN],tseg[MAXN];
ll yy[MAXN],zz[MAXN];
void construct(int l,int r,int p)
{
	t[p].l=l,t[p].r=r;
	t[p].cover=t[p].dy=0;
	if(l+1>=r) return ;
	int m=mid;
	construct(l,m,lch);
	construct(m,r,rch);
}
void modify(int l,int r,int op,int p)
{
	if(t[p].l==l&&t[p].r==r)
	{
		t[p].cover+=op;
		return ;
	}
	int m=mid;
	if(r<=m) modify(l,r,op,lch);
	else if(l>=m) modify(l,r,op,rch);
	else modify(l,m,op,lch),modify(m,r,op,rch);
}
void query(ll &ans,int p)
{
	if(t[p].cover>2)
	{
		ans+=yy[t[p].r]-yy[t[p].l];
	}
	else if(t[p].l+1<t[p].r)
	{
		t[lch].cover+=t[p].cover;
		t[rch].cover+=t[p].cover;
		query(ans,lch);
		query(ans,rch);
		t[lch].cover-=t[p].cover;
		t[rch].cover-=t[p].cover;
	}
}
map<ll,int> rec;
int main()
{
	int cas,n;
	ll x1,y1,z1,x2,y2,z2;
	scanf("%d",&cas);
	for(int kcas=1;kcas<=cas;kcas++)
	{
		scanf("%d",&n);
		int i;
		for(i=0;i<(n<<1);i+=2)
		{
			scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&x1,&y1,&z1,&x2,&y2,&z2);
			seg[i]=segment(x1,y1,y2,z1,z2,1);
			seg[i+1]=segment(x2,y1,y2,z1,z2,-1);
			yy[i]=y1,yy[i+1]=y2;
			zz[i]=z1,zz[i+1]=z2;
		}
		sort(seg,seg+i);
		sort(yy,yy+i);
		sort(zz,zz+i);
		int up=unique(yy,yy+i)-yy;
		int upz=unique(zz,zz+i)-zz;
		rec.clear();
		for(int i=0;i<up;i++)
			rec[yy[i]]=i;
		ll ans=0;
		for(int j=0;j<upz-1;j++)
		{
			ll area=0;
			int tn=0;
			for(int i=0;i<(n<<1);i++)
				if(seg[i].zd<=zz[j]&&seg[i].zu>=zz[j+1])
					tseg[tn++]=seg[i];
			if(tn==0) continue; 
			construct(0,up-1,1);
			for(int i=0;i<tn-1;i++)
			{
				modify(rec[tseg[i].yd],rec[tseg[i].yu],tseg[i].lr,1);
				ll tans=0;
				query(tans,1);
				area+=(tseg[i+1].x-tseg[i].x)*tans;
			}
			ans+=area*(zz[j+1]-zz[j]);
		}
		printf("Case %d: %I64d\n",kcas,ans);
	}
	return 0;
}