poj 1151 线段树求矩形面积的并

其实掌握了基本思路就没问题了，不过还是能学到一些思想的。比如说，离散化，扫描线，计算几何里的线段树建树的不同。

尤其是要注意离散化，许多线段树问题都要离散化才能做，否则超内存。这题的离散化就是把坐标的范围，用存放它们的下标来表示。

这样一来范围就小多了，然后每个节点表示的范围就是这两个下标 存放的对应值之间的范围。

扫描线的方法见这里：http://www.cnblogs.com/Booble/archive/2010/10/10/1847163.html

这个博客讲的虽然是周长，但是扫面线的思想是一样的。

本题代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAX 111
int n;
double ans=0.0;
struct Node
{
    int l,r;
    int c;
    double lf,rf,sum;
};
struct Line
{
    double x,y1,y2;
    int f;
};
Node node[MAX*6];
Line line[MAX*3];
double Y[MAX*3];
int cmp(Line a,Line b)
{
    return a.x<b.x;
}
void build_tree(int l,int r,int k)
{
    node[k].l=l;node[k].r=r;
    node[k].c=0;
    node[k].sum=0;
    node[k].lf=Y[l];
    node[k].rf=Y[r];
    if(l+1==r) return ;
    int mid=(l+r)/2;
    build_tree(l,mid,k*2);
    build_tree(mid,r,k*2+1);
}
void callen(int k)
{
    if(node[k].c>0)
        node[k].sum=node[k].rf-node[k].lf;
    else if(node[k].l+1==node[k].r) node[k].sum=0;
    else node[k].sum=node[k*2].sum+node[k*2+1].sum;
}
void update(Line L,int k)
{
    if(L.y1>=node[k].rf||L.y2<=node[k].lf)
        return ;
    if(L.y1<=node[k].lf&&node[k].rf<=L.y2)
    {
        node[k].c+=L.f;
        callen(k);
        return ;
    }
    update(L,k*2);
    update(L,k*2+1);
    callen(k);
}
int main()
{
    int i,j;
    int T=0;
    while(scanf("%d",&n),n)
    {
        T++;
        int t=1;
        double x1,y1,x2,y2;
        for(i=1;i<=n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            line[t].x=x1;
            line[t].y1=y1;
            line[t].y2=y2;
            line[t].f=1;
            Y[t++]=y1;
            line[t].x=x2;
            line[t].y1=y1;
            line[t].y2=y2;
            line[t].f=-1;
            Y[t++]=y2;
        }
        sort(line+1,line+t,cmp);
        sort(Y+1,Y+t);
        build_tree(1,t-1,1);
        ans=0;
        for(i=1;i<t-1;i++)
        {
            update(line[i],1);
            ans+=node[1].sum*(line[i+1].x-line[i].x);
        }
        printf("Test case #%d\n",T);
        printf("Total explored area: %0.2lf\n\n",ans);
    }
    return 0;
}