01背包问题
2014-11-13 18:33
197 查看
近期在看动态规划,了解到理解动态规划思想比较容易的一个题就是01背包问题,动手写了下程序。
重要的是要理解动态规划的思想
找到的一些关于01背包问题的求解思想:
首先介绍一下动态规划...
设计一个动态规划算法,通常可按照以下几个步骤进行:
(1) 找出最优解的性质,并刻画其结构特征。
(2) 递归地定义最优解的值
(3) 以自底而上的方式计算出最优值
(4) 根据计算最优值时得到的信息,构造一个最优解。
对于一个给定的问题,若具有以下两个性质,则可以考虑用动态规划法来求解。
(1) 最优子结构。如果一个问题的最优解中包含了其子问题的最优解,就说该问题具有最优子结构。当一个问题具有最优子结构时,提示我们动态规划法可能会适用,但是此时贪心策略可能也是适用的。
(2) 重叠子问题。指用来解原问题的递归算法可反复地解同样的子问题,而不是总在产生新的子问题。即当一个递归算法不断地调用同一个问题时,就说该问题包含重叠子问题。此时若用分治法递归求解,则每次遇到子问题都会视为新问题,会极大地降低算法的效率,而动态规划法总是充分利用重叠子问题,对于每个子问题仅计算一次,把解保存在一个在需要时就可以查看的表中,而每次查表的时间为常数。
问题:有n个物品,第i个物品价值为vi,重量为wi,其中vi和wi均为非负数,背包的容量为W,W为非负数。现需要考虑如何选择装入背包的物品,使装入背包的物品总价值最大。该问题以形式化描述如下:
目标函数为 :
![](https://img-blog.csdn.net/20141114135539102?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvaHVzdF96aG91eQ==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center)
约束条件为:
![](https://img-blog.csdn.net/20141114135558212?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvaHVzdF96aG91eQ==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center)
满足约束条件的任一集合(x1,x2,...,xn)是问题的一个可行解,问题的目标是要求问题的一个最优解。考虑一个实例,假设n=5,W=17, 每个物品的价值和重量如表9-1所示。可将物品1,2和5装入背包,背包未满,获得价值22,此时问题解为你(1,1,0,0,1)。也可以将物品4和5装入背包,背包装满,获得价值24,此时解为(0,0,0,1,1)。
![](https://img-blog.csdn.net/20141114135646260?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvaHVzdF96aG91eQ==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center)
下面根据动态规划的4个步骤求解该问题。
(1) 刻画0-1背包问题的最优解的结构。
可以将背包问题的求解过程看作是进行一系列的决策过程,即决定哪些物品应该放入背包,哪些物品不放入背包。如果一个问题的最优解包含了物品n,即xn=1,那么其余x1,x2,...,x(n-1)一定构成子问题1,2,...,n-1在容量W-wn时的最优解。如果这个最优解不包含物品n,即xn=0,那么其余x1,x2,...,x(n-1)一定构成子问题1,2,...,n-1在容量W时的最优解。
(2)递归定义最优解的值
根据上述分析的最优解的结构递归地定义问题最优解。设c[i,w]表示背包容量为w时,i个物品导致的最优解的总价值,得到下式。显然要求c[n,w]。
![](https://img-blog.csdn.net/20141114135714840?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvaHVzdF96aG91eQ==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center)
下面是写的程序:在设置物品的重量的价值时将数组下标为0的元素设置为了0,只是为了统一数组下标,也可不这么做,但是在后面这一步时,要注意下标需要减1。
程序如下
#include<stdio.h>#include<stdlib.h>
int Total = 17;
int max(int i , int j){ return (i>j)?i:j;}
int main(){ int v[6] = {0,4,5,10,11,13}; int w[6] = {0,3,4,7,8,9}; int c[6][18] = {0};
int i,j; for(i = 0 ; i < 6 ; ++i) c[i][0] = 0; for(j = 0 ; j < 18 ; ++j) c[0][j] = 0;
for(i = 1 ; i < 6 ; ++i) for(j = 1 ; j <= Total ; ++j) { if(w[i] > j) c[i][j] = c[i-1][j]; else c[i][j] = max(c[i-1][j - w[i]] + v[i],c[i-1][j]); } for(i = 0 ; i < 6 ; ++i) { for(j = 0 ; j < 18 ; ++j) printf("%3d",c[i][j]); printf("\n"); }}
运行结果如下所示(输出数组):
重要的是要理解动态规划的思想
找到的一些关于01背包问题的求解思想:
首先介绍一下动态规划...
设计一个动态规划算法,通常可按照以下几个步骤进行:
(1) 找出最优解的性质,并刻画其结构特征。
(2) 递归地定义最优解的值
(3) 以自底而上的方式计算出最优值
(4) 根据计算最优值时得到的信息,构造一个最优解。
对于一个给定的问题,若具有以下两个性质,则可以考虑用动态规划法来求解。
(1) 最优子结构。如果一个问题的最优解中包含了其子问题的最优解,就说该问题具有最优子结构。当一个问题具有最优子结构时,提示我们动态规划法可能会适用,但是此时贪心策略可能也是适用的。
(2) 重叠子问题。指用来解原问题的递归算法可反复地解同样的子问题,而不是总在产生新的子问题。即当一个递归算法不断地调用同一个问题时,就说该问题包含重叠子问题。此时若用分治法递归求解,则每次遇到子问题都会视为新问题,会极大地降低算法的效率,而动态规划法总是充分利用重叠子问题,对于每个子问题仅计算一次,把解保存在一个在需要时就可以查看的表中,而每次查表的时间为常数。
问题:有n个物品,第i个物品价值为vi,重量为wi,其中vi和wi均为非负数,背包的容量为W,W为非负数。现需要考虑如何选择装入背包的物品,使装入背包的物品总价值最大。该问题以形式化描述如下:
目标函数为 :
约束条件为:
满足约束条件的任一集合(x1,x2,...,xn)是问题的一个可行解,问题的目标是要求问题的一个最优解。考虑一个实例,假设n=5,W=17, 每个物品的价值和重量如表9-1所示。可将物品1,2和5装入背包,背包未满,获得价值22,此时问题解为你(1,1,0,0,1)。也可以将物品4和5装入背包,背包装满,获得价值24,此时解为(0,0,0,1,1)。
下面根据动态规划的4个步骤求解该问题。
(1) 刻画0-1背包问题的最优解的结构。
可以将背包问题的求解过程看作是进行一系列的决策过程,即决定哪些物品应该放入背包,哪些物品不放入背包。如果一个问题的最优解包含了物品n,即xn=1,那么其余x1,x2,...,x(n-1)一定构成子问题1,2,...,n-1在容量W-wn时的最优解。如果这个最优解不包含物品n,即xn=0,那么其余x1,x2,...,x(n-1)一定构成子问题1,2,...,n-1在容量W时的最优解。
(2)递归定义最优解的值
根据上述分析的最优解的结构递归地定义问题最优解。设c[i,w]表示背包容量为w时,i个物品导致的最优解的总价值,得到下式。显然要求c[n,w]。
下面是写的程序:在设置物品的重量的价值时将数组下标为0的元素设置为了0,只是为了统一数组下标,也可不这么做,但是在后面这一步时,要注意下标需要减1。
程序如下
#include<stdio.h>#include<stdlib.h>
int Total = 17;
int max(int i , int j){ return (i>j)?i:j;}
int main(){ int v[6] = {0,4,5,10,11,13}; int w[6] = {0,3,4,7,8,9}; int c[6][18] = {0};
int i,j; for(i = 0 ; i < 6 ; ++i) c[i][0] = 0; for(j = 0 ; j < 18 ; ++j) c[0][j] = 0;
for(i = 1 ; i < 6 ; ++i) for(j = 1 ; j <= Total ; ++j) { if(w[i] > j) c[i][j] = c[i-1][j]; else c[i][j] = max(c[i-1][j - w[i]] + v[i],c[i-1][j]); } for(i = 0 ; i < 6 ; ++i) { for(j = 0 ; j < 18 ; ++j) printf("%3d",c[i][j]); printf("\n"); }}
运行结果如下所示(输出数组):