编程之美 1.7 光影切割问题

光影切割问题

原书的题目表述一开始又没看懂，还是直接数学化吧，一张图胜就明白了。

问题：求区域A-B之间，被直线划分成多少区域，前提是不存在三条直线相交的情况。
解题思路：
1、问题化解
仔细观察发现：划分区域 = 直线数 + 交点数+1
问题化解为：区域内直线的交点数。
2、解法一：求区域内的直线交点
3、解法二：求序列的倒序。

仔细观察，左边交点为a>b>c，右边交点为c>b>a。则右边的序列逆序是3。刚好就是直线的交点数目。也就是存在一个逆序直线就有一个交点，这也很好想象。

逆序解释：如序列（2，4，3，1）逆序为某数之前比它大的数目：0+0+1+3=4。
示例代码：

#include <vector>
#include <iostream>
using namespace std;

/**********************************/
/*这里的直线都是截距式表示*/
/**********************************/

// 解法一：遍历直线，寻找所有直线的交点
typedef struct
{
	float x;
	float y;
}Point;

typedef struct
{
	float k;
	float b;
}Line;

typedef vector<Line>Lines;
typedef vector<Point>Points;

// 寻找两直线的交点
Point findIntersection(Line l1, Line l2)
{
	Point point;

	point.x = (l2.k - l1.k) / (l1.b - l2.b);
	point.y = l1.k * point.x + l1.b;

	return point;
}

// 寻找直线序列的交点
Points getIntersections(Lines lines)
{
	Points points;
	Point point;
	for(int i = 0; i < lines.size() - 1; ++i)
	{
		for (int j = i + 1; j < lines.size(); ++j)
		{
			if (lines.at(i).k != lines.at(j).k)
			{
				point = findIntersection(lines.at(i), lines.at(j));
				points.push_back(point);
			}
		}
	}

	return points;
}

// 寻找区域内直线的交点
int getRegionIntersections(Points points, float A, float B)
{
	int nums = 0;
	int size = points.size();
	for (int i = 0; i < size; ++i)
	{
		if (points.at(i).x >= A &&points.at(i).x <= B)
		{
			nums++;
		}
	}

	return nums;
}

int MethodOne(Lines lines, float A, float B)
{
	// 交点数
	Points points = getIntersections(lines);

	// 区域内交点数
	int nums = getRegionIntersections(points, A, B);
	return lines.size() + nums + 1;
}

/*****************************************************************/
// 解法二：
// 对第一个点进行从大到小排序
Point getMaxPoint(Points points, int* pos)
{
	float max = -INFINITY ;
	Point maxPoint = points.front();
	int size = points.size();
	for (int i = 0; i < size; ++i)
	{
		if (max < points.at(i).x)
		{
			max = points.at(i).x;
			maxPoint = points.at(i);
			*pos = i;
		}
	}
	return maxPoint;
}

// 对点序列第一元素x按大到小排序
Points sort(Points points)
{
	Points result;
	int pos;
	result.clear();
	Point maxPoint;
	while (!points.empty())
	{
		maxPoint = getMaxPoint(points, &pos);
		result.push_back(maxPoint);
		points.erase(pos+ points.begin());
	}

	return result;
}

// 求逆序数量
int getInvertedSorts(Points points)
{
	int nums = 0;
	int size = points.size();
	for (int i = 1; i < size; ++i)
	{
		for (int j = 0; j < i; ++j)
		{
			if (points.at(i).y > points.at(j).y)
			{
				nums++;
			}
		}
	}

	return nums;
}

int MethodTwo(Lines lines, float A, float B)
{
	// 这里的点是直线与A、B竖线的交点。
	Points points;
	Point point;

	// 求与竖线 A B 的交点
	int size = lines.size();
	for (int i = 0; i < size; ++i)
	{
		point.x = lines.at(i).k * A + lines.at(i).b;
		point.y = lines.at(i).k * B + lines.at(i).b;
		points.push_back(point);
	}

	 // 对竖线A交点的排序
	points = sort(points);

	// 对竖线B交点相对于竖线A交点的逆序
	int nums = getInvertedSorts(points);

	return size + nums+1;

	
}
const float A = 0;
const float B = 100;
int main()
{
	Lines lines;
	Line line;
	line.k = 1;
	line.b = 1;
	lines.push_back(line);

	line.k = 0;
	line.b = 2;
	lines.push_back(line);

	line.k = 0;
	line.b = 3;
	lines.push_back(line);

	line.k = 4;
	line.b = 6;
	lines.push_back(line);

	line.k = -1;
	line.b = 6;
	lines.push_back(line);

	cout << "方法一划分区域数量为： " << MethodOne(lines, A, B) << endl;
	cout << "方法二划分区域数量为： " << MethodTwo(lines, A, B) << endl;
	getchar();
	return 0;
}