题目要求的时统计覆盖过两次以上部分的面积,query判断修改下就行了,除此之外和Atlantis完全一样
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
#define MAXN 1005
#define ll p<<1
#define rr p<<1|1
struct node
{
int l,r;
int cover;
}t[MAXN*8];
struct edge
{
double x,yd,yu;
int lr;
bool operator < (const edge a) const {return x<a.x;}
}e[MAXN*2];
double yy[MAXN*2];
void construct(int l,int r,int p)
{
t[p].l=l,t[p].r=r;
t[p].cover=0;
if(r-l<=1) return ;
int m=(l+r)>>1;
construct(l,m,ll);
construct(m,r,rr);
}
void modify(int l,int r,int val,int p)
{
if(t[p].l==l&&t[p].r==r)
{
t[p].cover+=val;
return ;
}
int m=(t[p].l+t[p].r)>>1;
if(r<=m) modify(l,r,val,ll);
else if(l>=m) modify(l,r,val,rr);
else modify(l,m,val,ll),modify(m,r,val,rr);
}
void query(double &ans,int p)
{
if(t[p].cover>1)
{
ans+=yy[t[p].r]-yy[t[p].l];
return ;
}
else if(t[p].r-t[p].l>1)
{
t[ll].cover+=t[p].cover,t[rr].cover+=t[p].cover;
query(ans,ll);
query(ans,rr);
t[ll].cover-=t[p].cover,t[rr].cover-=t[p].cover;
}
}
int main()
{
map<double,int> hash;
int cas,n,i;
double x1,x2,yd,yu;
scanf("%d",&cas);
while(cas--)
{
hash.clear();
scanf("%d",&n);
for(i=0;i<n*2;i+=2)
{
scanf("%lf%lf%lf%lf",&x1,&yd,&x2,&yu);
e[i].x=x1,e[i].yd=yd,e[i].yu=yu,e[i].lr=1;
e[i+1].x=x2,e[i+1].yd=yd,e[i+1].yu=yu,e[i+1].lr=-1;
yy[i]=yd,yy[i+1]=yu;
}
sort(e,e+i);
sort(yy,yy+i);
int up=int (unique(yy,yy+i)-yy);
for(i=0;i<up;i++)
hash[yy[i]]=i;
construct(0,up-1,1);
double ans=0,dy;
for(i=0;i<n*2-1;i++)
{
modify(hash[e[i].yd],hash[e[i].yu],e[i].lr,1);
dy=0;
query(dy,1);
ans+=dy*(e[i+1].x-e[i].x);
}
printf("%.2lf\n",ans);
}
return 0;
}