Picture POJ - 1177(线段树+扫描线解决周长问题)

比解决面积的时候多了一个最左端点是否覆盖标记和最右端点覆盖标记,num是区间内线段段数,pushup注意线段合并的时候,当左子树最右端覆盖,右子树最左端覆盖,合并计算总线段段数时要减一,因为中间两个线段合并为一段

#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<stack>
#include<set>
#include<ctime>
#include<algorithm>
#define eps 1e-14
#define pi acos(-1)
#define ll long long
#define RD T*(rand()*2-RAND_MAX)
#define Drand (long double)rand()/RAND_MAX
#define INF 0x3f3f3f3f
using namespace std;
const int maxn=1e5+100;

int x[maxn],tot=1;

struct tree
{
    int len;
    int cov;
    int lcov;
    int rcov;
    int num;
}ans[maxn<<3];

struct line
{
    int l,r,h,f;
    line(){}
    line(int l,int r,int h,int f):l(l),r(r),h(h),f(f){}
};

int cmp(line a,line b)
{
    return a.h<b.h;
}

void pushup(int l,int r,int rt)
{
    if(ans[rt].cov>0){
        ans[rt].len=x[r+1]-x[l];
        ans[rt].lcov=ans[rt].rcov=1;
        ans[rt].num=1;
    }
    else{
        ans[rt].num=ans[rt<<1].num+ans[rt<<1|1].num;
        if(ans[rt<<1].rcov && ans[rt<<1|1].lcov)ans[rt].num--;
        ans[rt].len=ans[rt<<1].len+ans[rt<<1|1].len;
        ans[rt].lcov=ans[rt<<1].lcov;
        ans[rt].rcov=ans[rt<<1|1].rcov;
    }
}

void build(int l,int r,int rt)
{
    ans[rt].cov=0;
    ans[rt].lcov=0;
    ans[rt].rcov=0;
    ans[rt].len=0;
    ans[rt].num=0;
    if(l==r)return;
    int mid=(r+l)>>1;
    build(l,mid,rt<<1);
    build(mid+1,r,rt<<1|1);
}

void Update(int L,int R,int l,int r,int rt,int c)
{
    if(L<=l && r<=R)
    {
        ans[rt].cov+=c;
        pushup(l,r,rt);
        return;
    }
    int mid=(l+r)>>1;
    if(L<=mid)Update(L,R,l,mid,rt<<1,c);
    if(R>mid)Update(L,R,mid+1,r,rt<<1|1,c);
    pushup(l,r,rt);
}



int main()
{
//    freopen("in1.txt","r",stdin);
    int n;
    scanf("%d",&n);
    vector<line>v;
    for(int i=0;i<n;i++){
        int a,b,c,d;
        scanf("%d%d%d%d",&a,&b,&c,&d);
        v.push_back(line(a,c,b,1));
        v.push_back(line(a,c,d,-1));
        x[tot++]=a;
        x[tot++]=c;
    }
    sort(v.begin(),v.end(),cmp);
    sort(x+1,x+tot);
    tot=unique(x+1,x+tot)-x;
    build(1,tot-2,1);
    int all=0,last=0;
    for(int i=0;i<v.size()-1;i++){
        int l=lower_bound(x+1,x+tot,v[i].l)-x;
        int r=lower_bound(x+1,x+tot,v[i].r)-x-1;
        Update(l,r,1,tot-2,1,v[i].f);
//        cout<<ans[1].len<<" "<<last<<" "<<ans[1].num<<" "<<(v[i+1].h-v[i].h)<<endl;
        all+=abs(last-ans[1].len)+2*ans[1].num*(v[i+1].h-v[i].h);
        last=ans[1].len;
    }
    cout<<all+ans[1].len<<endl;
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值