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

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

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#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<
        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;
}