poj1177（扫描线+线段树）

题目：有n个可以相交的矩形在二维平面中，求这些矩形的总周长。
题解：转大佬，这篇题解讲得比较清楚，之前只会做面积，求周长也差不多的

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
#include<vector>
#define ll long long
#define ull unsigned long long
#define ld long double
#define rep1(i,n) for(int i=1;i<=n;i++)
#define rep0(i,n) for(int i=0;i<n;i++)
#define rep(i,l,r) for(int i=l;i<=r;i++)
#define fre freopen("in.txt","r",stdin)
using namespace std;
const ll mod=10007;
const int maxn=1e5+7,INF=0x3f3f3f3f;
#define lson p*2
#define rson p*2+1
int n;

struct data{
    int x,l,r,val;
}a[maxn];

struct S
{
    int l,r;
    int len,ls,rs,num,cov;
}t[maxn*4];
//cov覆盖次数,num连续段数,len被覆盖总长
void build(int p,int l,int r){
    t[p].l=l,t[p].r=r;
//    printf("%d %d\n",l,r);
    if(l==r)return;
    int mid=(l+r)>>1;//这边一定不能使用/2，会死循环的
    build(lson,l,mid);
    build(rson,mid+1,r);
}

void spread(int p){
    int l=t[p].l,r=t[p].r;
    if(t[p].cov){
        t[p].ls=t[p].rs=1;
        t[p].num=2;
        t[p].len=r-l+1;
    }
    else if(l==r)t[p].ls=t[p].rs=t[p].len=t[p].num=0;
    else{
        t[p].num=t[lson].num+t[rson].num;
        t[p].len=t[lson].len+t[rson].len;
        t[p].ls=t[lson].ls;t[p].rs=t[rson].rs;
        if(t[lson].rs&&t[rson].ls)t[p].num-=2;
    }
}

void change(int p,int l,int r,int val){
    if(t[p].l>=l&&t[p].r<=r){
        t[p].cov+=val;
        spread(p);
        return;
    }
//    printf("%d %d\n",t[p].l,t[p].r);
    int mid=(t[p].l+t[p].r)>>1;
    if(l<=mid)change(lson,l,r,val);
    if(r>mid)change(rson,l,r,val);
    spread(p);
}

bool cmp(data a,data b){
    return a.x<b.x;
}

int main()
{
    scanf("%d",&n);
    int m=0,l=INF,r=-INF;
    for(int x1,x2,y1,y2,i=1;i<=n;i++){
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        l=min(l,y1);r=max(r,y2);
        a[++m]=(data){x1,y1,y2,1};
        a[++m]=(data){x2,y1,y2,-1};
    }
    n=m;
    build(1,l,r-1);
    sort(a+1,a+1+n,cmp);
    int ans=0;
    for(int i=1;i<=n;i++){
        int tmp=t[1].len;
        if(i!=1)ans+=t[1].num*(a[i].x-a[i-1].x);
        change(1,a[i].l,a[i].r-1,a[i].val);
        ans+=abs(t[1].len-tmp);
    }
    printf("%d\n",ans);
    return 0;
}