hdu 1255 覆盖的面积

题目链接如下：

Problem - 1255

扫描线加离散化加线段树。

扫描线就是用来求面积的。感谢大佬：

HDU 1255（线段树&离散化&扫描线）详细解析_龙卡卡卡的博客-CSDN博客_hdu - 1255

离散化是因为这个有小数，不太好维护，于是直接映射成整数处理。

这里把一个区间当成叶子节点。核心是push_up操作。

void push_up(int l,int r,int rt){
    if (tree[rt].tag>=1){
        tree[rt].len1=x[r]-x[l];
    } else if (l+1==r){
        tree[rt].len1=0;
    } else{
        tree[rt].len1=tree[ls(rt)].len1+tree[rs(rt)].len1;
    }
    if (tree[rt].tag>=2){
        tree[rt].len2=x[r]-x[l];
    } else if (tree[rt].tag==1){
        tree[rt].len2=tree[ls(rt)].len1+tree[rs(rt)].len1;
    } else if (l+1==r){
        tree[rt].len2=0;
    } else{
        tree[rt].len2=tree[ls(rt)].len2+tree[rs(rt)].len2;
    }
}

这里的关键是第二个判断串的第二个else if里边为啥是两个len1，很简单，len1代表的是tag是1的节点的区间和，如果当前节点tag是1，那么len2自然而然就会包括子节点的len1和，因为两个加起来就是2。

这里用0开始建立完全二叉树，注意一下左右区间。

代码如下所示：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<ctime>
#include<vector>
#include<algorithm>
#include<cstring>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<string>
#include<cmath>
#include<climits>

using namespace std;
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define ls(rt) rt<<1|1
#define rs(rt) (rt<<1)+2
const int MAXN=1001;
struct edge{
    double l,r;
    double height;
    int f;
    bool operator<(const edge& e)const{
        return this->height<e.height;
    }
}g[MAXN<<1];

struct node{
    int tag;
    double len1,len2;
}tree[MAXN<<3];

double x[MAXN<<2];//离散化横坐标

int size;

// 用一个完全二叉树来建树
void init(int n){
    size=1;
    while (size<n) size<<=1;
    for (int i = 0; i < 2*size-1; ++i) {
        tree[i].tag=tree[i].len1=tree[i].len2=0;
    }
}

void push_up(int l,int r,int rt){
    if (tree[rt].tag>=1){
        tree[rt].len1=x[r]-x[l];
    } else if (l+1==r){
        tree[rt].len1=0;
    } else{
        tree[rt].len1=tree[ls(rt)].len1+tree[rs(rt)].len1;
    }
    if (tree[rt].tag>=2){
        tree[rt].len2=x[r]-x[l];
    } else if (tree[rt].tag==1){
        tree[rt].len2=tree[ls(rt)].len1+tree[rs(rt)].len1;
    } else if (l+1==r){
        tree[rt].len2=0;
    } else{
        tree[rt].len2=tree[ls(rt)].len2+tree[rs(rt)].len2;
    }
}

void update(int a,int b,int val,int l,int r,int rt){
    if (b<=l || a>=r) return;
    if (a<=l && r<=b){
        tree[rt].tag+=val;
        push_up(l,r,rt);
        return;
    }
    int mid=(l+r)>>1;
//    if (a<=mid){
//        update(a,b,val,lson);
//    }
//    if (b>mid){
//        update(a,b,val,rson);
//    }
    update(a,b,val,l,mid, ls(rt));
    update(a,b,val,mid,r,rs(rt));
    push_up(l,r,rt);
}

int main(){
    int t;
    cin>>t;
    while (t--){
        int n;
        cin>>n;
        double px,py,px1,py1;
        for (int i = 0; i < n; ++i) {
            cin>>px>>py>>px1>>py1;
            g[i<<1].f=-1,g[i<<1].height=py1,g[i<<1].l=px,g[i<<1].r=px1;
            g[i<<1|1].f=1,g[i<<1|1].height=py,g[i<<1|1].l=px,g[i<<1|1].r=px1;
            x[i<<1]=px,x[i<<1|1]=px1;
        }
        int xsize=(n<<1);
        sort(g,g+xsize);
        sort(x,x+xsize);
        xsize=unique(x,x+xsize)-x;
        init(xsize-1);
        for (int i = xsize; i <= size; ++i) {
            x[i]=x[xsize-1];
        }
        double last_height=g[0].height;
        double ans=0;
        for (int i = 0; i < (n<<1); ++i) {
            double h=g[i].height-last_height;
            last_height=g[i].height;
            ans+=h*(tree[0].len2);
            int a= lower_bound(x,x+xsize,g[i].l)-x;
            int b= lower_bound(x,x+xsize,g[i].r)-x;
            update(a,b,g[i].f,0,size,0);
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}