算法分析与设计（圆排列问题）

#include<bits/stdc++.h>
using namespace std;
int N;//定义圆的数量;
double r[10],x[10];//存储圆的半径和圆心x坐标
double minlen=10000;//最小圆排列长度
double best[10];//存取最小圆排列半径的方案
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 Calculate() {//计算整个圆排列的长度
	double min=0,max=0;
	for(int i=1; i<N; ++i) {
		if(x[i]-r[i]<min)//如果圆坐标减圆的半斤小于最小值
			min=x[i]-r[i];
		if(x[i]+r[i]>max)//如果圆坐标加圆的半径大于最大值
			max=x[i]+r[i];
	}
	if(max-min<minlen) {
		minlen=max-min;
		for(int i=1; i<N; ++i)
			best[i]=r[i];
	}
}
void PTree(int t) {//用排列树的结构排列所有圆排列的可能性
	if(t==N) {
		Calculate();//如果排列树的层数等于圆的数量则可以计算整个圆排列的长度
	} 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;
				PTree(t+1);
			}
			swap(r[t],r[j]);
		}
	}
}
int main() {

	printf("请输入你需要计算的圆的数量\n");
	scanf("%d",&N);
	printf("请输入%d个圆的半径\n",N);
	N=N+1;
	
	for(int i=1; i<N; i++) {
		scanf("%lf",&r[i]);
	}
	printf("\n");
	PTree(1);
	printf("最小圆排列长度%f\n",minlen);
	printf("最优圆排列半径顺序：\n");
	for(int i=1;i<N;i++){
		printf("%f ",best[i]);
	} 
	
}