HDU-1255-覆盖的面积

今天刚开始接触扫描线和线段树的做法，还有点不熟，直接上代码吧

代码：

#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;
}