大作业：圆排列问题

1. 问题 给定n个大小不等的圆c1,c2,…,cn，现要将这n个圆排进一个矩形框中，且要求各圆与矩形框的底边相切。圆排列问题要求从n个圆的所有排列中找出有最小长度的圆排列。
   
2. 解析主要部分分为三个函数：
   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.设计

1. 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