poj1151--计算几何+离散化

//把相交的矩形分成若干个小矩形
//vis[i][j]表示由标号为xi和yj的小矩形是不是原来矩形的一部分
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<queue>
#define dist(a,b) sqrt(1.0*(a.x-b.x)*(a.x-b.x)+1.0*(a.y-b.y)*(a.y-b.y))
#define cross(a,b,c) 1.0*(b.x-a.x)*(c.y-a.y)-1.0*(b.y-a.y)*(c.x-a.x)
#define dot(a,b,c) (b.x-a.x)*(c.x-a.x)+(b.y-a.y)*(c.y-a.y)
#define delt(a) fabs(a)<eps?0:a>0?1:-1
#define pi acos(-1.0)
#define eps 1e-8
#define inf 1e20
#define N 1005
#define M 1005
using namespace std;
typedef __int64 i64;
struct TPoint
{
	double x1,y1,x2,y2;
}pt[N],st;
double lx[N],ly[N];
bool vis[N][N];
int n,m,t;
bool scan()
{
	scanf("%d",&n);
	if(n==0) return 0;
	m=0;
	for(int i=0;i<n;i++)
	{
		scanf("%lf%lf%lf%lf",&pt[i].x1,&pt[i].y1,&pt[i].x2,&pt[i].y2);
		lx[m]=pt[i].x1,ly[m++]=pt[i].y1;
		lx[m]=pt[i].x2,ly[m++]=pt[i].y2;
	}
}
void solve(int cases)
{
	int i,k,j,i1,j1,i2,j2;
	double sumf=0;
	sort(lx,lx+m);
	sort(ly,ly+m);
	memset(vis,0,sizeof(vis));
	for(i=0;i<n;i++)
	{
		i1=lower_bound(lx,lx+m,pt[i].x1)-lx;
		i2=lower_bound(lx,lx+m,pt[i].x2)-lx;
		j1=lower_bound(ly,ly+m,pt[i].y1)-ly;
		j2=lower_bound(ly,ly+m,pt[i].y2)-ly;
		for(j=i1;j<i2;j++)
		{
			for(k=j1;k<j2;k++) vis[j][k]=1;
		}
	}
	for(i=0;i<m;i++)
	{
		for(j=0;j<m;j++)
		{
			sumf=sumf+vis[i][j]*(lx[i+1]-lx[i])*(ly[j+1]-ly[j]);
		}
	}
	printf("Test case #%d\nTotal explored area: %.2f\n\n",cases,sumf); 
}
int main()
{
	int cases=1;
	while(scan())
	{
		solve(cases++);
	}
    return 0;
}