POJ1151 离散化求矩形面积的并

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

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

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

 int main()
 {
  int k=;
  while(~scanf("%d",&n)&&n)
  {
   int cnt=;
   for(int i=;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=;
   memset(flag,false,sizeof(flag));
   solve();
   printf("Test case #%d\n",k++);
   printf("Total explored area: %.2lf\n\n",sum);
  }
 }