
2 5 1 1 4 2 1 3 3 7 2 1.5 5 4.5 3.5 1.25 7.5 4 6 3 10 7 3 0 0 1 1 1 0 2 1 2 0 3 1
7.63 0.00
#include<algorithm>
#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
const int INF=0x3f3f3f3f;
const int maxn=2510;
#define lson L,mid,ls
#define rson mid+1,R,rs
int n,m;
int cov[maxn<<2];
double len[maxn<<2],relen[maxn<<2],H[2500];
struct node
{
double x1,x2,h;
int v;
node(double a=0,double b=0,double c=0,int d=0):x1(a),x2(b),h(c),v(d){}
} seg[maxn];
bool cmp(node a,node b)
{
return a.h<b.h;
}
void init()
{
sort(H,H+m);
m=unique(H,H+m)-H;
}
int Hash(double x)
{
return lower_bound(H,H+m,x)-H;
}
void PushUp(int L,int R,int rt)
{
if(cov[rt]>1)//有标记肯定整块覆盖了
len[rt]=relen[rt]=H[R+1]-H[L];
else if(cov[rt]==1)
{
len[rt]=H[R+1]-H[L],relen[rt]=0;
if(L!=R)
relen[rt]=len[rt<<1]+len[rt<<1|1];
}
else if(L==R)//没有左右儿子了。
len[rt]=relen[rt]=0;
else//没有整块覆盖但是被部分覆盖了
len[rt]=len[rt<<1]+len[rt<<1|1],relen[rt]=relen[rt<<1]+relen[rt<<1|1];
}
void update(int L,int R,int rt,int l,int r,int d)
{
if(l<=L&&R<=r)
{
cov[rt]+=d;
PushUp(L,R,rt);
return;
}
int mid=(L+R)>>1,ls=rt<<1,rs=ls|1;
if(l<=mid)
update(lson,l,r,d);
if(r>mid)
update(rson,l,r,d);
PushUp(L,R,rt);
}
int main()
{
int t,i,ptr;
double x1,x2,y1,y2,ans;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(i=ptr=0;i<n;i++)
{
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
H[ptr]=x1;
seg[ptr++]=node(x1,x2,y1,1);
H[ptr]=x2;
seg[ptr++]=node(x1,x2,y2,-1);
}
m=ptr,ans=0;
init();
sort(seg,seg+ptr,cmp);
memset(len,0,sizeof len);
memset(relen,0,sizeof relen);
memset(cov,0,sizeof cov);
for(i=0,ptr--;i<ptr;i++)
{
update(0,m-1,1,Hash(seg[i].x1),Hash(seg[i].x2)-1,seg[i].v);//m个结点m-1条线段
ans+=(seg[i+1].h-seg[i].h)*relen[1];
}
printf("%.2lf\n",ans);
}
return 0;
}
hdu 1255 覆盖的面积(线段树&扫描线&重复面积),布布扣,bubuko.com
原文地址:http://blog.csdn.net/bossup/article/details/36190671