HDU-1255 覆盖的面积

最新推荐文章于 2026-07-14 12:05:14 发布

原创最新推荐文章于 2026-07-14 12:05:14 发布 · 437 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

本文介绍了一种使用线段树、扫描线与离散化技术解决矩形覆盖问题的方法，重点在于计算至少被覆盖两次的区域面积。通过具体实现代码展示了如何通过更新线段树来维护区间覆盖次数并计算目标面积。

线段树+扫描线+离散化

利用扫描线算法求解矩形覆盖的面积，和一般的套路一样

不同的是这道题只统计覆盖过两次及以上的面积，在初始化时cover为-1就可以了（或者其他的小修改）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>

using namespace std;
const int maxn=2000+10;

struct line
{
    double x,y1,y2;
    int f;
    bool operator < (const line &u) const
    {
        return x<u.x;
    }
}edge[maxn];
double len[maxn<<2];
int cover[maxn<<2];
double Y[maxn];

void build(int l,int r,int rt)
{
    len[rt]=0.0;
    cover[rt]=-1;
    if(r-l==1) return;
    int m=(l+r+1)>>1;
    build(l,m,rt<<1);
    build(m,r,rt<<1|1);
}

void PushUp(int rt,int l,int r)
{
    if(cover[rt]>0)
    {
        len[rt]=Y[r]-Y[l];
    }
    else if(r-l==1)
    {
        len[rt]=0.0;
    }
    else
    {
        len[rt]=len[rt<<1]+len[rt<<1|1];
    }
}

void update(double L,double R,int f,int l,int r,int rt)
{
    if(r-l==1)
    {
        cover[rt]+=f;
    }
    else
    {
        int m=(l+r+1)>>1;
        if(L<Y[m]) update(L,R,f,l,m,rt<<1);
        if(R>Y[m]) update(L,R,f,m,r,rt<<1|1);
    }
    PushUp(rt,l,r);
}

int main()
{
    //freopen("/home/zlwang/test.txt","r",stdin);
    int n,T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        double x1,x2,y1,y2;
        int cnt=0;
        for(int i=0;i<n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            Y[i]=y1;
            Y[n+i]=y2;
            edge[cnt].y1=y1;
            edge[cnt].y2=y2;
            edge[cnt].x=x1;
            edge[cnt].f=1;
            cnt++;
            edge[cnt].y1=y1;
            edge[cnt].y2=y2;
            edge[cnt].x=x2;
            edge[cnt].f=-1;
            cnt++;
        }
        sort(edge,edge+2*n);
        sort(Y,Y+2*n);
        cnt=1;
        for(int i=1;i<2*n;i++)
        if(Y[i]!=Y[i-1]) Y[cnt++]=Y[i];
        build(0,cnt-1,1);
        double ans=0.0;
        for(int i=0;i<2*n-1;i++)
        {
            update(edge[i].y1,edge[i].y2,edge[i].f,0,cnt-1,1);
            ans+=len[1]*(double)(edge[i+1].x-edge[i].x);
        }
        printf("%.2f\n",ans);
    }
    return 0;
}