扫描线求面积，求周长。

大佬的讲解，看完豁然开朗https://www.cnblogs.com/KonjakJuruo/p/6024266.html

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;

const int N=1100,INF=(int)1e9;
double z[N][4];
struct point{
    double d;
    int x,y;
}p[2*N];
struct node{
    int x0,x1,d;
    double y;
}a[N];
struct trnode{
    int l,r,lc,rc,cnt;
    double rl,len;
}t[2*N];
double num[N];
int n,tl,pl,al;

int minn(int x,int y){return x<y ? x:y;}
int maxx(int x,int y){return x>y ? x:y;}
bool cmp_d(point x,point y){return x.d<y.d;}
bool cmp_y(node x,node y){return x.y<y.y;}

int bt(int l,int r)
{
    int x=++tl;
    t[x].l=l;t[x].r=r;
    t[x].lc=t[x].rc=0;
    t[x].cnt=0;t[x].len=0;
    t[x].rl=num[r+1]-num[l];
    if(l<r)
    {
        int mid=(l+r)/2;
        t[x].lc=bt(l,mid);
        t[x].rc=bt(mid+1,r);
    }
    return x;
}

void upd(int x)
{
    int lc=t[x].lc,rc=t[x].rc;
    if(t[x].cnt>0) t[x].len=t[x].rl;
    else t[x].len=t[lc].len+t[rc].len;
}

void change(int x,int l,int r,int d)
{
    if(t[x].l==l && t[x].r==r)
    {
        t[x].cnt+=d;
        upd(x);
        return ;
    }
    int lc=t[x].lc,rc=t[x].rc,mid=(t[x].l+t[x].r)/2;
    if(r<=mid) change(lc,l,r,d);
    else if(l>mid) change(rc,l,r,d);
    else
    {
        change(lc,l,mid,d);
        change(rc,mid+1,r,d);
    }
    upd(x);
}

int main()
{
    freopen("a.in","r",stdin);
    int T=0;
    while(1)
    {
        scanf("%d",&n);
        if(n==0) break;
        int x,mx;pl=0;tl=0;al=0;t[0].len=0;
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<=3;j++)
            {
                scanf("%lf",&z[i][j]);
                if(j%2==0) p[++pl].d=z[i][j],p[pl].x=i,p[pl].y=j;
            }
        }
        sort(p+1,p+1+pl,cmp_d);
        mx=0;p[0].d=INF;
        for(int i=1;i<=pl;i++)
        {
            if(p[i].d!=p[i-1].d) mx++,num[mx]=p[i].d;
            z[p[i].x][p[i].y]=mx;
        }
        bt(1,mx);
        for(int i=1;i<=n;i++)
        {
            if(z[i][1]<z[i][3]) swap(z[i][1],z[i][3]);
            if(z[i][0]>z[i][2]) swap(z[i][0],z[i][2]);
            a[++al].x0=z[i][0];a[al].x1=z[i][2];a[al].y=z[i][1];a[al].d=-1;
            a[++al].x0=z[i][0];a[al].x1=z[i][2];a[al].y=z[i][3];a[al].d=1;
        }
        sort(a+1,a+1+al,cmp_y);
        // for(int i=1;i<=al;i++)
            // printf("%d -> %d  %lf %d\n",a[i].x0,a[i].x1,a[i].y,a[i].d);
        double w,h,ans=0;
        change(1,a[1].x0,a[1].x1-1,a[1].d);
        for(int i=2;i<=al;i++)
        {
            h=a[i].y-a[i-1].y;
            w=t[1].len;
            ans+=w*h;
            // printf("w = %lf  h = %lf\n",w,h);
            change(1,a[i].x0,a[i].x1-1,a[i].d);
        }
        printf("Test case #%d\nTotal explored area: %.2lf\n\n",++T,ans);
    }
    return 0;
}

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long LL;
const int N=10*5000,INF=(int)1e9;
struct node{
    int x1,x2,y,d;
}a[N];
struct point{
    int x,y;LL d;
}p[N];
struct trnode{
    int l,r,lc,rc,cnt;
    LL rl,len;
}t[2*N];
int n,mx,pl,al,tl,z[N][4];
LL num[N];

bool cmp_d(point x,point y){return x.d<y.d;}
bool cmp_y(node x,node y)
{
    if(x.y==y.y) return x.d>y.d;//处理重边
    else return x.y<y.y;
}
void inp(int x,int y,int d){p[++pl].x=x;p[pl].y=y;p[pl].d=d;}
void ina(int x1,int x2,int y,int d){a[++al].x1=x1;a[al].x2=x2;a[al].y=y;a[al].d=d;}
int myabs(int x){return x>0 ? x:-x;}
int maxx(int x,int y){return x>y ? x:y;}

void unique()
{
    int now=0;p[0].d=INF;
    sort(p+1,p+1+pl,cmp_d);
    for(int i=1;i<=pl;i++)
    {
        if(p[i].d!=p[i-1].d) now++,num[now]=(LL)p[i].d;
        z[p[i].x][p[i].y]=now;
    }
    mx=now;
}

void make_edge()
{
    al=0;mx=0;
    for(int i=1;i<=n;i++)
    {
        if(z[i][0]>z[i][2]) swap(z[i][0],z[i][2]);
        if(z[i][1]>z[i][3]) swap(z[i][1],z[i][3]);
        ina(z[i][0],z[i][2],z[i][1],1);
        ina(z[i][0],z[i][2],z[i][3],-1);
        mx=maxx(mx,z[i][2]);
    }
    // ina(1,mx-1,INF,-1);
}

int bt(int l,int r)
{
    int x=++tl;
    t[x].l=l;t[x].r=r;
    t[x].lc=t[x].rc=0;
    t[x].cnt=0;t[x].len=0;
    t[x].rl=num[r+1]-num[l];
    if(l<r)
    {
        int mid=(l+r)/2;
        t[x].lc=bt(l,mid);
        t[x].rc=bt(mid+1,r);
    }
    return x;
}

void upd(int x)
{
    int lc=t[x].lc,rc=t[x].rc;
    if(t[x].cnt>=1) t[x].len=t[x].rl;
    else t[x].len=t[lc].len+t[rc].len;
}

void change(int x,int l,int r,int d)
{
    if(t[x].l==l && t[x].r==r)
    {
        t[x].cnt+=d;
        upd(x);
        return;
    }
    int lc=t[x].lc,rc=t[x].rc,mid=(t[x].l+t[x].r)/2;
    if(r<=mid) change(lc,l,r,d);
    else if(l>mid) change(rc,l,r,d);
    else
    {
        change(lc,l,mid,d);
        change(rc,mid+1,r,d);
    }
    upd(x);
}

void output()
{
    for(int i=1;i<=tl;i++)
        printf("l = %d  r = %d  cnt = %d  len  = %lf  rl = %lf \n",t[i].l,t[i].r,t[i].cnt,t[i].len,t[i].rl);
}

int main()
{
    freopen("a.in","r",stdin);
    while(scanf("%d",&n)!=EOF)
    {
        LL pre,ans,now;ans=0;pl=0;
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<=3;j++)
            {
                scanf("%d",&z[i][j]);
                inp(i,j,z[i][j]);
            }
        }
        unique();
        for(int k=0;k<=1;k++)
        {
            if(k==1)
            {
                for(int i=1;i<=n;i++) swap(z[i][0],z[i][1]),swap(z[i][2],z[i][3]);
            }
            make_edge();
            sort(a+1,a+1+al,cmp_y);
            tl=0;bt(1,mx-1);
            pre=0;
            for(int i=1;i<=al;i++)
            {
                change(1,a[i].x1,a[i].x2-1,a[i].d);
                now=t[1].len;
                ans+=myabs(now-pre);
                pre=now;
            }
            // printf("ans = %d\n",ans);
        }
        printf("%lld\n",ans);
    }
    return 0;
}