今天刚开始接触扫描线和线段树的做法,还有点不熟,直接上代码吧
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=2010;
struct node
{
double x;
double y1;
double y2;
int flag;
bool operator <(const node &a)const
{
return x<a.x;
}
}p[maxn];
struct Tree
{
int l;
int r;
int cover;
double lena;
double lenb;
}t[maxn*3];
int n,m;
double py[maxn];
void Build(int l,int r,int index)
{
t[index].l=l;
t[index].r=r;
t[index].cover=0;
t[index].lena=t[index].lenb=0;
if(l+1==r)
return;
int mid=(l+r)>>1;
Build(l,mid,index<<1);
Build(mid,r,index<<1|1);
}
void Callen(int index)
{
if(t[index].cover>1)
t[index].lena=t[index].lenb=py[t[index].r]-py[t[index].l];
else if(t[index].cover==1)
{
t[index].lena=py[t[index].r]-py[t[index].l];
if(t[index].l+1==t[index].r)
t[index].lenb=0;
else
t[index].lenb=t[index<<1].lena+t[index<<1|1].lena;
}
else
{
if(t[index].l+1==t[index].r)
t[index].lena=t[index].lenb=0;
else
{
t[index].lenb=t[index<<1].lenb+t[index<<1|1].lenb;
t[index].lena=t[index<<1].lena+t[index<<1|1].lena;
}
}
}
void Update(int l,int r,int val,int index)
{
if(t[index].l==l&&t[index].r==r)
{
t[index].cover+=val;
Callen(index);
return;
}
int mid=(t[index].l+t[index].r)>>1;
if(r<=mid)
Update(l,r,val,index<<1);
else if(l>=mid)
Update(l,r,val,index<<1|1);
else
{
Update(l,mid,val,index<<1);
Update(mid,r,val,index<<1|1);
}
Callen(index);
}
int find(double val)
{
int l=1,r=m;
while(l<=r)
{
int mid=(l+r)>>1;
if(py[mid]==val)
return mid;
if(py[mid]<val)
l=mid+1;
else
r=mid-1;
}
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
m=1;
for(int i=1;i<=n;i++)
{
double x1,y1,x2,y2;
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
p[m].x=x1;
p[m].y1=y1;
p[m].y2=y2;
p[m].flag=1;
py[m++]=y1;
p[m].x=x2;
p[m].y1=y1;
p[m].y2=y2;
p[m].flag=-1;
py[m++]=y2;
}
m--;
sort(p+1,p+m+1);
sort(py+1,py+m+1);
Build(1,m,1);
int count=2;
double ans=0;
for(int i=1;i<m;i++)
{
int l=find(p[i].y1);
int r=find(p[i].y2);
Update(l,r,p[i].flag,1);
ans+=t[1].lenb*(p[i+1].x-p[i].x);
}
printf("%.2f\n",ans);
}
return 0;
}