POJ1151 离散化求矩形面积的并

/*第一道离散化的题目，虽然是水题，不过还是很高兴。。。*/

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
struct rect
{
 double x1,x2,y1,y2;
};
#define max 103<<1
rect a[max>>1];
double x[max],y[max];
bool flag[max][max];
int n;
double sum;

void solve()
{
 for(int i=0;i<n;i++)
 {
  for(int j=0;j<n<<1;j++)
  {
   if(x[j]>=a[i].x2)
    break;
   if(x[j]<a[i].x1)
    continue;
   for(int k=0;k<n<<1;k++)
   {
    if(y[k]>=a[i].y2)
     break;
    if(y[k]<a[i].y1)
     continue;
    flag[j][k]=true;
   }
  }
 }
 for(int i=0;i<n<<1;i++)
  for(int j=0;j<n<<1;j++)
  {
   if(flag[i][j])
    sum+=(x[i+1]-x[i])*(y[j+1]-y[j]);
  }
}

int main()
{
 int k=1;
 while(~scanf("%d",&n)&&n)
 {
  int cnt=0;
  for(int i=0;i<n;i++)
  {
   scanf("%lf%lf%lf%lf",&a[i].x1,&a[i].y1,&a[i].x2,&a[i].y2);
   x[cnt]=a[i].x1;y[cnt]=a[i].y1;cnt++;
   x[cnt]=a[i].x2;y[cnt]=a[i].y2;cnt++;
  }
  sort(x,x+cnt);
  sort(y,y+cnt);
  sum=0;
  memset(flag,false,sizeof(flag));
  solve();
  printf("Test case #%d\n",k++);
  printf("Total explored area: %.2lf\n\n",sum);
 }
}