矩形面积求并（codevs 3044）

/*
  将所有正方形沿线切割成一个个小正方形，然后离散化。 
*/ 
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#define M 210
using namespace std;
double x[M],y[M],s[M][4];
int xy[M][M];
int n,cas=0;
double sum;
int main()
{
    while(cin>>n)
    {   
        if(n==0)break;
        cas++;
        sum=0.0;
        memset(xy,0,sizeof(xy));
        for(int i=1;i<=n;i++)
        {
            cin>>s[i][0]>>s[i][1]>>s[i][2]>>s[i][3];
            x[i*2-2]=s[i][0];
            y[i*2-2]=s[i][1];
            x[i*2-1]=s[i][2];
            y[i*2-1]=s[i][3];
        }
        sort(x,x+2*n);
        sort(y,y+2*n);
        for(int k=1;k<=n;k++)
        {
            int x1,x2,y1,y2;
            for(x1=0;x1<2*n;x1++)
              if(x[x1]==s[k][0])break;
            for(x2=0;x2<2*n;x2++)
              if(x[x2]==s[k][2])break;
            for(y1=0;y1<2*n;y1++)
              if(y[y1]==s[k][1])break;
            for(y2=0;y2<2*n;y2++)
              if(y[y2]==s[k][3])break;
            for(int i=x1;i<x2;i++)
              for(int j=y1;j<y2;j++)
                xy[i][j]=1;
        }
        for(int i=0;i<2*n;i++)
          for(int j=0;j<2*n;j++)
            sum+=xy[i][j]*(x[i+1]-x[i])*(y[j+1]-y[j]);
        printf("%.2f\n",sum);
    }
    return 0;
}