HDU - 1542 - Atlantis （线段树扫描线求矩形并）

#include<bits/stdc++.h>
using namespace std;
typedef double ld;
const int N=207;
int n,m;
struct Line
{
int x1,x2;
int down;
ld y;
bool operator < (const Line & r) const
{
return y<r.y;
}
};
struct Shape
{
ld x1,x2,y1,y2;
Shape(){}
Shape(ld x1,ld y1,ld x2,ld y2) : x1(x1),y1(y1),x2(x2),y2(y2) { }
};
Shape s
;
vector<Line> line;
ld xx
,sum[N<<2];
int c[N<<2];
ld dis(int l,int r) { return xx[r]-xx[l-1]; }
void build(int rt,int l,int r)
{
c[rt]=0;sum[rt]=0;
if(l==r) return ;
int mid=(l+r)>>1;
build(rt<<1,l,mid);
build(rt<<1|1,mid+1,r);
}
inline void push_down(int rt,int l,int r)
{
if(c[rt]!=-1)
{
c[rt<<1]=c[rt<<1|1]=c[rt];
int m=(l+r)>>1;
sum[rt<<1]=(c[rt]?dis(l,m):0);
sum[rt<<1|1]=(c[rt]?dis(m+1,r):0);
}
}
inline void push_up(int rt)
{
c[rt]=(c[rt<<1]==c[rt<<1|1]&&c[rt<<1]!=-1?c[rt<<1]:-1);
sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void update(int rt,int l,int r,int ql,int qr,int v)
{
if(ql<=l&&qr>=r)
{
if(c[rt]!=-1)
{
c[rt]+=v;
sum[rt]=(c[rt]?dis(l,r):0);
return;
}
}
push_down(rt,l,r);
int mid=(l+r)>>1;
if(ql<=mid) update(rt<<1,l,mid,ql,qr,v);
if(qr>mid) update(rt<<1|1,mid+1,r,ql,qr,v);
push_up(rt);
}
int main()
{
int kase=1;
while(~scanf("%d",&n))
{
if(n==0) break;
for(int i=0;i<n;++i)
{
ld x1,x2,y1,y2;
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
s[i]=Shape(x1,y1,x2,y2);
xx[i<<1]=x1;xx[i<<1|1]=x2;
}
sort(xx,xx+n*2);
m=unique(xx,xx+n*2)-xx;
line.clear();
for(int i=0;i<n;++i)
{
int x1id=lower_bound(xx,xx+m,s[i].x1)-xx+1;
int x2id=lower_bound(xx,xx+m,s[i].x2)-xx;
line.push_back({x1id,x2id,1,s[i].y1});
line.push_back({x1id,x2id,-1,s[i].y2});
}
// 线段树上为[l,r]的点 表示[X[l-1],X[r]] 这段区间。
sort(line.begin(),line.end());
ld ans=0;
build(1,1,--m);
update(1,1,m,line[0].x1,line[0].x2,line[0].down);
for(int i=1;i<line.size();++i)
{
ans+=sum[1]*(line[i].y-line[i-1].y);
update(1,1,m,line[i].x1,line[i].x2,line[i].down);
}
printf("Test case #%d\nTotal explored area: %.2f\n\n",kase++,ans);
}
return 0;
}