P5490 【模板】扫描线

原题链接

题解请移步题解区，相当详细

code

#include<bits/stdc++.h>
using namespace std;

#define ll long long

inline void read(ll &x) {
    x = 0;
    ll flag = 1;
    char c = getchar();
    while(c < '0' || c > '9'){
        if(c == '-')flag = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        x = (x << 3) + (x << 1) + (c ^ 48);
        c = getchar();
    }
    x *= flag;
}

inline void write(ll x)
{
    if(x < 0){
        putchar('-');
        x = -x;
    }
    if(x > 9)
        write(x / 10);
    putchar(x % 10 + '0');
}

struct unit
{
    ll y1,y2,x,mark;
}seg[8000005];

struct
{
    ll val;
    ll left;
}tree[8010000];

vector<ll> y(4000005,0);

void pushup(ll node,ll l,ll r)
{
    tree[node].val=tree[node*2].val+tree[node*2+1].val;
}

void change(ll node,ll l,ll r,ll yl,ll yr,ll sit)
{

    if(yl>r||yr<l)return ;
    if(l>=yl&&r<=yr)
    {
        tree[node].left+=sit;
        if(tree[node].left)
        {
            tree[node].val=y[r+1]-y[l];//val代表当前区间生效的长度有多少
        }
        else
        {
            tree[node].val=0;
            pushup(node,l,r);//为什么可以直接不更新子节点后再累加呢？因为子节点当前线段没去过，如果往下更新会错误地改变值。同时子区间已经更新过了
            //说明去掉当前线段后恰好没有了这个区间
        }
        return;
    }

    //有交叉
    ll mid=(l+r)>>1;
    if(mid>=yl)change(node<<1,l,mid,yl,yr,sit);
    if((mid+1)<=yr)change((node<<1)|1,mid+1,r,yl,yr,sit);
    if(!tree[node].left)pushup(node,l,r);
}

bool cmp(unit a,unit b)
{
    return a.x<b.x;
}
int main()
{
    ll n;
    read(n);
    for(ll i=1;i<=n;i++)
    {
        ll x1,x2,y1,y2;
        read(x1); read(y1); read(x2); read(y2);

        if(y1>y2)swap(y1,y2);//y2代表上边界
        if(x1>x2)swap(x1,x2);//x2代表右边界

        y[(i<<1)-1]=y1;//离散化记录
        y[i<<1]=y2;

        seg[(i<<1)-1]={y1,y2,x1,1};
        seg[i<<1]={y1,y2,x2,-1};//记录矩形线段对竖区间的修改
    }

    n<<=1;
    sort(y.begin()+1,y.begin()+n+1);
    ll len=unique(y.begin()+1,y.begin()+1+n)-y.begin()-1;//代表有多少个y点

    sort(seg+1,seg+1+n,cmp);

    ll ans=0;
    for(ll i=1;i<n;i++)
    {
        ll y1=lower_bound(y.begin()+1,y.begin()+len+1,seg[i].y1)-y.begin();
        ll y2=lower_bound(y.begin()+1,y.begin()+len+1,seg[i].y2)-y.begin();//当前线段上下边界的离散值

        change(1,1,len-1,y1,y2-1,seg[i].mark);//减一是因为点与区间的对应关系
        if(seg[i+1].x-seg[i].x>0)  ans+=tree[1].val*(seg[i+1].x-seg[i].x);//就像是公交车，等位于这条竖线上的线段修改全部完成之后才会进入下一步
    }

    write(ans);
    putchar('\n');

    return 0;
}