线段树要求的被覆盖的线段的总长
由于两个性质:
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);
}
}