数据结构——扫描线

扫描线算法基础

扫描线刚刚学完，自己总结总结巩固一下，如果有不足或不对的地方，请大佬们及时指出。
1. 知识点
首先，学习扫描线，我们要有一定线段树和离散化的知识
2. 作用
它可以帮助我们解决一些重叠矩形的面积和周长的问题
3. 具体思路
我的做题思路一般是先根据样例画出图形和树的结构，以便观察更新规律，找到“pushup”子函数的写法
下面，讲一下扫描线的具体实现方法：

如图，我们从下往上便利，第1条边长为20-10=10，并且第1条边是入边，区间[10,20]+1；直到到达第2条边，我们得到第一部分的宽15-10=5，所以第一部分的S=510=50;

然后，第二条边长为25-10=15，且发现第2条边也是入边，区间[10,25]+1;到达第三条边，得宽为20-15=5，S=515=75;

第三条边长为25-15=10，同时，发现此边为靠左方块得出边，所以区间[10,20]-1;第三边到第四边的宽为25.5-20=5.5，S=5.5*10=55;
第四条边为出边，后面没有更多的区间，遍历结束。
注： 区间的加1减1作用是计算这个区间覆盖的边的条数
代码实现

#include <iostream>
#include<stdio.h>
#include<algorithm>
#define me(a,b) memset(a,b,sizeof(a))
typedef long long ll;
const int maxn=1e3+10;
using namespace std;
double xx[maxn],tree[maxn<<1];
int cover[maxn<<1];
struct node{
    double lx,rx,y;
    int val;
    node(double a=0,double b=0,double c=0,int d=0):lx(a),rx(b),y(c),val(d){}
    friend bool operator<(node a,node b)
    {
        return a.y<b.y;
    }//重载小于号写在结构体里面，sort可以不用写成sort(a,a+n,cmp);
}line[maxn<<1];
void pushup(int l,int r,int rt)
{
    if(cover[rt]) tree[rt]=xx[r+1]-xx[l];
    else if(l==r) tree[rt]=0;
    else tree[rt]=tree[rt<<1]+tree[rt<<1|1];//将结果往上更新
}
void update(int rt,int l,int r,int lx,int rx,int val)
{
    if(lx<=l&&rx>=r)
    {
        cover[rt]+=val;
        pushup(l,r,rt);
        return;
    }//在查询区间内，进入条件，更新区间覆盖的边的条数，更新区间长度
    int mid=(l+r)>>1;
    if(rx<=mid)
        update(rt<<1,l,mid,lx,rx,val);
    else if(lx>mid)
        update(rt<<1|1,mid+1,r,lx,rx,val);
    else
    {
        update(rt<<1,l,mid,lx,rx,val);
        update(rt<<1|1,mid+1,r,lx,rx,val);
    }
    pushup(l,r,rt);
}
int main()
{
    int t,k=0;
    while(scanf("%d",&t)&&t)
    {
        me(line,0);
        me(xx,0);
        me(cover,0);
        me(tree,0);
        double ans=0;
        int cnt=0;//赋初值
        while(t--)
        {
            double x1,x2,y1,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            line[++cnt]=node(x1,x2,y1,1);//赋入边的值，后面的1为入边区间加1
            xx[cnt]=x1;//获得所有x的值
            line[++cnt]=node(x1,x2,y2,-1);//赋出边区间的值，出边区间减1
            xx[cnt]=x2;//同上
        }
        sort(xx+1,xx+cnt+1);//给xx排序，以便后面离散化
        sort(line+1,line+cnt+1);//按y从小到大排序，因为上图是从小到大遍历的
        int len=unique(xx+1,xx+cnt+1)-(xx+1);//离散化操作，将大数映射成小数，减小内存
        for(int i=1;i<cnt;i++)
        {
            int lx=(int)(lower_bound(xx+1,xx+len+1,line[i].lx)-xx);//找到第i条边对应的离散化后的左横坐标
            int rx=(int)(lower_bound(xx+1,xx+len+1,line[i].rx)-xx-1);//找到第i条边对应的离散化后的右横坐标
            int val=line[i].val;
            update(1,1,len,lx,rx,val);
            ans+=tree[1]*(line[i+1].y-line[i].y);
        }
        printf("Test case #%d\nTotal explored area: %.2f\n",++k,ans);
        printf("\n");
    }
    return 0;
}