扫描线

线段树要求的被覆盖的线段的总长
由于两个性质:
1.只关心线段树根节点上被矩形覆盖的长度
2.修改操作总是成对出现。
所以不需要lazy标记

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
const int N=1e4+100;
struct Node{
	double y,x1,x2;int id;	
}node[N*2];
bool cmp(Node a,Node b)
{
	return a.y<b.y;
}
struct tree{
	int l,r,cnt;double len;
}t[N*8];
vector<double>nums;
int find(double x)
{
	return lower_bound(nums.begin(),nums.end(),x)-nums.begin();
}
void build(int x,int l,int r)
{
	t[x].l=l;t[x].r=r;t[x].len=0;t[x].cnt=0;
	if(l==r)
	{
		return;
	}
	int mid=(l+r)/2;
	build(x*2,l,mid);
	build(x*2+1,mid+1,r);
}
void pushup(int x)//x可能是根节点 
{
	if(t[x].cnt)
	{
		t[x].len=nums[t[x].r+1]-nums[t[x].l];
	}
	else if(t[x].l!=t[x].r) 
	{
		t[x].len=t[x*2].len+t[x*2+1].len;
	}
	else
	{
		t[x].len=0;
	}
}
void insert(int x,int l,int r,int k)
{
	if(t[x].l>=l&&t[x].r<=r)
	{
		t[x].cnt+=k;
		pushup(x);
		return;
	}
	int mid=(t[x].l+t[x].r)/2;
	if(l<=mid) insert(x*2,l,r,k);
	if(r>mid) insert(x*2+1,l,r,k);
	pushup(x);
}

int main()
{
	int n;int h=0;
	while(scanf("%d",&n)&&n)
	{
		nums.clear();//一定要初始化 
		for(int i=1;i<=n;i++)
		{
			double x1,y1,x2,y2;
			scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
			nums.push_back(x1);
		 	nums.push_back(x2);
			node[i].id=1;node[i].y=y1;node[i].x1=x1;node[i].x2=x2;
			node[i+n].id=-1;node[i+n].y=y2;node[i+n].x1=x1;node[i+n].x2=x2;
		}
		sort(nums.begin(),nums.end());
		nums.erase(unique(nums.begin(),nums.end()),nums.end());
		build(1,0,nums.size()-2);
		sort(node+1,node+2*n+1,cmp);
		double ans=0;
		for(int i=1;i<=2*n;i++)
		{
		
			if(i!=1)
			{
				ans+=(node[i].y-node[i-1].y)*t[1].len;			
			}
			insert(1,find(node[i].x1),find(node[i].x2)-1,node[i].id);
		}
		printf("Test case #%d\n",++h);
		printf("Total explored area: %.2lf\n\n",ans);
	}
	

}

石油大:Black and White
这道题不满足上边第二个性质需要lazy(cnt)标记

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1e4+100;
struct segment{
	int x,y1,y2;
}seg[N*2];
struct tree{
	int l,r,len,cnt;
}t[N*4];
bool cmp(segment a,segment b)
{
	return a.x<b.x;
}
void build(int x,int l,int r)
{
	t[x].l=l;t[x].r=r;t[x].cnt=t[x].len=0;
	if(l==r)
	{
		return;
	}
	int mid=(l+r)/2;
	build(x*2,l,mid);
	build(x*2+1,mid+1,r);
}
void pushup(int x)
{
	t[x].len=t[x*2].len+t[x*2+1].len;
}
void pushdown(int x)
{
	if(t[x].cnt)
	{
		t[x].cnt=0;
		t[x*2].cnt^=1;t[x*2+1].cnt^=1;
		t[x*2].len=t[x*2].r-t[x*2].l+1-t[x*2].len;
		t[x*2+1].len=t[x*2+1].r-t[x*2+1].l+1-t[x*2+1].len;
	}
}
void insert(int x,int l,int r)
{
	if(t[x].l>=l&&t[x].r<=r)
	{
		t[x].len=t[x].r-t[x].l+1-t[x].len;
		t[x].cnt^=1;
		return;
	}
	pushdown(x);
	int mid=(t[x].l+t[x].r)/2;
	if(l<=mid) insert(x*2,l,r);
	if(r>mid) insert(x*2+1,l,r);
	pushup(x);
}
int main()
{
	int n,k;
	int h;
	scanf("%d",&h);
	while(h--)
	{
		scanf("%d%d",&n,&k);
		build(1,1,n);
		for(int i=1;i<=k;i++)
		{
			int x1,x2,y1,y2;
			scanf("%d%d%d%d",&x1,&x2,&y1,&y2);
			seg[i].x=x1;seg[i].y1=y1;seg[i].y2=y2;
			seg[i+k].x=x2+1;seg[i+k].y1=y1;seg[i+k].y2=y2;
		}
		sort(seg+1,seg+2*k+1,cmp);

		int ans=0;
		for(int i=1;i<=2*k;i++)
		{
			if(i!=1)
			{
				ans+=(seg[i].x-seg[i-1].x)*t[1].len;
			}
			insert(1,seg[i].y1,seg[i].y2);
		}
		printf("%d\n",ans);
	}

}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值