分情况讨论
1.cnt>1 : 说明该区间被覆盖两次或以上,那么长度就可以直接计算,就是该区间的长度
剩下的情况就是cnt=1或cnt=0
2.先看叶子节点,因为是叶子没有孩子了,所以被覆盖两次货以上的长度就是0(无论cnt=1或cnt=0都是0,因为是叶子。。。)
3.不是叶子节点 ,且cnt=1.注意这里,cnt=1确切的意义是什么,应该是,可以确定,这个区间被完全覆盖了1次,而有没有被完全覆盖两次或以上则不知道无法确定,那么怎么怎么办了,只要加上t[lch].s + t[rch].s 即,看看左右孩子区间被覆盖了一次或以上的长度,那么叠加在双亲上就是双亲被覆盖两次或以上的长度
3.不是叶子节点,且cnt=0,确切的意义应该是不完全不知道被覆盖的情况(不知道有没有被覆盖,被覆盖了几次,长度是多少都不知道),这种情况,只能由其左右孩子的信息所得
t[lch].ss + t[rch].ss , 即直接将左右孩子给覆盖了两次或以上的长度加起来,这样才能做到不重不漏
#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <algorithm>
#include <math.h>
#include <map>
#include <vector>
#include <string.h>
#include <string>
#include <queue>
#include <set>
using namespace std;
typedef long long ll;
const int N = 1100;
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) {}
}e[N<<1];
bool cmp(Edge a, Edge b)
{
return a.h < b.h;
}
struct Node {
int l, r; //横坐标的区间,是横坐标的下标
int s; //该节点被覆盖的情况(是否完全被覆盖)
double len1; //该区间被覆盖1次及以上的长度
double len2; //该区间被覆盖2次及以上的长度
}q[N<<3];//每个矩形有两个横坐标所以需要2*4N的节点
#define ls i<<1
#define rs i<<1|1
#define m(i) ((q[i].l+q[i].r)>>1)
double x[N<<1];
void build(int i, int l, int r)
{
q[i].r = r, q[i].l = l;
q[i].s = q[i].len1 = q[i].len2 = 0;
if(l == r) return;
int mid = m(i);
build(ls, l, mid);
build(rs, mid+1, r);
}
void pushup(int i)
{
if(q[i].s) //非零,已经被整段覆盖
{
q[i].len1 = x[q[i].r+1] - x[q[i].l]; //左闭右开解决区间缺失问题
}
else if(q[i].l == q[i].r) //这是一个点不是线段
{
q[i].len1 = 0;
}
else
{
q[i].len1 = q[ls].len1 + q[rs].len1;
}
if(q[i].s>1) q[i].len2 = x[q[i].r+1] - x[q[i].l]; //该区间被覆盖2次及以上
else if(q[i].l == q[i].r) q[i].len2 = 0;
else if(q[i].s == 1) q[i].len2 = q[ls].len1 + q[rs].len1;
else q[i].len2 = q[ls].len2 + q[rs].len2;
}
void update(int i, int l, int r, int xx)
{
if(q[i].l == l && q[i].r == r)
{
q[i].s += xx;
pushup(i);
return;
}
int mid = m(i);
if(r <= mid) update(ls, l, r, xx);
else if(l > mid) update(rs, l, r, xx);
else
{
update(ls, l, mid, xx);
update(rs, mid+1, r, xx);
}
pushup(i);
}
int main()
{
int n, T;
scanf("%d", &T);
while(T--)
{
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,cmp);
//离散化坐标
sort(x,x+tot);
//去重
int k = unique(x, x + tot) - x;
build(1, 0, k-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(1, l, r, e[i].f);
ans += (e[i+1].h - e[i].h)*q[1].len2;
}
printf("%.2f\n",ans);
}
return 0;
}