编程题：求多个长方形组成的集合能否覆盖某个长方形

基本问题如下，有一个长方形组成的集合，且长方形的边平行于x轴或y轴，故每个长方形均可以用左上角定点坐标加右下角坐标来确定。给定一个新输入的长方形，问长方形集合能否覆盖这个长方形？例如下图，长方形A、B和C组成了集合，输入的为长方形D,需要判断D是否能被A、B、C覆盖，显然这个例子D是无法被完全覆盖的。

此题其实是一个典型的计算几何问题，解决的问题的方法也比较固定，通常使用平面扫描的方法，我们先考虑一个其他问题：如何计算长方形集合所覆盖的面积，即上图A、B和C三个长方形所覆盖的面积，这个问题其实是类似的。以上图为例，平面扫描的基本思路如下：

从最左边的边开始首先扫到的是A的左边，此时从上往下扫描每个长方形的上下边的y坐标，设长方形A的上下边y坐标分别为和,其他长方形以类似方式表示。此时依次扫描每个y坐标，统计所有在(包含两个坐标)间的y坐标。针对上图示例，我们推断下整个算法过程:

1)扫描到，此时其包含的y坐标为、和，因此在扫描到下一个边即之前，此段覆盖面积为;

2)扫描到,此时其包含的y坐标为、、和，而上一步包含的y坐标再不删除之前仍然起作用（扫描到右边才会删除指定的y坐标，否则一直保留），故此段覆盖面积为;

3)扫描到,此时其包含的y坐标为、、和，这里需要注意这里、、和之前均出现过一次，故这里出现次数为2，要通过两次扫描右边删除才会删除，否则仍然会被用来计算面积，故此段覆盖面积为：;

4)扫描到,此时其包含的y坐标为、和，被删掉，不在用于计算覆盖面积，和出现次数变为1，但仍然用来计算覆盖面积，故此段覆盖面积为：

5）扫描到,此时其包含的y坐标为、和，被删掉，不在用于计算覆盖面积，和出现次数变为1，但仍然用来计算覆盖面积，故此段覆盖面积为：。

至此整个过程完成，所有覆盖面积求和，可以得到总的覆盖面积，下面给出C++代码。

#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
struct Rectangle
{
	double x1, y1, x2, y2;
};
double unionArea(const vector<Rectangle>& rects) {
	if (rects.empty())
		return 0.0;
	//事件信息：（x坐标，左侧还是右侧(1和-1)，长方形的序号）。
	typedef pair<double, pair<int, double> > Event;
	vector<Event> events;
	vector<double> ys;
	//遍历各长方形并找出y坐标和事件的集合。
	for (int i = 0; i < rects.size(); ++i)
	{
		//依次将上边和下边坐标压入vector。
		ys.push_back(rects[i].y1);
		ys.push_back(rects[i].y2);
		//依次将事件压入vector。
		events.push_back(Event(rects[i].x1,make_pair(1,i)));
		events.push_back(Event(rects[i].x2, make_pair(-1, i)));
	}
	//对y坐标集合进行排序并去重。
	sort(ys.begin(), ys.end());
	ys.erase(unique(ys.begin(),ys.end()));
	//对事件目录排序，按照x坐标从小到大排序。
	sort(events.begin(),events.end());
	int ret = 0;
	//count[i]=ys[i]~ys[i+1]区间中重叠的四边形个数。
	vector<int> count(ys.size() - 1, 0);
	for (int i = 0; i < events.size(); ++i)//从最左边的边开始循环。
	{
		//x表示x轴的坐标，delta：1表示为长方形的左边，-1表示长方形的右边。
		int x = events[i].first, delta = events[i].second.first;
		//rectangle表示长方形为原始输入中的第几个。
		int rectangle = events[i].second.second;
		//更新 count[]
		//y1:当前边（可能左边也可能右边）对应的长方形的上边坐标。
		//y2:当前边（可能左边也可能右边）对应的长方形的下边坐标。
		int y1 = rects[rectangle].y1, y2 = rects[rectangle].y2;
		for (int j = 0; j < ys.size(); j++)
		{
			if (y1 <= ys[j] && ys[j] < y2)
				count[j] += delta;
		}
		//计算cutLength
		double cutLength = 0.0;
		for (int j = 0; j < ys.size() - 1; j++)
		{
			if (count[j] > 0)
				cutLength += ys[j + 1] - ys[j];
		}
		if (i + 1 < events.size())
			ret += cutLength * (events[i + 1].first - x);
	}
	return ret;
}
int main()
{

	return 0;
}

在解决了这一问题之后，再求我们题目的问题就比较简单了，同样以上图为例我们咋回到长方形D开始边的x坐标。检测此位置能否被全覆盖，然后同样扫描各个边，直到长方形D的右边，一旦在某个边的位置无法全覆盖那么就返回False,如果扫描到最右边依然没问题，那么返回True。