hdu1828线段树扫描线求周长

思想和扫描线求面积差不多

#include<iostream>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<cmath>
using namespace std;

struct node
{
    int l;
    int r;
    int len;   //该区间可与下一个将要插入的线段组成并面积的长度
    int cover; //该区间覆盖了几根线段
    int num;   //记录该区间由几个子区间组成了可与下一个将要插入的线段组成并面积,即可并子区间的数目
    int r_cover; //记录该区间右节点是否被可并区间覆盖
    int l_cover;
};

node tree[100000];
int n;
int yy[10004],len;

struct Line
{
    int l;
    int r;
    int x;
    int cover;
};

Line line[10004];

int cmp(Line a,Line b)
{
    return a.x<b.x;
}

int find(int x)
{
    int l=0,r=len,mid;
    while(l<=r)
    {
        mid=(l+r)/2;
        if(yy[mid]==x)
            return mid;
        if(yy[mid]<x)
            l=mid+1;
        else
            r=mid-1;
    }
    return l;
}

void build(int i,int l,int r)
{
    tree[i].l=l;
    tree[i].r=r;
    tree[i].cover=tree[i].l_cover=tree[i].r_cover=tree[i].len=tree[i].num=0;
    if(l+1==r)
        return;
    int mid=(l+r)/2;
    build(2*i,l,mid);
    build(2*i+1,mid,r);
}

void cal(int i)
{
    if(tree[i].cover) //整个被覆盖
    {
        tree[i].len=yy[tree[i].r]-yy[tree[i].l]; //可用长度为整个区间
        tree[i].l_cover=tree[i].r_cover=1; //左右节点都被覆盖了
        tree[i].num=1;  //由一个区间组成,即该区间
    }
    else if(tree[i].l+1==tree[i].r) //叶子区间
    {
        tree[i].len=0;
        tree[i].l_cover=tree[i].r_cover=0;
        tree[i].num=0;
    }
    else
    {
        tree[i].len=tree[2*i].len+tree[2*i+1].len; //dp
        tree[i].l_cover=tree[2*i].l_cover; //左节点是否覆盖,取决于左子区间的左节点是否被覆盖
        tree[i].r_cover=tree[2*i+1].r_cover; //同理
        tree[i].num=tree[2*i].num+tree[2*i+1].num-tree[2*i].r_cover*tree[2*i+1].l_cover; //该线段可用长度的区间组成数等于左右子区间的组成数
    }                                                                                    //但是,当左右子区间是连续的时候,结果要减一
}

void update(int i,int l,int r,int w)
{
    if(tree[i].l>r || tree[i].r<l)
        return;
    if(tree[i].l>=l && tree[i].r<=r)
    {
        tree[i].cover+=w;
        cal(i);
        return;
    }
    update(2*i,l,r,w);
    update(2*i+1,l,r,w);
    cal(i);
}

int main()
{
    int i,x1,x2,y1,y2,m,a,b;
    while(scanf("%d",&n)!=EOF)
    {
        m=0;
        for(i=0;i<n;i++)
        {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            yy[m]=y1;
            line[m].cover=1;
            line[m].x=x1;
            line[m].l=y1;
            line[m++].r=y2;

            yy[m]=y2;
            line[m].cover=-1;
            line[m].x=x2;
            line[m].l=y1;
            line[m++].r=y2;
        }
        sort(yy,yy+m);
        sort(line,line+m,cmp);
        len=1;
        for(i=1;i<m;i++)//quchong
        {
            if(yy[i-1]!=yy[i])
                yy[len++]=yy[i];
        }
        len--;
        build(1,0,len);
        int ans=0,pre=0;
        for(i=0;i<m;i++)
        {
            a=find(line[i].l);
            b=find(line[i].r);
            update(1,a,b,line[i].cover);
            ans+=abs((tree[1].len-pre)); //加上y坐标上的周长
            if(i==m-1)
            break;
            pre=tree[1].len;
            ans+=tree[1].num*(line[i+1].x-line[i].x)*2;//加上x坐标上的周长,因为一个y边连着两个x边,所以乘二
        }
        printf("%d\n",ans);
    }
    return 0;
}
    原文作者:B树
    原文地址: https://blog.csdn.net/aonaigayiximasi/article/details/50958163
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞

发表评论

电子邮件地址不会被公开。 必填项已用*标注