Poj1151线段树求面积并

简单模拟一下，如图三个矩形，分为将其上下底边由下至上标号为1,2...6，每扫到下底边，该底边所覆盖的范围cover+1，每扫描到下底边，其覆盖的范围cover-1

过程如下图

这样依次计算，当扫描到第n-1条边的时候就计算出来面积的并

其中还是要用到离散化的知识

代码

#include <stdio.h>
#include <algorithm>

using namespace std;

const int maxn=205;

struct node
{
    int l,r,cover;
    double ly,ry,len;
}segTree[maxn*3];

struct line
{
    double x,y1,y2;
    int tag;//表示是底边(1)还是顶边(-1)
    bool operator < (const line &l)const {
        return x<l.x;
    }

}lines[maxn];

double y[maxn];

void build(int i,int l,int r)
{
    segTree[i].l=l;
    segTree[i].r=r;
    segTree[i].len=0;
    segTree[i].cover=0;
    segTree[i].ly=y[l];
    segTree[i].ry=y[r];
    if(l+1==r)return;
    int mid=(l+r)/2;
    build(i<<1,l,mid);
    build(i<<1|1,mid,r);
}

void cal(int i)
{
    if(segTree[i].cover>0){
        segTree[i].len=segTree[i].ry-segTree[i].ly;
        return;
    }
    if(segTree[i].l+1==segTree[i].r)segTree[i].len=0;
    else segTree[i].len=segTree[i<<1].len+segTree[i<<1|1].len;
}

void update(int i,line &l)
{
    if(l.y1==segTree[i].ly&&l.y2==segTree[i].ry){
        segTree[i].cover+=l.tag;
        cal(i);
        return;
    }
    if(l.y2<=segTree[i<<1].ry)update(i<<1,l);
    else if(l.y1>=segTree[i<<1|1].ly)update(i<<1|1,l);
    else{
        int y=l.y2;
        l.y2=segTree[i<<1].ry;
        update(i<<1,l);
        l.y2=y;
        y=l.y1;
        l.y1=segTree[i<<1|1].ly;
        update(i<<1|1,l);
        l.y1=y;
        /*line tmp(l);
        tmp.y2=segTree[i<<1].ry;
        update(i<<1,tmp);
        tmp=l;
        tmp.y1=segTree[i<<1|1].ly;
        update(i<<1|1,tmp);*/
    }
    cal(i);//相当于一个push_up
}

int main()
{
    //freopen("in.txt","r",stdin);
    int T=1,n,tot;
    double x1,y1,x2,y2;
    while(scanf("%d",&n)!=EOF&&n){
        tot=1;
        for(int i=0;i<n;i++){
            scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
            lines[tot].x=x1;
            lines[tot].y1=y1;
            lines[tot].y2=y2;
            lines[tot].tag=1;
            y[tot]=y1;
            ++tot;
            lines[tot].x=x2;
            lines[tot].y1=y1;
            lines[tot].y2=y2;
            lines[tot].tag=-1;
            y[tot]=y2;
            ++tot;
        }
        sort(lines+1,lines+tot);
        sort(y+1,y+tot);
        build(1,1,tot-1);
        update(1,lines[1]);
        double area=0;
        for(int i=2;i<tot;i++){
            area+=segTree[1].len*(lines[i].x-lines[i-1].x);
            update(1,lines[i]);
        }
        printf("Test case #%d\nTotal explored area: %.2f\n\n",T++,area);
    }
    return 0;
}