矩形面积求并（codevs 3044）

题目描述 Description

输入n个矩形，求他们总共占地面积（也就是求一下面积的并）

输入描述 Input Description

可能有多组数据，读到n=0为止(不超过15组)

每组数据第一行一个数n，表示矩形个数(n<=100)

接下来n行每行4个实数x1,y1,x2,y1(0 <= x1 < x2 <= 100000;0 <= y1 < y2 <= 100000)，表示矩形的左下角坐标和右上角坐标

输出描述 Output Description

每组数据输出一行表示答案

样例输入 Sample Input

2
10 10 20 20
15 15 25 25.5
0

样例输出 Sample Output

180.00

数据范围及提示 Data Size & Hint

无

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