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

代码如下：

View Code

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <algorithm>
 4 #include <vector>
 5 #include <iostream>
 6 
 7 using namespace std;
 8 
 9 typedef vector<int> VI;
10 typedef pair<int, int> PII;
11 typedef vector<PII> VPII;
12 typedef pair<double, double> PDD;
13 typedef vector<PDD> VPDD;
14 
15 #define PB push_back
16 #define MPR make_pair
17 #define FI first
18 #define SE second
19 #define REP(i, n) for (int i = 0;  i < (n); i++)
20 const int N = (1 << 6) + 10;
21 const double eps = 1e-10;
22 VPDD rec[N];
23 double ttWt[N], wt[10];
24 
25 int main() {
26 //    freopen("in", "r", stdin);
27 //    freopen("out", "w", stdout);
28     int T, n, x;
29     double res;
30     cin >> T;
31     while (T-- && cin >> res) {
32         cin >> n;
33         REP(i, 1 << n) rec[i].clear();
34         REP(i, n) {
35             cin >> wt[i];
36             rec[1 << i].PB(MPR(0.0, 0.0));
37         }
38         REP(i, 1 << n) {
39             ttWt[i] = 0.0;
40             REP(j, n) {
41                 if (i & (1 << j)) ttWt[i] += wt[j];
42             }
43 //            cout << i << "= =" << ttWt[i] << endl;
44             for (int t = (i - 1) & i; t > 0; t = (t - 1) & i) {
45                 int op = t ^ i;
46 //                cout << i << ' ' << t << ' ' << op << endl;
47                 int sz1 = rec[t].size(), sz2 = rec[op].size();
48                 double cl = ttWt[op] / ttWt[i], cr = ttWt[t] / ttWt[i];
49                 REP(ii, sz1) {
50                     REP(jj, sz2) {
51                         double l1 = rec[t][ii].FI, r1 = rec[t][ii].SE, l2 = rec[op][jj].FI, r2 = rec[op][jj].SE;
52                         rec[i].PB(MPR(max(l1 + cl, l2 - cr), max(r1 - cl, r2 + cr)));
53 //                        REP(j, n) cout << ((t & (1 << j)) != 0);
54 //                        cout << ' ';
55 //                        REP(j, n) cout << ((op & (1 << j)) != 0);
56 //                        cout << ' ';
57 //                        cout << ttWt[t] << ' ' << ttWt[op];
58 //                        cout << "= =" << max(l1 + cl, l2 - 1.0) << ' ' << max(r1 - 1.0, r2 + cr) << endl;
59                     }
60                 }
61             }
62         }
63         int fullSet = (1 << n) - 1, sz = rec[fullSet].size();
64         double mx = -1.0;
65 //        cout << sz << endl;
66         REP(i, sz) {
67             if (rec[fullSet][i].FI + rec[fullSet][i].SE + eps < res) mx = max(mx, rec[fullSet][i].FI + rec[fullSet][i].SE);
68 //            cout << rec[fullSet][i].FI + rec[fullSet][i].SE << endl;
69         }
70         printf("%.16f\n", mx);
71     }
72     return 0;
73 }

——written by Lyon

posted @ 2013-05-04 15:17 LyonLys 阅读(508) 评论(0) 收藏举报

刷新页面返回顶部

poj 2743 && LA 3403 Mobile Computing (mid-easy Search)