HDU1542 线扫描求矩形面积并

原题的链接为：http://acm.hdu.edu.cn/showproblem.php?pid=1542

这道题是最简单的先扫描的例题，适合初学线扫描的认真仔细理解的一道题。其他的线扫描的题都是在这上面的加深，比如矩形面积并，周长等，思想都是一样的，只是更新部分有的差别而已。说来也就，建树，去重，最后加一条边，计算一次矩形的面积。线扫描的题最主要的就是考虑好状态的更新。

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
#define nn 110

struct node
{
int l,r,v;
double len;
}tree[nn<<3];

struct seg
{
double l,r,h;
int f;//-1表示上边，1表示下边
}s[nn<<1];
double pos[nn<<1];

bool cmp(seg a,seg b)
{
return a.h<b.h;
}
int binary(double c,int low,int top)//二分查找该点的离散值
{
while(low<=top)
{
int mid=(low+top)>>1;
if(pos[mid]==c)     return mid;
else if(pos[mid]>c) top=mid;
else low=mid+1;
}
return -1;
}
void build(int ll,int rr,int rt)//建树
{
tree[rt].l=ll;  tree[rt].r=rr;  tree[rt].v=0;
if(ll==rr)  return;
int mid=(ll+rr)>>1;
build(ll,mid,rt<<1);
build(mid+1,rr,rt<<1|1);
}
void getlen(int rt)//长度更新
{
if(tree[rt].v)//该区间被覆盖
{
tree[rt].len=pos[tree[rt].r+1]-pos[tree[rt].l];
}
else if(tree[rt].l==tree[rt].r)//更新的区间为一个区间
{
tree[rt].len=0;
}
else//更新的区间不为一个点，由他的两个儿子节点来更新他的长度
{
tree[rt].len=tree[rt<<1].len+tree[rt<<1|1].len;
}
}
void updata(int ll,int rr,int c,int rt)//更新的区间查找
{
if(tree[rt].l==ll && tree[rt].r==rr)
{
tree[rt].v+=c;
getlen(rt);
return ;
}
int mid=(tree[rt].l+tree[rt].r)>>1;
if(rr<=mid)
{
updata(ll,rr,c,rt<<1);
}
else if(ll>mid)
{
updata(ll,rr,c,rt<<1|1);
}
else
{
updata(ll,mid,c,rt<<1);
updata(mid+1,rr,c,rt<<1|1);
}
getlen(rt);
}
int main()
{
int n,t=1;
while(scanf("%d",&n)&&n)
{
int num=0;
for(int i=0;i<n;i++)
{
double x1,x2,y1,y2;
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
s[num].l=x1;   s[num].r=x2;    s[num].h=y1;    s[num].f=1;
s[num+1].l=x1;   s[num+1].r=x2;    s[num+1].h=y2;    s[num+1].f=-1;
pos[num]=x1;   pos[num+1]=x2;
num+=2;
}
sort(pos,pos+num);//离散化
sort(s,s+num,cmp);
int m=1;
for(int i=1;i<num;i++)//去重
{
if(pos[i]!=pos[i-1])
pos[m++]=pos[i];
}
build(0,m-1,1);

double ans=0;//记录面积
for(int i=0;i<num;i++)
{

int l=binary(s[i].l,0,m-1);
int r=binary(s[i].r,0,m-1)-1;
printf("%d %d\n",l,r);
updata(l,r,s[i].f,1);
ans+=(s[i+1].h-s[i].h)*tree[1].len;

}
printf("Test case #%d\n",t++);
printf("Total explored area: %.2lf\n\n",ans);
}
return 0;
}