扫描线算法

洛谷P5490

参考

• https://oi-wiki.org/geometry/scanning/

题解：

• 设想有一条直线从下到上扫描整个图形，如图所示，我们把图形分割成若干个矩形，每次计算一个颜色的矩形面积。
• 对输入矩形的y坐标进行升序排序，每扫描一条边，其下一条边的位置就知道了，即高是固定的，唯一要求的就是宽，知道高和宽就可以求解面积。对于每一条扫描线我们用flag标记，矩形的下边flag=1，上边flag=-1。
• 如何求宽？可以将x轴按照x坐标分割成若干段，每一段做一个标记c。每当从下到上扫描到一个边，每一段都加上一个flag增量，如果c＞0表示这一段有效并求解宽的总长度，在结合下一条边的位置求总面积。可以用线段树维护x轴的每一段。
• 线段树维护区间有一种特殊情况需要处理，比如区间为[1,6]，需要求区间[1,4]的长度，[1,6]被分成[1,3]和[4,6]，显然我们只能求出[1,3]的长度，而[3,4]却无法求得。对于这种情况，我们只需对右结点的编号r加1，得len=x[r+1]-x[l]。相应的在查询的时候r要减1。
• 
• 
• 
• 
• 
• 
• 
•  

代码

#include <bits/stdc++.h>
using namespace std;
int const N = 1e6 + 10;
typedef long long ll;
int n;
ll x[N<<1];
struct Edge{
    ll y,x1,x2;
    int flag;
    Edge(){};
    Edge(ll a,ll b,ll c,int d):y(a),x1(b),x2(c),flag(d){};
    bool operator < (const Edge& e)const{
        return y < e.y;
    }
}e[N<<1];
struct Node{
    int id,l,r,cover;
    ll len;
}node[N<<2];
void build(int id,int l,int r){
    node[id].l = l, node[id].r = r;
    if(l == r)  return;
    else{
        int mid = (l + r) >> 1;
        build(id << 1,l,mid);
        build(id << 1 | 1,mid + 1,r);
    }
}
void query(int id){ //询问整个区间cover＞0的区间长度
    int l = node[id].l, r = node[id].r;
    if(node[id].cover)
        node[id].len = x[node[id].r + 1] - x[node[id].l];
    else
        node[id].len = node[id<<1].len + node[id<<1|1].len;
}
void update(int id,int L,int R,int flag){   //自顶向下
    int l = node[id].l, r = node[id].r;
    if(L <= l && r <= R){
        node[id].cover += flag;
    }else{
        int mid = (l + r) >> 1;
        if(L <= mid)    update(id << 1,L,R,flag);
        if(mid < R)     update(id << 1 | 1,L,R,flag);
    }
    query(id);  //自底向上
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        ll x1,y1,x2,y2;
        scanf("%lld%lld%lld%lld",&x1,&y1,&x2,&y2);
        x[i] = x1,  x[i+n] = x2;
        e[i] = Edge(y1,x1,x2,1);
        e[i+n] = Edge(y2,x1,x2,-1);
    }
    sort(e + 1,e + 1 + 2 * n);  //对y轴进行排序
    sort(x + 1,x + 1 + 2 * n);
    int m = unique(x + 1,x + 1 + 2 * n) - (x + 1);  //离散化，去掉重复的x点
    build(1,1,m);
    ll ans = 0;
    for(int i=1;i<=2*n-1;i++){  //扫描y轴,每个点维护一个区间
        int l = lower_bound(x + 1,x + 1 + m,e[i].x1) - x;
        int r = lower_bound(x + 1,x + 1 + m,e[i].x2) - x;
        update(1,l,r - 1,e[i].flag);
        ans += node[1].len * (e[i+1].y - e[i].y);
    }
    cout<<ans<<endl;
    return 0;
}