圆排列问题

圆排列问题

1. 问题

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

2. 解析

圆排列问题的解空间是一棵排列树。按照回溯法搜索排列树的算法框架，设开始时a=[r1,r2,……rn]是所给的n个元的半径，则相应的排列树由a[1:n]的所有排列构成。
1.center计算圆在当前圆排列中的横坐标，由x^2 = sqrt((r1+r2)^2-(r1-r2)2)推导出x = 2sqrt(r1r2)。

2.Compute计算当前圆排列的长度。变量lenmin记录当前最小圆排列长度。数组r存储所有圆的半径。数组x则记录当前圆排列中各圆的圆心横坐标。

3.在递归算法Backtrack中，当i>n时，算法搜索至叶节点，得到新的圆排列方案。此时算法调用Compute计算当前圆排列的长度，适时更新当前最优值。当i<n时，当前扩展节点位于排列树的i-1层。此时算法选择下一个要排列的圆，并计算相应的下界函数。

·举例：
当n=3，且所给的3个圆的半径分别为1，1，2时，这三个圆有两种可能：


显然第二种情况的长度要小于第一种情况下的长度

3. 设计

double center(int t){
double temp=0;
for(int j=1;j<t;++j) {
		double xvalue=x[j]+2.0*sqrt(r[t]*r[j]); 
		if(xvalue>temp)
			temp=xvalue;
}
	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=high-low;
		for(int i=1;i<N;++i)
			bestr[i]=r[i];
	}
}
void backtrack(int t){
	if(t==N){
		compute();
	}else{
		for(int j=t;j<N;++j){
			swap(r[t],r[j]);
			double centerx=center(t);
			if(centerx+r[t]+r[1]<minlen){
				x[t]=centerx;
				backtrack(t+1);
			}
			swap(r[t],r[j]);
		}
	}
}

4. 分析

时间复杂度：由回溯法的排列树可知，搜索子结点的时间复杂度是O(n!)次，即全排列的时间复杂度，而Backtrack()函数每次计算需要O(n)计算时间，从而整个算法的计算时间复杂性为O((n+1)!)。虽然理论上时间复杂度很大，但实际的消耗实际由于增加了剪枝条件，会比O((n+1)!)小很多。
空间复杂度：由于r[]数组，x[]数组的大小都为n，故空间复杂度为O(n)

5. 源码

https://github.com/Marshmello11/Algorithm/tree/master/FinalExperiment