超大01背包问题(二分)

  1. #include   
  2. #include   
  3. #include   
  4. using namespace std;  
  5. const int N = 50;  
  6. const long long INF = 0x3fffffff;  
  7. typedef long long LL;  
  8. int n;  
  9. LL w[N], v[N];  
  10. LL W;  
  11. pair  pi[1 << (N / 2)];  
  12.   
  13. void solve() {  
  14.     int n2 = n / 2;  
  15.     for(int i = 0; i < (1 << n2); i++) {  
  16.         LL sw = 0, sv = 0;  
  17.         for(int j = 0; j < n2; j++) {  
  18.             if((i >> j) & 1) {  
  19.                 sw += w[j];  
  20.                 sv += v[j];  
  21.             }  
  22.         }  
  23.         pi[i] = make_pair(sw, sv);  
  24.     }  
  25.   
  26.     sort(pi, pi + (1 << n2));  
  27.     int m = 1;  
  28.     for(int i = 1; i < (1 << n2); i++) {  
  29.         if(pi[m-1].second < pi[i].second) {  
  30.             pi[m++] = pi[i];  
  31.         }  
  32.     }  
  33.   
  34.     LL res = 0;  
  35.     for(int i = 0; i < (1 << (n - n2)); i++) {  
  36.         LL sw = 0, sv = 0;  
  37.         for(int j = 0; j < n - n2; j++) {  
  38.             if((i >> j) & 1) {  
  39.                 sw += w[n2 + j];  
  40.                 sv += v[n2 + j];  
  41.             }  
  42.         }  
  43.         if(sw <= W) {  
  44.             LL tv = (lower_bound(pi, pi + m, make_pair(W - sw, INF)) - 1)->second;  
  45.             res = max(res, sv + tv);  
  46.         }  
  47.     }  
  48.     printf("%lld\n", res);  
  49. }  
  50.   
  51. int main() {  
  52.     while(~scanf("%d%lld", &n, &W)) {  
  53.         for(int i = 0; i < n; i++) {  
  54.             scanf("%lld%lld", &w[i], &v[i]);  
  55.         }  
  56.         solve();  
  57.     }  
  58.     return 0;  
  59. }  

你可能感兴趣的:(大三ACM,-,DP,挑战程序设计)