HDU 1255 ( 覆盖的面积 ) 矩形面积并

最新推荐文章于 2021-01-19 17:58:25 发布

原创最新推荐文章于 2021-01-19 17:58:25 发布 · 446 阅读

0 ·

CC 4.0 BY-SA版权

线段树专栏收录该内容

13 篇文章

订阅专栏

本文介绍了一种解决矩形面积并问题的算法，通过线段树数据结构记录覆盖次数，实现了对覆盖两次以上区域面积的高效计算。适用于ACM竞赛等场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1255

矩形面积并，不过题目要求的是覆盖两次以上的面积，这样的话，就在线段树里增加一个变量表示覆盖两次以上的长度，更新的时候注意一下这个变量就可以了

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
const int maxn=2005;
struct line{
    double h,l,r;
    int f;
    line(){}
    line(double a,double b,double c,int d):h(a),l(b),r(c),f(d){}
}lines[maxn];
double xs[maxn];
struct seg{
    double len,dlen;
    int cover;
}sum[maxn<<2];
int linecnt,xcnt;
inline bool cmp(line a,line b){return a.h<b.h;}
inline void addline(double x1,double y1,double x2,double y2){
    xs[++linecnt]=x1;lines[linecnt]=line(y1,x1,x2,1);
    xs[++linecnt]=x2;lines[linecnt]=line(y2,x1,x2,-1);
}
int Search(double k,int l,int r){
    while(l<=r){
        int m=(l+r)>>1;
        if(xs[m]==k)return m;
        if(xs[m]>k)r=m-1;
        else l=m+1;
    }
    return l;
}
inline void pushup(int rt,int l,int r){
    if(sum[rt].cover>=2){
        sum[rt].len=sum[rt].dlen=xs[r+1]-xs[l];
    }
    else if(sum[rt].cover==1){
        sum[rt].len=xs[r+1]-xs[l];
        if(l==r)sum[rt].dlen=0;
        else sum[rt].dlen=sum[rt<<1].len+sum[rt<<1|1].len;
    }
    else{
        if(l==r)sum[rt].len=sum[rt].dlen=0;
        else{
            sum[rt].len=sum[rt<<1].len+sum[rt<<1|1].len;
            sum[rt].dlen=sum[rt<<1].dlen+sum[rt<<1|1].dlen;
        }
    }
}
void update(int L,int R,int c,int l,int r,int rt){
    if(L<=l&&r<=R){
        sum[rt].cover+=c;
        pushup(rt,l,r);
        return;
    }
    int m=(l+r)>>1;
    if(L<=m)update(L,R,c,lson);
    if(R>m)update(L,R,c,rson);
    pushup(rt,l,r);
}
int main(){
//    freopen("D://input.txt","r",stdin);
    int T;scanf("%d",&T);
    while(T--){
        int n;scanf("%d",&n);
        linecnt=xcnt=0;
        while(n--){
            double xx[2],yy[2];
            scanf("%lf%lf%lf%lf",&xx[0],&yy[0],&xx[1],&yy[1]);
            addline(xx[0],yy[0],xx[1],yy[1]);
        }
        sort(xs+1,xs+linecnt+1);sort(lines+1,lines+linecnt+1,cmp);
        int i;memset(sum,0,sizeof(sum));
        for(i=1;i<=linecnt;i++){
            if(xs[i+1]!=xs[i])
                xs[++xcnt]=xs[i];
        }
        double area=0;
        for(i=1;i<linecnt;i++){
            int L=Search(lines[i].l,1,xcnt),R=Search(lines[i].r,1,xcnt)-1;
            if(L<=R){
                update(L,R,lines[i].f,1,xcnt,1);
                area+=sum[1].dlen*(lines[i+1].h-lines[i].h);
            }
        }
        printf("%.2f\n",area);
    }
    return 0;
}