//0 1背包
#include<iostream>
#include<stdio.h>
#include<string.h>
using namespace std;
int main()
{
int i,j,n,v,f[1100],w[1100],p[1100];
scanf("%d",&t);
while(t--)
{
memset(f,0,sizeof(f));
scanf("%d%d",&n,&v);
for(i=0;i<n;i++)
scanf("%d",&p[i]);
for(i=0;i<n;i++)
scanf("%d",&w[i]);
for(i=0;i<n;i++)
{
for(j=v;j>=w[i];j--)
f[j]=f[j]>(f[j-w[i]]+p[i])?f[j]:(f[j-w[i]]+p[i]);
}
printf("%d\n",f[v]);
}
return 0;
}
//完全背包
#include<stdio.h>
#define max 999999
int main()
{
int i,j,n,v,f[10010],w[10010],p[10010],t,a,b;
scanf("%d",&t);
while(t--)
{
scanf("%d%d%d",&a,&b,&n);
v=b-a;
for(i=1;i<=v;i++)
f[i]=max;
f[0]=0;
for(i=0;i<n;i++)
scanf("%d%d",&p[i],&w[i]);
for(i=0;i<n;i++)
for(j=w[i];j<=v;j++)
if(f[j]>f[j-w[i]]+p[i])
f[j]=f[j-w[i]]+p[i];
if(f[v]==max) printf("This is impossible.\n");
else
printf("The minimum amount of money in the piggy-bank is %d.\n",f[v]);
}
return 0;
}
01背包是指每个物品只能用一次
完全背包指的是每个物品都能无限使用
没什么技巧，DP类的明确了状态后就写起来快了，写多了程序自然也就有经验了，有经验就能快速的推出状态，总之就是多写题
程序的不同很简单
完全背包
for i:=1 to n do //枚举1-N的物品
for j:=a[i] to m do //枚举1-M的能背出来的范围
f[j]:=f[j] or f[j-a[i]];
01背包
for i:=1 to n do //枚举1-N的物品
for j:=m downto a[i] do //枚举1-M的能背出来的范围
f[j]:=f[j] or f[j-a[i]];
注意，之前要先f[0]:=true；
可以看出，区别只在第二个循环的正倒，F是一个布尔数组，F[i]表示i这个数字能够被组合出来
为什么循环的正倒会导致这样的区别呢？
我们举个例子：
假设我们现在组合出了3，现在讨论的物品体积是1，那么,即a[i]=1, f[3]:=true;
当正循环时
当枚举到4的时候，f[4]=f[4] or f[4-1]=true，因为f[3]=true;
当枚举到5的时候，f[5]=f[5] or f[5-1]=true，因为f[4]=true;
显然，这个“1”会被无限的使用下去直到到达M的上限
逆循环则相反
复制的
标注：下文n即为物品件数，c[i]表示第i件物品的耗费（体积），V为背包容量，a[i]表示第i件物品的价值 dp[]数组存放的即为最优解。
一、01背包
最简单的背包，每件物品选或者不选。
1 for(int i = 1; i <= n; ++i){
2 for(int j = V; j >= c[i]; --j){
3 dp[j] = max(dp[j],dp[j-c[i]]+a[i]);
4 }
5 }
二、完全背包
每件物品可以选无限次
1 for(int i = 1; i <= n; ++i){
2 for(int j = c[i]; j <= V; ++j){
3 dp[j] = max(dp[j],dp[j-c[i]] + a[i];
4 }
5 }
注意：初始化方面的细节，如果要求恰好装满背包，那么在初始化时除了dp[0]为0其它dp[1..V]均设为-∞（求的是最大解，如果求的是最小解，则为∞） 如果并没有要求必须把背包装满，而是只希望价格尽量大，初始化时应该将f[0..V]全部设为0。
三、多重背包
每件物品可以选有限次 设每件物品的个数为count[i];
O(V*Σlog count[i])写法（二进制优化）
1 for(int i = 1; i <= n; ++i){
2 int k;
3 for(k = 1; k*2 < count[i] + 1; k *= 2){
4 for(int j = V; j >= k*c[i] ; --j){
5 dp[j] = max(dp[j],dp[j-k*c[i]] + k*a[i]);
6 }
7 }
8 k = count[i] + 1 - k; //k等于前面项数之和+1,（k-1）就是前面项数之和，1-k为其相反数
9 for(int j = V; j >= k*c[i]; --j){
10 dp[j] = max(dp[j],dp[j-k*c[i]] + k*a[i]);
11 }
12
13 }
O(VN)的写法
/*此种方法求最大值不确定可否，但判断是否能将1-V内的背包装满可以*/
----------------------------------------------------------------------
1 int used[MAXN];
2 for(int i = 1; i <= N; ++i){
3 memset(used,0,sizeof(used));
4 for(int j = c[i]; j <= V; ++j){
5 if(used[j-c[i]] < count[i] && dp[j] < dp[j-c[i]] + a[i]){
6 dp[j] = dp[j-c[i]] + a[i];
7 used[j] = used[j-c[i]] + 1;
8 }
9 }
10 } //有待验证，估计不对，有兴趣的可以验证下
----------------------------------------------------------------------
现在dp[i]用来表示容量为i的背包能否被所给物品恰好装满 上面的01背包和完全背包以及上面一种O(V*Σlog count[i])的写法也可以做这样的改变
1 memset(dp,false,sizeof(dp));
2 dp[0] = true;
3 for(int i = 1; i <= N; ++i){
4 memset(used,0,sizeof(used));
5 for(int j = c[i]; j <= V; ++j){ //注意这里是顺序的（和完全背包相同），其实多重背包也可以理解成被限制了的完全背包
6 if(!dp[j] && dp[j-c[i]] && used[j-c[i]] < count[i]){ //注意这里的!dp[j]不能少
7 dp[j] = true;
8 used[j] = used[j-c[i]] + 1;
9 }
10 }
11 }
四、混合背包
以上三种背包的混合，有的物品只能取一次，有的物品能取无限次，有的物品能取有限次
算法伪代码：
for i=1..n
if 第i件物品是01背包
for j=V...c[i]
dp[j]=max(); //同上面01背包代码
else if 第i件物品是完全背包
for j=c[i]...V
dp[j]=max(); //同上面完全背包
else if 第i件物品是多重背包
MultiplePack(c[i],a[i],count[i]) //同上多重背包
五、二维费用的背包
如问题：每件物品都有各自的体积v[i]和重量w[i]，一个背包最多可以装的体积为V，最多可以装的重量为W,问最多能装的价值。即有两种属性限制的背包问题。再比如，一个一般的01背包问题，再加上个数限制，即最多拿K件物品，也可以理解为每件物品的个数耗费为1，背包再加一个容量为K的属性，也是二维费用背包。对于二维费用的背包，只需对dp数组加一维即可，01背包倒序，完全背包正序
以01背包为例：
for i=1...n
for j=V..v[i]; //若为完全背包则写为for j=v[i]..V for k=w[i]..W;
for k=W..w[i];
dp[j][k] = max(dp[j][k],dp[j-v[i]][k-w[i]] + a[i]);
六、分组背包的问题
将n件物品分为若干组，每组中只能选择一件物品，问可以得到的最大价值。
for 所有的组k
for v=V..0
for 所有的i属于组k
dp[v]=max{dp[v],dp[v-c[i]]+a[i]}
七、有依赖的背包问题
要购买一件物品，必须先购买另一件物品。这就是有依赖关系的背包问题 背包九讲里讲的很清楚了，我就不多说了。
下面两个部分暂且放下。以后遇到具体的问题时再补充。
八、泛化物品的背包
九、背包问题问法的变化
下面就说几个比较经典的题目：01背包 完全背包就不说了，混合背包也没必要说
一、【POJ1014】【TOJ1034】【HDU1059】【ZOJ1149】Dividing——多重背包问题
一般在几个OJ上同时出现并且长久不衰的题都是经典题目，比如这题，是一道关于多重背包的题目
【题目大意】6种价值分别为1 2 3 4 5 6的石头，给出每种石头的数目，问能不能将这些石头分为价值相等的两堆。
【解析】经典的多重背包问题，a[]数组即为1 2 3 4 5 6，count数组题目给出，怎么做呢？上面提到一点，设dp[i]表示容量为i的背包能不能被所给物体装满。如果dp[sum/2]==1，那么就说明可以分为两堆（sum是所有石头的总价值），直接套用上面的模板即可。
1 #include <cstdio>
2 #include <cstring>
3 using namespace std;
4
5 int a[10];
6 int dp[120100];
7 int main(){
8 int cas = 0;
9 while( ++ cas ){
10 int sum = 0;
11 for(int i = 1; i <= 6; ++i){
12 scanf("%d",&a[i]);
13 sum += a[i]*i;
14 }
15 if(sum == 0)break;
16
17 memset(dp,0,sizeof(dp));
18 printf("Collection #%d:\n",cas);
19
20 if(sum % 2 == 1){
21 puts("Can't be divided.");
22 puts("");
23 continue;
24 }
25
26 dp[0] = 1;
27 //O（NV）的方法
28 int dpt[120000];
29 for(int i = 1; i <= 6; ++i){
30 memset(dpt,0,sizeof(dpt));
31 for(int j = i; j <= sum/2; ++j){
32 if (!dp[j] && dp[j-i]&& dpt[j-i] < a[i]) {
33 dpt[j] = dpt[j-i] + 1;
34 dp[j] = 1;
35 }
36 }
37 }
38 /*
39 //二进制优化的方法
40
41 for(int i = 1;i <= 6;++i){
42 int k;
43 for(k = 1;k*2 < a[i]+1;k *= 2)
44 for(int v = sum/2;v >= k*i;--v)
45 if(dp[v-k*i])
46 dp[v]=true;
47 k = a[i]+1-k;
48 for(int v = sum/2;v > k*i;--v)
49 if([v-k*i])
50 dp[v] = true;
51
52 }
53 */
54 if(dp[sum/2]){
55 puts("Can be divided.");
56 }
57 else {puts("Can't be divided.");};
58 puts("");
59
60 }
61 return 0;
62 }
二、【TOJ3540】Consumer有依赖关系的背包
【题目大意】给若干组物品，每组物品都有一个箱子（箱子自身也有cost），然后就是物品的cost和value，要买某个物品必须也要买装这个物品的箱子，给一定钱数，问能获得的最大价值。
【题目解析】可以参照背包九讲第六讲和第七讲里的内容解决该问题（相当对口，这题就是根据背包九讲来的）
有一种比较简单的写法，两次背包即可。
1 #include<iostream>
2 using namespace std;
3 const int MAXN = 100000;
4 int main()
5 {
6 int n,total,boxcost,goodnum,cost,value,i,j,t;
7 int dpbox[MAXN],dptotal[MAXN];
8 while(scanf("%d%d",&n,&total) != EOF)
9 {
10 memset(dptotal,0,sizeof(dptotal));
11 for(i = 0; i < n; i++)
12 {
13 scanf("%d%d",&boxcost,&goodnum);
14 memcpy(dpbox,dptotal,sizeof(dptotal));
15 for(j = 0; j < goodnum; j++)
16 {
17 scanf("%d%d",&cost,&value);
18 for(t = total - boxcost; t >= cost; t--)
19 {
20 if(dpbox[t] < dpbox[t - cost] + value)
21 dpbox[t] = dpbox[t-cost] + value;
22 }
23 }
24 for(t = total; t >= boxcost; t--)
25 if(dptotal[t] < dpbox[t-boxcost])
26 dptotal[t] = dpbox[t-boxcost];
27 }
28 printf("%d\n",dptotal[total]);
29 }
30 return 0;
31 }
这题源自多校联合赛，是一个比较难的背包综合题
【题目大意】一个人有T分钟去做工作，共有n组工作，每一组有若干工作，每个工作消耗的时间和得到的快乐值已知，每一个组有一个标识位，标识位为0，则该组中至少应该选取一样工作完成，标识位为1，则该组中至多有一个工作被选择，标识位为2，无限制。
【分析】对标识位为2的工作组，因为无限制条件，每个工作可做可不做，即用简单的01背包即可处理。对标识位为1的工作组，因为至多选择一个工作，所以本次操作只对dp数组选取最优值更新一次即可。对至少选择一个的工作组，还不是太理解，只是看别人代码写的，还请大家集思广益，有想法可以留言评论，以后理解透彻了以后会做补充。
1 #include <cstdio>
2 #include <algorithm>
3 #include <cstring>
4 #include <cmath>
5 #include <queue>
6 #include <vector>
7 #include <map>
8
9 using namespace std;
10
11 const int maxt = 100 + 10;
12
13 int dp[maxt];
14 int td[maxt], n, m, t;
15
16 void dp0() { //至少取一个
17 memcpy (td, dp, sizeof(td));
18 memset (dp, -1, sizeof(dp));
19 for (int i = 0, a, b; i < m; ++i) {
20 scanf ("%d%d", &a, &b);
21 for (int j = t; j >= a; --j) {//注意怎样保证至少取一个？？
22 if (dp[j-a] != -1) {
23 dp[j] = max(dp[j], dp[j-a] + b);
24 }
25 if (td[j-a] != -1) {
26 dp[j] = max(dp[j], td[j-a] + b);
27 }
28 }
29 }
30 }
31
32 void dp1() { //最多取一个,用上次得到的dp值更新最优值，只更新一次
33 memcpy(td, dp, sizeof(td));
34 for (int i = 0, a, b; i < m; ++i) {
35 scanf ("%d%d", &a, &b);
36 for (int j = t; j >= a; --j) {
37 if (td[j-a] != -1) { //用-1初始化，这种情况能够被装入的都是恰好被装入
38 dp[j] = max(dp[j], td[j-a] + b);
39 }
40 }
41 }
42 }
43
44 void dp2() { //无条件限制，01背包，更新多次
45 for (int i = 0, a, b; i < m; ++i) {
46 scanf ("%d%d", &a, &b);
47 for (int j = t; j >= a; --j) {
48 if (dp[j-a] != -1) {
49 dp[j] = max(dp[j], dp[j-a] + b);
50 }
51 }
52 }
53 }
54
55 void init() {
56 memset (dp, -1, sizeof(dp));
57 dp[0] = 0;
58 for (int i = 0, s; i < n; ++i) {
59 scanf ("%d%d", &m, &s);
60 if (s == 0) dp0();
61 else if (s == 1) dp1();
62 else if (s == 2) dp2();
63 }
64 }
65
66 void solve() {
67 int ans = -1;
68 for (int i = 0; i <= t; ++i)
69 ans = max(ans, dp[i]);
70 printf ("%d\n", ans);
71 }
72
73 int main() {
74 //freopen("ar.txt","r",stdin);
75 while (scanf ("%d%d", &n, &t) != EOF) {
76 init();
77 solve();
78 }
79 return 0;
80 }
融云，免费为你的App加入IM功能——让你的App“聊”起来！！
· 初创企业：融资有风险 额度需谨慎
· 安卓用户不满，Adobe将推出更多安卓类应用
· 微软发布Win10预览版build 10041
· Moto 360寻谷歌替代品 加速入华
· 这家面部表情追踪公司获欧盟投资
» 更多新闻...
· 在团队中如何带领新手
· 微博背后的那些算法
· 谈“测试驱动的开发”
· 好团队不可能凭空出现，赢在Leader的可行规划