蛮力法0/1背包
蛮力法
蛮力法是一种简单直接解决问题的方法,常常直接基于问题的描述,所以蛮力法也是最容易应用的方法。
蛮力法所依赖 的基本技术是遍历,即采用一定的策略依次处理待求解问题的所有元素,从而找出问题的解。由于其需要依次穷举待处理的元素,因此蛮力法是一种典型的指数级时间算法。
问题
给定n个重量为{w1,w2,···,wn}、价值为{v1,v2,···,vn}的物品和一个容量为C的背包,0/1背包是一个求解这些物品中的一个最有价值的子集,并且能够装入到背包中。
应用实例
有n项可投资的项目,每个项目需要投入资金si,可获利润为vi,现有可用资金M,应选择那些项目来投资才能获得最大利润。
想法
用蛮力法解决0/1背包问题,需要考虑给定n个物品集合的所有子集,找出所有重量不超过背包容量的子集,计算每个可能子集的总价值,然后找到价值最大的子集。例如,给定4个物品的重量为{7,3,4,5},价值为{42,12,40,25},和一个容量为10的背包,下表为求解的过程。
伪代码
输入:重量{w1,w2,···,wn},价值{v1,v2,···,vn},容量C
输出:装入背包的物品编号
初始化最大价值maxValue=0;结果子集S=∅;对集合{1,2,···,n}的每一个子集T,执行如下操作:
2.1 初始化背包的价值value=0;背包的重量weight=0;
2.2 对子集T的每一个元素j
2.2.1 如果weight+wj<C,则weight=weight+wj;value=value+vj;
2.2.2 否则,转入步骤2考察下一个子集;
2.3 如果maxValue<value,则maxValue=value;S=T;输出S中的各元素;
源代码
#include<iostream>#include<cstdio>#include<cstdlib>//用于计算程序运算时间#include<ctime>using namespace std;#define N 100struct goods{int sign;//物品序号int wight;//物品重量int value;//物品价值};int n, maxValue, cv, cw, C;//物品数量,价值最大,当前价值,当前重量,背包容量int X[N],cx[N];//最终存储状态,当前存储状态struct goods goods[N];int Force(int i){if(i > n-1){if(maxValue < cv && cw + goods[i].wight <= C){for(int k = 0;k < n;k++)X[k] = cx[k];//存储最优路径maxValue = cv;}return maxValue;}cw = cw + goods[i].wight;cv = cv + goods[i].value;cx[i] = 1;//装入背包Force(i+1);cw = cw-goods[i].wight;cv = cv-goods[i].value;cx[i] = 0;//不装入背包Force(i+1);return maxValue;}int main(){clock_t starttime,finish; //用于计算时间double time;starttime=clock();//用于计算程序运行时间printf("物品种类n:");scanf("%d",&n);printf("背包容量C:");scanf("%d",&C);for(int i = 0; i < n; i++){printf("物品%d的重量w[%d]及其价值v[%d]:",i+1,i+1,i+1);scanf("%d%d",&goods[i].wight,&goods[i].value);}int sum1 = Force(0);printf("蛮力法求解0/1背包问题:\nX=[");for(int i = 0; i < n; i++){cout << X[i]<<" ";}printf("] 装入总价值%d\n",sum1);finish=clock();time = (double)(finish-starttime)/1000;//时间单位为秒starttime=clock();//用于计算程序运行时间cout<<"time is:"<<time<<endl;return 0;}
算法分析
对于一个具有n个元素的集合,其子集数量是2n,所以无论生成子集的算法效率有多高,蛮力法求解0/1背包都会是一个Ω(2n)的算法。