圆排列问题

1、问题

圆排列问题：给定n个圆的半径序列，将它们放到矩形框中，各圆与矩形底边相切，求具有最小排列长度的圆排列。
例如：当n=3，且所给的3个圆半径分别为1,1,2时，这3个圆的最小长度的圆排列 最小长度为2+4根号2.

2、解析

圆的摆放位置不同，得到的结果可能不同，所以这是一个排列问题，在回溯里面需要添加全排列。

要得到圆排列的长度，需要知道最右边圆心的横坐标,+它的半径+第一个圆的半径【根据算法的约定，第一个圆的圆心横坐标为0，代码中有标记】

根据递归的思想，第n个圆的圆心横坐标=与他相切的圆的横坐标+2*sqrt（r[n]*r[n-1]）[(r1+r2)^2 -(r1-r2)^2 ]

用函数compute()来计算当前圆的排列长度，找到排列的最左端 if(x[i]-r[i]<left) left=x[i]-r[i]; 找到排列的最右端if(x[i]+r[i]>right)//计算出最右端 right=x[i]+r[i]; right-left就得到了圆排列的长度

3、设计

double center(int t){
	double length=0;
	for(int i=1;i<t;i++){
	   double temp=x[i]+2*sqrt((r[t]*r[i]));
		if(temp>length)
		length=temp;
	}
	return length;
} 
void compute(){//计算当前圆排列的长度 
	double left=0;
	double right=0;
	for(int i=1;i<N;i++){
		if(x[i]-r[i]<left)
		left=x[i]-r[i];
		if(x[i]+r[i]>right)
		right=x[i]+r[i];
	} 
	if((right-left)<Min){
	Min=right-left;
	for(int i=1;i<N;i++) 
	best[i]=r[i]; 
	} 
} 
void dfs(int t){
	if(t==N)//最后一个圆 
	compute();
	else{
		for(int i=t;i<N;i++){
			swap(r[t],r[i]);
			double cen=center(t);
		if(r[1]+cen+r[t]<Min){
			x[t]=cen;
			dfs(t+1); 	
		}
		swap(r[t],r[i]);
		}	
	}
} 

4、分析

时间复杂性为O((n+1)!)