poj 2743 && LA 3403 Mobile Computing (mideasy Search)

2743 -- Mobile Computing

https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1404

　　中等偏简单的一道搜索题，用状态压缩可以直接计算出来。因为mobile的一条rod只有两端，所以我们只要考虑两边分别是放置哪个/哪些重物就可以了。这个可以用集合的形式来表示，然后就是一个子集分割，但是这里不允许空集，所以初始集合有所不同。sample真实神奇，我的mobile的力矩公式都用反了居然还可以过了sample。然后换过来了以后，又发现坐标偏移的调整写错了。错了几次，终于A了。- -

代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <iostream>

using namespace std;

typedef vector<int> VI;
typedef pair<int, int> PII;
typedef vector<PII> VPII;
typedef pair<double, double> PDD;
typedef vector<PDD> VPDD;

#define PB push_back
#define MPR make_pair
#define FI first
#define SE second
#define REP(i, n) for (int i = 0;  i < (n); i++)
const int N = (1 << 6) + 10;
const double eps = 1e-10;
VPDD rec[N];
double ttWt[N], wt[10];

int main() {
//    freopen("in", "r", stdin);
//    freopen("out", "w", stdout);
    int T, n, x;
    double res;
    cin >> T;
    while (T-- && cin >> res) {
        cin >> n;
        REP(i, 1 << n) rec[i].clear();
        REP(i, n) {
            cin >> wt[i];
            rec[1 << i].PB(MPR(0.0, 0.0));
        }
        REP(i, 1 << n) {
            ttWt[i] = 0.0;
            REP(j, n) {
                if (i & (1 << j)) ttWt[i] += wt[j];
            }
//            cout << i << "= =" << ttWt[i] << endl;
            for (int t = (i - 1) & i; t > 0; t = (t - 1) & i) {
                int op = t ^ i;
//                cout << i << ' ' << t << ' ' << op << endl;
                int sz1 = rec[t].size(), sz2 = rec[op].size();
                double cl = ttWt[op] / ttWt[i], cr = ttWt[t] / ttWt[i];
                REP(ii, sz1) {
                    REP(jj, sz2) {
                        double l1 = rec[t][ii].FI, r1 = rec[t][ii].SE, l2 = rec[op][jj].FI, r2 = rec[op][jj].SE;
                        rec[i].PB(MPR(max(l1 + cl, l2 - cr), max(r1 - cl, r2 + cr)));
//                        REP(j, n) cout << ((t & (1 << j)) != 0);
//                        cout << ' ';
//                        REP(j, n) cout << ((op & (1 << j)) != 0);
//                        cout << ' ';
//                        cout << ttWt[t] << ' ' << ttWt[op];
//                        cout << "= =" << max(l1 + cl, l2 - 1.0) << ' ' << max(r1 - 1.0, r2 + cr) << endl;
                    }
                }
            }
        }
        int fullSet = (1 << n) - 1, sz = rec[fullSet].size();
        double mx = -1.0;
//        cout << sz << endl;
        REP(i, sz) {
            if (rec[fullSet][i].FI + rec[fullSet][i].SE + eps < res) mx = max(mx, rec[fullSet][i].FI + rec[fullSet][i].SE);
//            cout << rec[fullSet][i].FI + rec[fullSet][i].SE << endl;
        }
        printf("%.16f\n", mx);
    }
    return 0;
}

——written by Lyon