线段树--扫描线

引入

扫描线一般是用在求图形面积上，由于图形会重合，而数据范围又太大，这时就要用扫描线了。
由于扫描线需要修改区间，所以要使用线段树–lazy标记还有一些小技巧。

例题

题目描述
桌面上放了N个平行于坐标轴的矩形，这N个矩形可能有互相覆盖的部分，求它们组成的图形的面积。
输入
有多组测试数据。

每组测试数据输入第一行为一个数N，表示矩形的数量。

下面N行，每行四个整数，分别表示每个矩形的左下角和右上角的坐标，坐标范围为–10^8到108之间的整数。

当N==0时，测试文件结束。

输出
对于每个测试用例，在单独一行中输出一个整数，表示图形的面积。
样例输入
1
1 2 3 4
3
1 1 4 3
2 -1 3 2
4 0 5 2
0
样例输出
4
10

分析

求重叠图形面积，也就是引入中提到的扫描线。下面具体说说怎么实现。
如果有如下丑图：

其实矩形可以化简成两条竖线，我们可以定义一个结构体Line{x,y,yy,flag},表示一条边，x是横坐标，y和yy是边的两个端点，flag==1表示是矩形的左边，flag ==-1表示是矩形的右边。
如下丑图：

那么每条线可以映射在一条主线上，我们可以先把端点去重排序，然后以端点划分出区间。
如下丑图：

那么问题就变成了维护一段区间，可以用线段树实现。
在主线上加上一条边：找到这条边的两个端点在哪个区间，然后update即可。
因为端点排过一次序，所以找端点所在的区间可以用二分查找。
最后是求面积，面积=长*宽=两条扫描线的距离 *主线上每个区间的长总和。

代码

#include<bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define full(a,b) memset(a,b,sizeof a)
#define N 2000+5
int n;
long long ans;
int t,tmp[2*N],dy[2*N];//两倍大小好存两条边 
struct Line//存边 
{
    int x,y,yy,flag;
} line[2*N];
struct Tree//线段树 
{
    int cnt;
    long long len;
} tree[10*N];
void init()//初始化 
{
    ans=0;t=0;
}
int cmp(Line a,Line b)
{
    return a.x<b.x;
}
void dist()//去重排序 
{
    sort(tmp+1,tmp+1+n);
    dy[++t]=tmp[1];
    for(int i=2; i<=n; i++)
        if(tmp[i]!=tmp[i-1]) dy[++t]=tmp[i];
    --t;//t个点，t-1个区间 
    return;
}
void built(int k,int l,int r)//建树 
{
    tree[k].len=0;//长为0 
    tree[k].cnt=0;//覆盖次数=0 
    if(l==r) return;
    int mid=(l+r)>>1;
    built(2*k,l,mid);
    built(2*k+1,mid+1,r);
}
int half(int x)//二分查找 
{
    int l=1,r=t+1;//找点要t+1 
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(dy[mid]==x) return mid;
        if(dy[mid]<x) l=mid+1;
        else r=mid-1;
    }
}
void pushdown(int k,int l,int r)
{
    if(tree[k].cnt) tree[k].len=dy[r+1]-dy[l];//覆盖过直接求 
    else tree[k].len=tree[2*k].len+tree[2*k+1].len;//找儿子 
}
void update(int k,int l,int r,int x,int y,int v)
{
    if(x>r||y<l) return;
    if(x<=l&&r<=y)
    {
        tree[k].cnt+=v;
        pushdown(k,l,r);//懒惰 
        return;
    }
    int mid=(l+r)>>1;
    update(2*k,l,mid,x,y,v);
    update(2*k+1,mid+1,r,x,y,v);
    pushdown(k,l,r);//回溯 
}
int main()
{
//  freopen("「VOJ1056」图形面积 加强版.in","r",stdin);
//  freopen("「VOJ1056」图形面积 加强版.out","w",stdout);
    while(~scanf("%d",&n)&&n)
    {
        init();
        for(int i=1; i<=n; i++)
        {
            int a,b,c,d;
            scanf("%d%d%d%d",&a,&b,&c,&d);
            line[2*i-1]=Line{a,b,d,1};//存边 
            line[2*i]=Line{c,b,d,-1};
            tmp[2*i-1]=b;//记录端点 
            tmp[2*i]=d;
        }
        n*=2;//边是矩形的两倍 
        dist();
        built(1,1,t);
        sort(line+1,line+1+n,cmp);
        for(int i=1; i<=n; i++)
        {
            if(i>1) ans+=(line[i].x-line[i-1].x)*tree[1].len;//算面积 
            int x=half(line[i].y),y=half(line[i].yy)-1;//找端点所在区间 
            update(1,1,t,x,y,line[i].flag);//映射到主线上 
        }
        printf("%lld\n",ans);
    }
    return 0;
}