poj1151 矩形面积并

解题报告

题目
：http://poj.org/problem?id=1151

题目大意
：给定平面上n个矩阵（n<100）求总面积

算法
：离散化 +

染色

思路
：如果点都是整数，而且矩阵左下角和右上角的点的坐标很小，那么我们可以开一个二维数组area[MAX_X][MAX_Y]，然后如果有矩阵（1，1，3，4）我们就将x 1到3，y 1到4标记为1，当有矩阵

（2，3，4，5）那我们在把2到4，3到5标记为1.这样就解决的矩阵间重复的问题。最后看一下area里有几个1.
现在矩阵规模没有改变，但是点的坐标是浮点型而且最大100000，我们把它转化为上面的模型，将点离散化为1，2，3.。。。。。。然后染色。然后看为1的小格area[i][j],那么这个格对应于原来的面积就是(ori[i + 1] – ori[i]) * (ori[j + 1] *
ori[j])
ori[i]表示i对应着的原来浮点型的点。

离散化是注意判重。离散化用map实现就行。

提交情况
：runtime error
1次



Wrong answer 1次



Presentation error 1次



Accepted 1次

经验与收获
：练习了map的同时，学习了矩形面积并的解决方法。

AC code
:

#include <cstdio>

#include <map>

#include <algorithm>



using namespace std;



#define MAXN 410



struct POINT{


double x, y;

};



struct _Matrix{


POINT left, right;

} matrix[MAXN];



double point[MAXN * 2], ori[MAXN * 2];



map<double,int> Pt;

map<double,int>::iterator it;

bool area[MAXN][MAXN];



bool comp(const double &A, const double &B){


return A < B;

}



int main(){


int n, i, j, k, m, key, CASE = 0;


while(~scanf("%d", &n), n){


m = 0;


for(i = 0; i < n; i ++){


scanf("%lf %lf %lf %lf", &matrix[i].left.x, &matrix[i].left.y, &matrix[i].right.x, &matrix[i].right.y);


point[m ++] = matrix[i].left.x;


point[m ++] = matrix[i].left.y;


point[m ++] = matrix[i].right.x;


point[m ++] = matrix[i].right.y;


}


sort(point, point + m, comp);


key = 1;


Pt[point[0]] = key;


ori[key] = point[0];


for(i = 1; i < m; i ++){


if(point[i] != point[i - 1]){


Pt[point[i]] = ++key;


ori[key] = point[i];


}


}


memset(area, 0, sizeof(area));


for(i = 0; i < n; i ++){


for(j = Pt[matrix[i].left.x]; j < Pt[matrix[i].right.x]; j ++)


for(k = Pt[matrix[i].left.y]; k < Pt[matrix[i].right.y]; k ++)


area[j][k] = 1;


}


double ans = 0;


for(i = 1; i <= key; i ++)


for(j = 1; j <= key; j ++)


if(area[i][j]){


ans += (ori[i + 1] - ori[i]) * (ori[j + 1] - ori[j]);


}


printf("Test case #%d\n", ++CASE);


printf("Total explored area: %.2lf\n\n", ans);


}


return 0;

}