yuan

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

2. 解析
   使用分支界限解决，遍历所有排序情况，通过剪枝减去排列比当下长的序列的情况及其子结点，从求出具有最小排列长度的圆序列。
   假定给n个圆的半径，那么排列情况共有n!个情况，计算着n种情况的排列长度，比较得出最小值以及圆排列。可以利用回溯法搜索树的框架。
   要求出最小圆半径排列，相邻两个圆必定相切或者相离。

3. 设计
   
   

#include "stdio.h"
#include "math.h"
#define MAX 10

int n;  //圆的个数
float min = 10000;  //当前圆排列的最小长度
int r[MAX];  //各个圆的半径
int bestr[MAX];   //最优排列下圆的半径
float x[MAX];  //各个圆的圆心

void swap(int &a, int &b)
{
	int temp = a;
	a = b;
	b = temp;
}

//计算当前圆排列的长度
void compute()
{
	int i;
	float low = 0;  //圆位于原点右侧
	float heigh = 0;
	for (i = 1; i <= n; i++)
	{
		if (x[i] - r[i] < low)  //寻找最左边的坐标，当圆位于原点左侧时才会更新
			low = x[i] - r[i];
		if (x[i] + r[i] > heigh)  //寻找最右边的坐标
			heigh = x[i] + r[i];
	}
	if (heigh - low < min)  //如果当前圆排列的长度小于之前的最优值
		min = heigh - low;  //更新
}

//计算当前所选择圆的圆心
float center(int i)
{
	int j;
	float temp = 0.0;
	float valuex;
	for (j = 1; j < i; j++)
	{
		valuex = x[j] + 2.0*sqrt(r[j] * r[i]);  //寻找当前层下最右边的圆心
		if (valuex > temp)
			temp = valuex;
	}
	return temp;
}

void backtrack(int i)
{
	if (i > n)  //如果到最后一个圆
	{
		compute();  //计算当前圆排列的最小长度
		int  j;
		for (j = 1; j <= n; j++)
			bestr[j] = r[j];
	}
	else
	{
		int j;
		for (j = i; j <= n; j++)  //选择要排列的圆
		{
			swap(r[i], r[j]); //选择第j个圆
			float centerx = center(i);  //计算圆心
			if (centerx + r[i] < min)  //下界约束，当前圆的长度<最优值，才能产生最优值
			{
				x[i] = centerx;  //记录圆心
				backtrack(i + 1);  //搜索下一层
			}
			swap(r[i], r[j]);  //返回上一层
		}
	}
}

void circle(int n1, int r1[])
{
	n = n1;
	int i;
	for (i = 1; i <= n; i++)
		r[i] = r1[i];
	backtrack(1);
}

int main()
{

	int rr[] = { 0, 1, 1, 2 };
	int nn = 3;

	circle(nn, rr);
	printf("圆的个数为：%d\n", n);
	printf("半径分别为：\n");
	int i;
	for (i = 1; i <= n; i++) {
		printf("%d ", rr[i]);
	}
	printf("\n");
	printf("最优排列下，圆的半径为：\n");
	for (i = 1; i <= n; i++) {
		printf("%d ", bestr[i]);
	}
	printf("\n");
	printf("最小长度为：%.3f\n", min);
	return 0;
}

算法时间复杂度为O(n！)