编程算法 - 背包问题(三种动态规划) 代码(C)

背包问题(三种动态规划) 代码(C)


本文地址: http://blog.csdn.net/caroline_wendy


题目参考: http://blog.csdn.net/caroline_wendy/article/details/37912949


可以用动态规划(Dynamic Programming, DP)求解, 可以通过记忆化搜索推导出递推式, 可以使用三种不同的方向进行求解.

动态规划主要是状态转移, 需要理解清晰.


代码:

[cpp]  view plain copy
  1. /* 
  2.  * main.cpp 
  3.  * 
  4.  *  Created on: 2014.7.17 
  5.  *      Author: spike 
  6.  */  
  7.   
  8. /*eclipse cdt, gcc 4.8.1*/  
  9.   
  10. #include <stdio.h>  
  11. #include <memory.h>  
  12. #include <limits.h>  
  13.   
  14. #include <utility>  
  15. #include <queue>  
  16. #include <algorithm>  
  17.   
  18. using namespace std;  
  19.   
  20. class Program {  
  21.     static const int MAX_N = 100;  
  22.   
  23.     int n=4, W=5;  
  24.     int w[MAX_N] = {2,1,3,2}, v[MAX_N]={3,2,4,2};  
  25.     int dp[MAX_N+1][MAX_N+1]; //默认初始化为0  
  26. public:  
  27.     void solve() {  
  28.         for (int i=n-1; i>=0; i--) {  
  29.             for (int j=0; j<=W; ++j) {  
  30.                 if (j<w[i]) {  
  31.                     dp[i][j] = dp[i+1][j];  
  32.                 } else {  
  33.                     dp[i][j] = max(dp[i+1][j], dp[i+1][j-w[i]] + v[i]);  
  34.                 }  
  35.             }  
  36.         }  
  37.         printf("result = %d\n", dp[0][W]);  
  38.     }  
  39.   
  40.     void solve1() {  
  41.         for (int i=0; i<n; ++i) {  
  42.             for (int j=0; j<=W; ++j) {  
  43.                 if (j<w[i]) {  
  44.                     dp[i+1][j] = dp[i][j];  
  45.                 } else {  
  46.                     dp[i+1][j] = max(dp[i][j], dp[i][j-w[i]]+v[i]);  
  47.                 }  
  48.             }  
  49.         }  
  50.         printf("result = %d\n", dp[n][W]);  
  51.     }  
  52.   
  53.     void solve2() {  
  54.         for (int i=0; i<n; i++) {  
  55.             for (int j=0; j<=W; ++j) {  
  56.                 dp[i+1][j] = max(dp[i+1][j], dp[i][j]);  
  57.                 if (j+w[i]<=W) {  
  58.                     dp[i+1][j+w[i]] = max(dp[i+1][j+w[i]], dp[i][j]+v[i]);  
  59.                 }  
  60.             }  
  61.         }  
  62.         printf("result = %d\n", dp[n][W]);  
  63.     }  
  64. };  
  65.   
  66.   
  67. int main(void)  
  68. {  
  69.     Program P;  
  70.     P.solve2();  
  71.     return 0;  
  72. }  

输出:

[plain]  view plain copy
  1. result = 7  


节省空间, 可以使用1维数组的动态规划.

代码:

[cpp]  view plain copy
  1. /* 
  2.  * main.cpp 
  3.  * 
  4.  *  Created on: 2014.7.17 
  5.  *      Author: spike 
  6.  */  
  7.   
  8. /*eclipse cdt, gcc 4.8.1*/  
  9.   
  10. #include <stdio.h>  
  11. #include <memory.h>  
  12. #include <limits.h>  
  13.   
  14. #include <utility>  
  15. #include <queue>  
  16. #include <algorithm>  
  17.   
  18. using namespace std;  
  19.   
  20. class Program {  
  21.     static const int MAX_N = 100;  
  22.   
  23.     int n=4, W=5;  
  24.     int w[MAX_N] = {2,1,3,2}, v[MAX_N]={3,2,4,2};  
  25.     int dp[MAX_N+1];  
  26. public:  
  27.     void solve() {  
  28.         memset(dp, 0, sizeof(dp));  
  29.         for (int i=0; i<n; ++i) {  
  30.             for (int j=W; j>=w[i]; --j) {  
  31.                 dp[j] = max(dp[j], dp[j-w[i]]+v[i]);  
  32.             }  
  33.         }  
  34.         printf("result = %d\n", dp[W]);  
  35.     }  
  36. };  
  37.   
  38.   
  39. int main(void)  
  40. {  
  41.     Program P;  
  42.     P.solve();  
  43.     return 0;  
  44. }  

输出:

[plain]  view plain copy
  1. result = 7  

你可能感兴趣的:(算法,动态规划)