- 问题 给定n个大小不等的圆c1,c2,…,cn,现要将这n个圆排进一个矩形框中,且要求各圆与矩形框的底边相切。圆排列问题要求从n个圆的所有排列中找出有最小长度的圆排列。
- 解析主要部分分为三个函数:
1.center求每个圆的圆心坐标,利用直角三角形的定理,可以得出横坐标x的计算公式x^2 = sqrt((r1+r2)2-(r1-r2)2)推导出x = 2sqrt(r1r2)。如下图,三个圆半径分别为1,1,2cm,可得总长为4根号2+2,只要大小合适,目标圆就有可能与排列中的任意一个圆相切
所以只要大小合适,目标圆就有可能与排列中的任意一个圆相切如上图二所示,第二个圆和第三个圆并没有相交,所以在求x[3]的坐标的时候,x[1]+a的值和x[2]+b的值要进行比较,哪一个大则把哪一个作为最小的长度,之后利用for循环不断重复此动作
2.Compute计算当前圆排列的长度。变量minlen顾名思义就是表示最小圆排列的长度。数组r存储所有圆的半径。数组x则记录当前圆排列中各圆的圆心横坐标。如图所示,当其中一个圆的半径远远大于其他的圆时,因为已知所有圆的x[]和r[],很容易求出每个圆的左右坐标,通过比较找出最小的左部坐标和最大的右部坐标,一减就是该圆排列的长度,然后把每次不同的排列长度相比较,找到更小的minlen就更新。
3.在递归算法Backtrack中,当i>n时,算法搜索至叶节点,得到新的圆排列方案。此时算法调用Compute计算当前圆排列的长度,适时更新当前最优值。当i<n时,当前扩展节点位于排列树的i-1层。此时算法选择下一个要排列的圆,并计算相应的下界函数。
如图所示,当半径为1,1,2时,共有6种排列办法,ABC,ACB,BAC,BCA,CAB,CBA
3.设计
- center函数:求每个圆的圆心坐标
double center(int m){//计算所有圆的圆心横坐标
double a=0;
for(int j=1;j<m;++j){//判断目标圆是否与排在它之前的任一圆相切
double b=x[j]+2.0*sqrt(r[m]*r[j]);//通过直角三角形的简单定理,就可以计算圆的横坐标
if(b>a)
a=b;
}
return a;
}
2.compute函数:计算当前圆排列的长度
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];//不断通过for循环更新左边的low与右边的high
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)
bestsort[i]=r[i];
}
}
3、backtrack函数:回溯算法
void backtrack(int m){
if(m==N)
compute();
else{
for(int j=m;j<N;++j){
swap(r[m],r[j]);//交换两个位置
double center_x=center(m);
if(center_x+r[m]+r[1]<minlen){//先判断当前层是否在范围内,如果是那么继续不断搜索下一层,如果不是直接回溯。
x[m]=center_x;
backtrack(m+1);
}
swap(r[m],r[j]);
}
}
}
4.分析时间复杂度:
在回溯算法中,搜索子结点的时间复杂度是O(n!)次,而Backtrack()函数每次计算需要O(n)计算时间,从而整个算法的计算时间复杂性为O((n+1)!)
空间复杂度:
由于r[]数组,x[]数组的大小都为n,故空间复杂度为O(n)
5.源码
https://github.com/hackkkkkk/calculate