hdu 1542/1255 Atlantis/覆盖的面积

Posted 2021-02-03 fengxunling

两道扫描线+线段树的入门题。
基本没有什么区别，前者是模板，后者因为是求覆盖次数至少在两次以上的，这个同样是具有并集性质的，所以把cover的判断条件更改一下就可以了qwq。

hdu1542 代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 200010
using namespace std;
int n,c,cnt;
double y[MAXN];
struct Node{
    double x,l,r;
    int cover;
    bool flag;
}node[MAXN];
struct Line
{
    double x,y_up,y_down;
    int flag;
}line[MAXN];
inline bool cmp(struct Line x,struct Line y){return x.x<y.x;}
inline int ls(int x){return x<<1;}
inline int rs(int x){return x<<1|1;}
inline void build(int x,int l,int r)
{
    node[x].l=y[l],node[x].r=y[r],node[x].x=-1,node[x].flag=false,node[x].cover=0;
    if(l+1==r){node[x].flag=true; return;}
    int mid=(l+r)>>1;
    build(ls(x),l,mid);
    build(rs(x),mid,r);
}
inline double q_update(int x,double pos,double l,double r,int flag)
{
    if(l>=node[x].r||r<=node[x].l) return 0;
    if(node[x].flag)
    {
        if(node[x].cover<=0) 
        {
            node[x].x=pos;
            node[x].cover+=flag;
            return 0;
        }
        double pre=node[x].x;
        double ans=(pos-pre)*(node[x].r-node[x].l);
        node[x].x=pos;
        node[x].cover+=flag;
        return ans;
    }
    return q_update(ls(x),pos,l,r,flag)+q_update(rs(x),pos,l,r,flag);
}
int main()
{
    scanf("%d",&n);
    while(n!=0)
    {
        double x1,x2,y1,y2;
        cnt=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            line[++cnt].x=x1,line[cnt].y_down=y1,line[cnt].y_up=y2,line[cnt].flag=1,y[cnt]=y1;
            line[++cnt].x=x2,line[cnt].y_down=y1,line[cnt].y_up=y2,line[cnt].flag=-1,y[cnt]=y2;
        }
        sort(&line[1],&line[cnt+1],cmp);
        sort(&y[1],&y[cnt+1]);
        build(1,1,cnt);
        double ans=0;
        for(int i=1;i<=cnt;i++)
            ans+=q_update(1,line[i].x,line[i].y_down,line[i].y_up,line[i].flag);
        printf("Test case #%d
Total explored area: %.2lf

",++c,ans);
        scanf("%d",&n);
    }
    return 0;
}

hdu1255 代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 200010
using namespace std;
int n,c,cnt,t;
double y[MAXN];
struct Node{
    double x,l,r;
    int cover;
    bool flag;
}node[MAXN];
struct Line
{
    double x,y_up,y_down;
    int flag;
}line[MAXN];
inline bool cmp(struct Line x,struct Line y){return x.x<y.x;}
inline int ls(int x){return x<<1;}
inline int rs(int x){return x<<1|1;}
inline void build(int x,int l,int r)
{
    node[x].l=y[l],node[x].r=y[r],node[x].x=-1,node[x].flag=false,node[x].cover=0;
    if(l+1==r){node[x].flag=true; return;}
    int mid=(l+r)>>1;
    build(ls(x),l,mid);
    build(rs(x),mid,r);
}
inline double q_update(int x,double pos,double l,double r,int flag)
{
    if(l>=node[x].r||r<=node[x].l) return 0;
    if(node[x].flag)
    {
        if(node[x].cover<=1) 
        {
            node[x].x=pos;
            node[x].cover+=flag;
            return 0;
        }
        double pre=node[x].x;
        double ans=(pos-pre)*(node[x].r-node[x].l);
        node[x].x=pos;
        node[x].cover+=flag;
        return ans;
    }
    return q_update(ls(x),pos,l,r,flag)+q_update(rs(x),pos,l,r,flag);
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        double x1,x2,y1,y2;
        cnt=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            line[++cnt].x=x1,line[cnt].y_down=y1,line[cnt].y_up=y2,line[cnt].flag=1,y[cnt]=y1;
            line[++cnt].x=x2,line[cnt].y_down=y1,line[cnt].y_up=y2,line[cnt].flag=-1,y[cnt]=y2;
        }
        sort(&line[1],&line[cnt+1],cmp);
        sort(&y[1],&y[cnt+1]);
        build(1,1,cnt);
        double ans=0;
        for(int i=1;i<=cnt;i++)
            ans+=q_update(1,line[i].x,line[i].y_down,line[i].y_up,line[i].flag);
        printf("%.2lf
",ans);
    }
    return 0;
}