算法分析与实践-大作业-圆排序问题

1.问题

圆排列问题：给定n个圆的半径序列，将它们放到矩形框中，各圆与矩形底边相切，求具有最小排列长度的圆排列。

2.解析

圆排列问题主要用到了回溯法。问题的解空间是一棵排列数。按照回溯法遍历每种排列，求得最小排列长度的最优解。
变量n为圆的个数
数组r[n]存储每个圆的半径，x[n]存储每个圆的圆心坐标，bestr[n]存储最小圆排序的半径顺序。
minlen存储最小排列长度。
center函数计算该圆在排列中的横坐标（首个圆的圆心坐标为0）。
compute函数计算当前排列长度。
Backtrack函数为核心函数，通过递归遍历每种排列。

第一次遍历得到ABC的排列，然后通过不断递归回溯，得到最优解。

3.设计

double center(int t)//得到每个圆的圆心坐标
{
	double temp = 0;
	for (int j = 1; j < t; ++j)//存在跨圆相切的情况，所以循环判断
	{
		double value = x[j] + 2.0*sqrt(r[t] * r[j]);
		if (value > temp)
			temp = value;
	}
	return temp;
}
void compute()//计算最小排列长度
{
	double low = 0, high = 0;
	for (int i = 1; i <= n; ++i)//求出圆排序后的矩形左边框和右边框坐标
	{
		if (x[i] - r[i] < low)
			low = x[i] - r[i];
		if (x[i] + r[i] > high)
			high = x[i] + r[i];
	}
	if (high - low < minlen)//更新minlen和bestr数组
	{
		minlen = high - low;
		for (int i = 1; i <= n; ++i)
			bestr[i] = r[i];
	}
}
void backtrack(int t)
{
	if (t == n+1)
	{
		compute();
	}
	else
	{
		for (int j = t; j <= n; ++j)
		{
			double temp;//全排列
			temp = r[t];
			r[t] = r[j];
			r[j] = temp;
			double centerx = center(t);
			if (centerx + r[t] + r[1] < minlen)//因为第一个圆的圆心坐标为0，所以要加上r[1]
			{
				x[t] = centerx;
				backtrack(t + 1);
			}
			temp = r[t];
			r[t] = r[j];
			r[j] = temp;
		}
	}
}

4.分析

backtrack函数的复杂度为O(n!)，每次排列都需要调用一次compute函数，而compute函数复杂度为O(n)，所以整个算法复杂度为O(n*n!)

5.源码

https://github.com/LonelyTaker/Algorithm-analysis

参考：https://blog.csdn.net/qq_37373250/article/details/81477394