圆排列问题

Ø给定n个大小不等的圆c1,c2,…,cn，现要将这n个圆排进一个矩形框中，且要求各圆与矩形框的底边相切。圆排列问题要求从n个圆的所有排列中找出有最小长度的圆排列。

Ø例如，当n=3，且所给的3个圆的半径分别为1，1，2时，这3个圆的最小长度的圆排列如图所示。其最小长度为2+4√2。

void Circle::Backtrack(int t)

{

if (t>n) Compute();

else

for (int j = t; j <= n; j++) {

Swap(r[t], r[j]);

float centerx=Center(t);

if (centerx+r[t]+r[1]<min)

{//下界约束

x[t]=centerx;

Backtrack(t+1);}

Swap(r[t], r[j]);}

}

float Circle::Center(int t)

{// 计算当前所选择圆的圆心横坐标

float temp=0;

for (int j=1;j<t;j++) {

float valuex=x[j]+2.0*sqrt(r[t]*r[j]);

if (valuex>temp) temp=valuex;

}

return temp;

}

void Circle::Compute(void)

{// 计算当前圆排列的长度

float 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<min) min=high-low;

}

Ø上述算法尚有许多改进的余地。例如，象1,2,…,n-1,n和n,n-1, …,2,1这种互为镜像的排列具有相同的圆排列长度，只计算一个就够了，可减少约一半的计算量。

Ø另一方面，如果所给的n个圆中有k个圆有相同的半径，则这k个圆产生的k! 个完全相同的圆排列，只计算一个就够了。