圆排列问题

一、问题

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

二、解析

首先，已知圆的个数n以及记录各圆半径的数组r[i],i=1~n ，记录各圆圆心横坐标的数组x[i],i=1~n。要记录最短排列长度minlen，最终求出排列顺序。要注意，只要大小合适，目标圆就有可能与排列中的任意一个圆相切。要求x[n]时，要从前往后的一一比较，先得到x[1]+a1的值，再得到x[2]+a2的值……一直到x[n-1]+an-1的值，与上一次的值相比较，若距离更大则更新，否则不做处理。
圆排列的解空间是一棵排列树。假设有n个圆排列：1，2，3……n。只要计算出n！/2种排列：

所以，只要知道第一个和第n个位置的数，那么排列时，就只剩下第2个到第n-1个位置。
另外，如果有m个圆的半径相同，那么交换后总长度也相同，所以只要计算一个。
如果有n个圆，前m个已经排好位置，那么还有n-m个待排序定位。当搜索到第m+1个圆时：

两个长度相加就是新的最短长度。

三、设计

1、求每个圆的圆心坐标 ：

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、计算当前圆排列的长度

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];
        if(x[i]+r[i]>high)
            high=x[i]+r[i];
    }
    if(high-low<minlen){//把每次不同的排列长度进行比较，找到更小的minlen
        minlen=high-low;
        for(int i=1;i<N;++i)
            bestsort[i]=r[i];
    }
}

3、回溯

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]);
        }
    }
}

四、分析

时间复杂度：
在回溯算法中，搜索子结点的时间复杂度是O(n!)次，而Backtrack()函数每次计算需要O(n)计算时间，从而整个算法的计算时间复杂性为O((n+1)!)。
空间复杂度：
由于r[]数组，x[]数组的大小都为n，故空间复杂度为O(n)。
源码在这