UVa 1354 天平难题（Mobile Computing）

题目比较长，需要题目的参照其他文章。
看了很多解法，觉得都不是很容易看懂。
尤其是刘汝佳自身的版本。

下面写一种易懂的版本

首先直接看这幅图。
实际上就是枚举二叉树，
每次任选两个节点，然后当成一个新的结点，并且这个节点的质量就是任选的两个节点质量之和。
这样一直反复枚举下去，就可以得到最后节点的头。
这里的二叉树是自底向上形成的。

然后每个头节点的最左值显然等于左边的值 + 左边节点已经有的左值 或是右边节点有的左值，减去右边的值。
这里也比较难理解一点。
一开始其实每个节点的左右值都为0，因此节点的左右值是一个相对属性。
都是相对于根节点的左右值，因此某个根节点的左值只会更左， 而另一个节点的左值，需要减掉这个根节点的右值，因为这个根节点的右节点的左是相对于这个根节点的右节点的，要让他相对于这个根节点，那么应该减去这个根节点的右值。这里也比较难理解一点。

刘汝佳版本快的不行， 我执行5w组程序可能要1分钟
刘汝佳版本只要1s不到
如果想要深入算法的推荐看刘汝佳的版本。

下面是易懂版本。

#include<bits/stdc++.h>
#define LL long long
using namespace std;

double maxWid;
double nowWid;
int num;

struct Node{
    double left, right;
    double w;
}n[1<<6];

bool vis[1<<6];

void dfs(int cur) {
    if(cur == num - 1) {
        double width = n[2 * cur].left + n[2 * cur].right;
        if(width < maxWid && width > nowWid) {
            nowWid = width;
        } 
        return;
    }
    for(int i = 0; i < cur + num; i++) {
        for(int j = 0; j < cur + num; j++) {
            if(i == j) continue;
            if(!vis[i] && !vis[j]) {
                vis[i] = vis[j] = 1;
                n[cur + num].w = n[i].w + n[j].w;
                double leftV = n[j].w / n[cur + num].w;
                double rightV = n[i].w / n[cur + num].w;
                n[cur + num].left = max(n[i].left + leftV,
                n[j].left - rightV);
                n[cur + num].right = max(n[i].right - leftV,
                n[j].right + rightV);
                dfs(cur + 1);
                vis[i] = vis[j] = 0;
            }
        }
    }
}

int main() {
    //freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
    ios::sync_with_stdio(false);
    cin.tie(0);
    int T;
    cin >> T;
    while(T--) {
        memset(vis, 0, sizeof(vis));
        memset(n, 0, sizeof(n));
        cin >> maxWid >> num;
        nowWid = 0;
        for(int i = 0; i < num; i++) 
        cin >> n[i].w;
        dfs(0);
        if(nowWid == 0 && num != 1) {
            cout << "-1" << endl;
        } else {
            cout << fixed << setprecision(15) << nowWid << endl;
        }
    }
    return 0;
}