一个连接,还不错:https://www.cnblogs.com/scau20110726/archive/2013/04/12/3016765.html
hdu1542:http://acm.hdu.edu.cn/showproblem.php?pid=1542
题意:给你每个长方形的左下点的坐标和右上点的坐标,求他们的面积和。
hdu 1255:http://acm.hdu.edu.cn/showproblem.php?pid=1255
题意:给你每个长方形的左下点的坐标和右上点的坐标,求出被这些矩形覆盖过至少两次的区域的面积
与上一题议题差别在第一题求总的面积,第二题求覆盖两次以上的面积。就len2表示下。
1.q[rt].s>1 : 说明该区间被覆盖两次或以上,那么长度就可以直接计算,就是该区间的长度
剩下的情况就是q[rt].s=1或q[rt].s=0
2.先看叶子节点,因为是叶子没有孩子了,所以被覆盖两次货以上的长度就是0(无论q[rt].s=1或q[rt].s=0都是0,因为是叶子。。。)
3.不是叶子节点 ,且q[rt].s=1.注意这里,q[rt].s=1确切的意义是什么,应该是,可以确定,这个区间被完全覆盖了1次,而有没有被完全覆盖两次或以上则不知道无法确定,那么怎么怎么办了,只要加上q[ls].len1+q[rs].len1 即,看看左右孩子区间被覆盖了一次或以上的长度,那么叠加在双亲上就是双亲被覆盖两次或以上的长度
3.不是叶子节点,且q[rt].s=0,确切的意义应该是不完全不知道被覆盖的情况(不知道有没有被覆盖,被覆盖了几次,长度是多少都不知道),这种情况,只能由其左右孩子的信息所得
q[ls].len2+q[rs].len2 , 即直接将左右孩子给覆盖了两次或以上的长度加起来,这样才能做到不重不漏
#include <cstdio>
#include <queue>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <math.h>
#include <map>
#include <set>
#include <bitset>
#include <stack>
#define ull unsigned long long
#define ls rt<<1
#define rs rt<<1|1
#define mems(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long ll;
const ll mod=1000000007;
const int N=2210;
const double pi=acos(-1);
const int inf=0x3f3f3f3f;
const int M=50000+5;
struct edge
{
double l,r;//这条线的左右端点的横坐标
double h;//线段纵坐标
int f;//矩形的上边 1 ,下边 -1
edge() {}
edge(double l,double r,double h,int f):l(l),r(r),h(h),f(f) {}
bool operator <(const edge &a) const
{
return h<a.h;
}
} e[N<<1];
bool cmp(edge a,edge b)
{
return a.h<b.h;
}
struct node
{
int l,r;//横坐标区间
int s;//该节点被覆盖的情况
double len1;//区间被覆盖一次以上的长度
double len2;//区间被覆盖两次的长度
} q[N<<3];//每个矩形有两个横坐标,所以是2*4N的节点
//double sum[N<<2];
//int cnt[n<<2];
double x[N<<2];
int binary(double c,int low,int top)//二分查找该点的离散值
{
while(low<=top)
{
int mid=(low+top)>>1;
if(x[mid]==c) return mid;
else if(x[mid]>c) top=mid;
else low=mid+1;
}
return -1;
}
void build(int l,int r,int rt)
{
q[rt].l=l;
q[rt].r=r;
q[rt].s=q[rt].len1=q[rt].len2=0;
if(l==r)
return;
int mid=(l+r)>>1;
build(l,mid,ls);
build(mid+1,r,rs);
}
void pushup(int rt)
{
if(q[rt].s)//非零,已经被整段覆盖
{
q[rt].len1=x[q[rt].r+1]-x[q[rt].l];
}
else if(q[rt].l==q[rt].r)//这是一个点
{
q[rt].len1=0;
}
else
{
q[rt].len1=q[ls].len1+q[rs].len1;
}
// if(q[rt].s>1)//该区间被覆盖两次以上
// q[rt].len2=x[q[rt].r+1]-x[q[rt].l];
// else if(q[rt].l==q[rt].r)
// q[rt].len2=0;
// else if(q[rt].s==1)
// q[rt].len2=q[ls].len1+q[rs].len1;
// else
// q[rt].len2=q[ls].len2+q[rs].len2;
}
void update(int l,int r,int c,int rt)
{
if(q[rt].l==l&&r==q[rt].r)
{
q[rt].s+=c;
pushup(rt);
return;
}
int mid=(q[rt].l+q[rt].r)>>1;
if(r<=mid)
update(l,r,c,ls);
else if(l>mid)
update(l,r,c,rs);
else
{
update(l,mid,c,ls);
update(mid+1,r,c,rs);
}
pushup(rt);
}
int main()
{
int n;
int p=1;
while(scanf("%d",&n))
{
int tot=0;
for(int i=0; i<n; i++)
{
double x1,x2,y1,y2;
scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
e[tot]=edge(x1,x2,y1,1);
e[tot+1]=edge(x1,x2,y2,-1);
x[tot]=x1;
x[tot+1]=x2;
tot+=2;
}
sort(e,e+tot);
//离散化坐标
sort(x,x+tot);
int k=unique(x,x+tot)-x;//去重
build(0,k-1,1);
double ans=0.0;
for(int i=0;i<tot;i++)
{
int l=lower_bound(x,x+k,e[i].l)-x;
int r=lower_bound(x,x+k,e[i].r)-x-1;
update(l,r,e[i].f,1);
ans+=(e[i+1].h-e[i].h)*q[1].len1;
}
printf("Test case #%d\n",p++);
printf("Total explored area: %.2f\n\n",ans);
}
return 0;
}