【POJ1151】Atlantis-线段树+离散化+扫描线

测试地址：Atlantis
题目大意：平面上有n个边与坐标轴平行的矩形，求这些矩形的面积并。
做法：本题需要用到线段树+离散化+扫描线。
矩形面积并是线段树的一个经典应用，然而我一直拖到现在才了解了解法，由于不会插图，我在这里讲不清楚，大家可以看这篇文章。
我在这里只讲一下处理的技巧：这里维护动态线段并时涉及到区间操作，然而并不需要使用Lazy-Tag，因为线段插入和删除是一一对应的，所以直接把信息存放在询问到的那些节点，然后pushup即可。具体的处理请参看代码。
以下是本人代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int n,y[110][2];
double x1[110],y1[110],x2[110],y2[110];
struct forsort
{
    double val;
    int id,idd;
}f[1010];
struct oper
{
    double x;
    int y1,y2,add;
}p[1010];
struct segnode
{
    int cov;
    double sum;
}seg[4010];

bool cmp(forsort a,forsort b)
{
    return a.val<b.val;
}

bool cmp0(oper a,oper b)
{
    return a.x<b.x;
}

void pushup(int no,int l,int r)
{
    if (seg[no].cov) seg[no].sum=f[r+1].val-f[l].val;
    else if (l==r) seg[no].sum=0.0;
         else seg[no].sum=seg[no<<1].sum+seg[no<<1|1].sum;
}

void buildtree(int no,int l,int r)
{
    seg[no].cov=0,seg[no].sum=0.0;
    if (l==r) return;
    int mid=(l+r)>>1;
    buildtree(no<<1,l,mid);
    buildtree(no<<1|1,mid+1,r); 
}

void modify(int no,int l,int r,int s,int t,int d)
{
    if (l>=s&&r<=t)
    {
        seg[no].cov+=d;
        pushup(no,l,r);
        return;
    }
    int mid=(l+r)>>1;
    if (s<=mid) modify(no<<1,l,mid,s,t,d);
    if (t>mid) modify(no<<1|1,mid+1,r,s,t,d);
    pushup(no,l,r);
}

int main()
{
    int t=0;
    while(scanf("%d",&n)&&n)
    {
        t++;
        for(int i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1[i],&y1[i],&x2[i],&y2[i]);
            f[2*i-1].val=y1[i],f[2*i-1].id=i,f[2*i-1].idd=0;
            f[2*i].val=y2[i],f[2*i].id=i,f[2*i].idd=1;
        }
        sort(f+1,f+2*n+1,cmp);
        for(int i=1;i<=2*n;i++)
            y[f[i].id][f[i].idd]=i;
        for(int i=1;i<=n;i++)
        {
            p[2*i-1].x=x1[i],p[2*i].x=x2[i];
            p[2*i-1].y1=p[2*i].y1=y[i][0];
            p[2*i-1].y2=p[2*i].y2=y[i][1];
            p[2*i-1].add=1,p[2*i].add=-1;
        }
        sort(p+1,p+2*n+1,cmp0);

        buildtree(1,1,2*n);
        double ans=0.0;
        for(int i=1;i<2*n;i++)
        {
            modify(1,1,2*n,p[i].y1,p[i].y2-1,p[i].add);
            ans+=(p[i+1].x-p[i].x)*seg[1].sum;
        }
        if (t>1) printf("
");
        printf("Test case #%d
Total explored area: %.2f
",t,ans);
    }

    return 0;
}