pku 1151 （简单离散化求矩形覆盖总面积）

    这是道简单计算几何题。只用了离散化，没有线段树。离散化是这样一种思想：n个矩形分布在平面中他们的坐标可能分布很广。考虑这样两矩形(1,3,2,4),(1000000,1000000,2000000,2000000)要是按模拟的方法来做在时间和空间上会遭成浪费，我们只是想知道矩形的长和宽，要是用两个1维数组记录下矩形的坐标的，我们只要记录下4条边的覆盖情况。于是，将连续的空间离散成来x[0]=1,x[1]=2,x[2]=10000000,x[3]=2000000.y[0]=3,y[1]=4,y[2]=10000000,y[3]=200000000.在用个辅助2维数组
记录横纵线的覆盖情况(x，y数组的点看成横纵线)。

1. /*pku1151
2.   Name: Atlantis
3.   Date: 28-07-08 18:54
4.   Description: 离散化
5. */
6. #include<stdio.h>
7. #include<math.h>
8. #include<stdlib.h>
9. #include<algorithm>
10. #define pr printf
11. double x[202],y[202];//最多100个矩形，所以最多只有200条横线或纵线
12. double a[102][4];//矩形实际坐标
13. bool cover[202][202];
14. const double EPS = 1e-7;
15. int cmp(const void * b,const  void * a){//b>a 返回1
16.          if(fabs(*(double* )b-*(double *)a)<EPS)return 0;
17.             else if(*(double* )b-*(double *)a>0)return 1;
18.              else return -1;
19. }
20. int main(){
21.     int n,i,k,j,h1,h2,v1,v2,index=1;
22.     while(scanf("%d",&n)){
23.         if(n==0)break;
24.     for(i=0,k=0;i<n;i++,k++){
25.         scanf("%lf%lf%lf%lf",&a[i][0],&a[i][1],&a[i][2],&a[i][3]);
26.         x[k]=a[i][0],y[k]=a[i][1],x[++k]=a[i][2],y[k]=a[i][3];
27.         }
28.         qsort(x,2*n,sizeof(x[0]),cmp);
29.         qsort(y,2*n,sizeof(y[0]),cmp);//将横线，竖线排好序
30.         memset(cover,0,sizeof(cover));//初始时无覆盖区
31.         //对每个矩形，进行一次横纵线的扫描，覆盖范围
32.         for(i=0;i<n;i++)
33.            {
34.                 k=0;
35.                 while(fabs(x[k]-a[i][0])>EPS)k++;h1=k;
36.                 k=0;
37.                 while(fabs(y[k]-a[i][1])>EPS)k++;v1=k;
38.                 k=0;
39.                 while(fabs(x[k]-a[i][2])>EPS)k++;h2=k;
40.                 k=0;
41.                 while(fabs(y[k]-a[i][3])>EPS)k++;v2=k;
42.                 //标记该矩形覆盖的区域
43.                 for(j=h1;j<h2;j++)      //不考虑最右
44.                     for(k=v1;k<v2;k++)  //不考虑最下
45.                         cover[j][k]=true;
46.            }      
47.            double sum=0;           
48.            //计算总面积
49.          for(i=0;i<2*n-1;i++)
50.             for(j=0;j<2*n-1;j++)
51.                     sum+=(cover[i][j]*(x[i+1]-x[i])*(y[j+1]-y[j]));
52.          pr("Test case #%d/n",index++);
53.          pr("Total explored area: %.2lf/n/n",sum);
54.     }
55. }