POJ 1151 矩形面积并

题解：
把横坐标离散化，然后把矩形拆成上下两条边，上边标成-1，下边标成1，然后按情况增减面积即可~

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <algorithm>

#define N 222

using namespace std;

struct LI
{
    double y,sx,tx;
    int nsx,ntx,fg;
}li[N<<1];

int gs,cs,cnt,n,cas;
int c[N<<1];
double a[N<<1],ans;

inline bool cmp(const LI &a,const LI &b)
{
    return a.y<b.y;
}

inline void read()
{
    gs=cs=0;double tx,ty,sx,sy;
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%lf%lf%lf",&sx,&sy,&tx,&ty);
        a[++gs]=sx; a[++gs]=tx;
        li[++cs].y=sy; li[cs].sx=sx; li[cs].tx=tx; li[cs].fg=-1;
        li[++cs].y=ty; li[cs].sx=sx; li[cs].tx=tx; li[cs].fg=1;
    }
    sort(li+1,li+1+cs,cmp);
    sort(a+1,a+1+gs);
    cnt=unique(a+1,a+1+gs)-a-1;
    for(int i=1;i<=cs;i++)
    {
        li[i].nsx=lower_bound(a+1,a+1+cnt,li[i].sx)-a;
        li[i].ntx=lower_bound(a+1,a+1+cnt,li[i].tx)-a;
    }
}

inline void go()
{
    memset(c,0,sizeof c);
    ans=0.0;
    for(int i=1;i<=cs;i++)
        for(int j=li[i].nsx;j<li[i].ntx;j++)
        {
            if(li[i].fg==-1)
            {
                c[j]--;
                if(c[j]==-1) ans-=(a[j+1]-a[j])*li[i].y;
            }
            else
            {
                c[j]++;
                if(c[j]==0) ans+=(a[j+1]-a[j])*li[i].y;
            }
        }
    printf("Test case #%d\n",++cas);
    printf("Total explored area: %.2lf\n\n",ans);
}

int main()
{
    while(scanf("%d",&n),n) read(),go();
    return 0;
}