超大背包问题 折半枚举法

1. #include <cstdio>  
2. #include <algorithm>  
3. #include <map>  
4. #include <vector>  
5. using namespace std;  
6. typedef long long ll;  
7. const int maxn = 40 + 5;  
8. const int INF = 10000000;  
9.   
10. int n;  
11. ll w[maxn], v[maxn];  
12. ll W;  
13.   
14. pair<ll, ll> ps[1 << (maxn / 2)];    //（重量， 价值）  
15. void solve()  
16. {  
17.     //枚举前半部分  
18.     int n2 = n / 2;  
19.     for (int i = 0; i < 1 << n2; i++){  
20.         ll sw = 0, sv = 0;  
21.         for (int j = 0; j < n2; j++){  
22.             if (i >> j & 1){  
23.                 sw += w[i];  
24.                 sv += v[i];  
25.             }  
26.         }  
27.         ps[i] = make_pair(sw, sv);  
28.     }  
29.   
30.     //去除多余的元素  
31.     sort(ps, ps + (1 << n2));  
32.     int m = 1;  
33.     for (int i = 1; i < 1 << n2; i++){  
34.         if (ps[m - 1].second < ps[i].second){  
35.             ps[m++] = ps[i];  
36.         }  
37.     }  
38.   
39.     //枚举后半部分并求解  
40.     ll res = 0;  
41.     for (int i = 0; i < 1 << (n - n2); i++){  
42.         ll sw = 0, sv = 0;  
43.         for (int j = 0; j < n - n2; j++){  
44.             if (i >> j & 1){  
45.                 sw += w[n2 + j];  
46.                 sv += v[n2 + j];  
47.             }  
48.         }  
49.         if (sw <= W){  
50.             ll tv = (lower_bound(ps, ps + m, make_pair(W - sw, INF)) - 1) -> second;  
51.             res = max(res, sv + tv);  
52.         }  
53.     }  
54.     printf("%lld\n", res);  
55. }