版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、一、实验内容:分别用蛮力法、动态规划法、回溯法和分支限界法求解0/1背包问题。注:0/1背包问题:给定种物品和一个容量为的背包,物品的重量是,其价值为,背包问题是如何使选择装入背包内的物品,使得装入背包中的物品的总价值最大。其中,每种物品只有全部装入背包或不装入背包两种选择。二、所用算法的基本思想及复杂度分析:1.蛮力法求解0/1背包问题:1)基本思想:对于有n种可选物品的0/1背包问题,其解空间由长度为n的0-1向量组成,可用子集数表示。在搜索解空间树时,深度优先遍历,搜索每一个结点,无论是否可能产生最优解,都遍历至叶子结点,记录每次得到的装入总价值,然后记录遍历过的最大价值。2)代码:#i
2、nclude<iostream>#include<algorithm>using namespace std;#define N 100/最多可能物体数struct goods/物品结构体int sign;/物品序号int w;/物品重量int p;/物品价值aN;bool m(goods a,goods b)return (a.p/a.w)>(b.p/b.w);int max(int a,int b)return a<b?b:a;int n,C,bestP=0,cp=0,cw=0;int XN,cxN;/*蛮力法求解0/1背包问题*/int Force(i
3、nt i)if(i>n-1)if(bestP<cp&&cw+ai.w<=C)for (int k=0;k<n;k+)Xk=cxk;/存储最优路径bestP=cp;return bestP;cw=cw+ai.w;cp=cp+ai.p;cxi=1;/装入背包Force(i+1);cw=cw-ai.w;cp=cp-ai.p;cxi=0;/不装入背包Force(i+1);return bestP;int KnapSack1(int n,goods a,int C,int x)Force(0);return bestP;int main()goods bN;pri
4、ntf("物品种数n: ");scanf("%d",&n);/输入物品种数printf("背包容量C: ");scanf("%d",&C);/输入背包容量for (int i=0;i<n;i+)/输入物品i的重量w及其价值vprintf("物品%d的重量w%d及其价值v%d: ",i+1,i+1,i+1);scanf("%d%d",&ai.w,&ai.p);bi=ai;int sum1=KnapSack1(n,a,C,X);/调用蛮力法求
5、0/1背包问题printf("蛮力法求解0/1背包问题:nX= ");for(i=0;i<n;i+)cout<<Xi<<" "/输出所求Xn矩阵printf("装入总价值%dn",sum1);bestP=0,cp=0,cw=0;/恢复初始化3)复杂度分析:蛮力法求解0/1背包问题的时间复杂度为:。2.动态规划法求解0/1背包问题:1)基本思想:令表示在前个物品中能够装入容量为的背包中的物品的最大值,则可以得到如下动态函数:按照下述方法来划分阶段:第一阶段,只装入前1个物品,确定在各种情况下的背包能够得到的
6、最大价值;第二阶段,只装入前2个物品,确定在各种情况下的背包能够得到的最大价值;以此类推,直到第个阶段。最后,便是在容量为的背包中装入个物品时取得的最大价值。 2)代码: #include<iostream> #include<algorithm> using namespace std; #define N 100/最多可能物体数 struct goods/物品结构体 int sign;/物品序号 int w;/物品重量 int p;/物品价值 aN; bool m(goods a,goods b) return (a.p/a.w)>(b.p/b.w); int
7、 max(int a,int b) return a<b?b:a; int n,C,bestP=0,cp=0,cw=0; int XN,cxN; int KnapSack2(int n,goods a,int C,int x) int VN10*N; for(int i=0;i<=n;i+)/初始化第0列 Vi0=0; for(int j=0;j<=C;j+)/初始化第0行 V0j=0; for(i=1;i<=n;i+)/计算第i行,进行第i次迭代 for(j=1;j<=C;j+) if(j<ai-1.w) Vij=Vi-1j; else Vij=max(V
8、i-1j,Vi-1j-ai-1.w+ai-1.p); j=C;/求装入背包的物品 for (i=n;i>0;i-) if (Vij>Vi-1j) xi-1=1; j=j-ai-1.w; elsexi-1=0; return VnC;/返回背包取得的最大价值 int main() goods bN; printf("物品种数n: "); scanf("%d",&n);/输入物品种数 printf("背包容量C: "); scanf("%d",&C);/输入背包容量 for (int i=0
9、;i<n;i+)/输入物品i的重量w及其价值v printf("物品%d的重量w%d及其价值v%d: ",i+1,i+1,i+1); scanf("%d%d",&ai.w,&ai.p); bi=ai; int sum2=KnapSack2(n,a,C,X);/调用动态规划法求0/1背包问题 printf("动态规划法求解0/1背包问题:nX= "); for(i=0;i<n;i+) cout<<Xi<<" "/输出所求Xn矩阵 printf("装入总价值
10、%dn",sum2); for (i=0;i<n;i+) ai=bi; /恢复aN顺序 3)复杂度分析:动态规划法求解0/1背包问题的时间复杂度为:。3.回溯法求解0/1背包问题:1)基本思想:回溯法:为了避免生成那些不可能产生最佳解的问题状态,要不断地利用限界函数(bounding function)来处死那些实际上不可能产生所需解的活结点,以减少问题的计算量。这种具有限界函数的深度优先生成法称为回溯法。对于有n种可选物品的0/1背包问题,其解空间由长度为n的0-1向量组成,可用子集数表示。在搜索解空间树时,只要其左儿子结点是一个可行结点,搜索就进入左子树。当右子树中有可能包
11、含最优解时就进入右子树搜索。 2)代码: #include<iostream> #include<algorithm> using namespace std; #define N 100/最多可能物体数 struct goods/物品结构体 int sign;/物品序号 int w;/物品重量 int p;/物品价值 aN; bool m(goods a,goods b) return (a.p/a.w)>(b.p/b.w); int max(int a,int b) return a<b?b:a; int n,C,bestP=0,cp=0,cw=0; i
12、nt XN,cxN; int BackTrack(int i) if(i>n-1) if(bestP<cp) for (int k=0;k<n;k+)Xk=cxk;/存储最优路径 bestP=cp; return bestP; if(cw+ai.w<=C)/进入左子树 cw=cw+ai.w; cp=cp+ai.p; cxai.sign=1;/装入背包 BackTrack(i+1); cw=cw-ai.w; cp=cp-ai.p;/回溯,进入右子树 cxai.sign=0;/不装入背包 BackTrack(i+1); return bestP; int KnapSack3
13、(int n,goods a,int C,int x) for(int i=0;i<n;i+) xi=0; ai.sign=i; sort(a,a+n,m);/将各物品按单位重量价值降序排列 BackTrack(0); return bestP; int main() goods bN; printf("物品种数n: "); scanf("%d",&n);/输入物品种数 printf("背包容量C: "); scanf("%d",&C);/输入背包容量 for (int i=0;i<n;
14、i+)/输入物品i的重量w及其价值v printf("物品%d的重量w%d及其价值v%d: ",i+1,i+1,i+1); scanf("%d%d",&ai.w,&ai.p); bi=ai; int sum3=KnapSack3(n,a,C,X);/调用回溯法求0/1背包问题 printf("回溯法求解0/1背包问题:nX= "); for(i=0;i<n;i+) cout<<Xi<<" "/输出所求Xn矩阵 printf("装入总价值%dn",su
15、m3); for (i=0;i<n;i+) ai=bi; /恢复aN顺序3)复杂度分析:最不理想的情况下,回溯法求解0/1背包问题的时间复杂度为:。由于其对蛮力法进行优化后的算法,其复杂度一般比蛮力法要小。空间复杂度:有个物品,即最多递归层,存储物品信息就是一个一维数组,即回溯法求解0/1背包问题的空间复杂度为。4.分支限界法求解背包问题:1)基本思想:分支限界法类似于回溯法,也是在问题的解空间上搜索问题解的算法。一般情况下,分支限界法与回溯法的求解目标不同。回溯法的求解目标是找出解空间中满足约束条件的所有解,而分支限界法的求解目标则是找出满足约束条件的一个解,或是在满足约束条件的解中找
16、出使某一目标函数值达到极大或极小的解,即在某种意义下的最优解。首先,要对输入数据进行预处理,将各物品依其单位重量价值从大到小进行排列。在下面描述的优先队列分支限界法中,节点的优先级由已装袋的物品价值加上剩下的最大单位重量价值的物品装满剩余容量的价值和。算法首先检查当前扩展结点的左儿子结点的可行性。如果该左儿子结点是可行结点,则将它加入到子集树和活结点优先队列中。当前扩展结点的右儿子结点一定是可行结点,仅当右儿子结点满足上界约束时才将它加入子集树和活结点优先队列。当扩展到叶节点时为问题的最优值。2)代码:#include<iostream>#include<algorithm&
17、gt;using namespace std;#define N 100/最多可能物体数struct goods/物品结构体int sign;/物品序号int w;/物品重量int p;/物品价值aN;bool m(goods a,goods b)return (a.p/a.w)>(b.p/b.w);int max(int a,int b)return a<b?b:a;int n,C,bestP=0,cp=0,cw=0;int XN,cxN;struct KNAPNODE/状态结构体bool s1N; /当前放入物体int k;/搜索深度int b;/价值上界int w;/物体重量
18、int p;/物体价值;struct HEAP/堆元素结构体KNAPNODE *p;/结点数据int b; /所指结点的上界;/交换两个堆元素void swap(HEAP &a, HEAP &b)HEAP temp = a;a = b;b = temp;/堆中元素上移void mov_up(HEAP H, int i)bool done = false;if(i!=1)while(!done && i!=1)if(Hi.b>Hi/2.b)swap(Hi, Hi/2);elsedone = true;i = i/2;/堆中元素下移void mov_down(
19、HEAP H, int n, int i)bool done = false;if(2*i)<=n)while(!done && (i = 2*i) <= n)if(i+1<=n && Hi+1.b > Hi.b)i+;if(Hi/2.b<Hi.b)swap(Hi/2, Hi);elsedone = true;/往堆中插入结点void insert(HEAP H, HEAP x, int &n)n+;Hn = x;mov_up(H,n);/删除堆中结点void del(HEAP H, int &n, int i)HE
20、AP x, y;x = Hi; y = Hn;n -;if(i<=n)Hi = y;if(y.b>=x.b)mov_up(H,i);elsemov_down(H, n, i);/获得堆顶元素并删除HEAP del_top(HEAP H, int &n)HEAP x = H1;del(H, n, 1);return x;/计算分支节点的上界void bound( KNAPNODE* node, int M, goods a, int n)int i = node->k;float w = node->w;float p = node->p;if(node-&
21、gt;w>M) / 物体重量超过背包载重量node->b = 0; / 上界置为0elsewhile(w+ai.w<=M)&&(i<n) w += ai.w; / 计算背包已装入载重p += ai+.p; / 计算背包已装入价值if(i<n)node->b = p + (M - w)*ai.p/ai.w;elsenode -> b = p;/用分支限界法实现0/1背包问题int KnapSack4(int n,goods a,int C, int X)int i, k = 0; / 堆中元素个数的计数器初始化为0int v;KNAPNO
22、DE *xnode, *ynode, *znode;HEAP x, y, z, *heap;heap = new HEAPn*n; / 分配堆的存储空间for( i=0; i<n; i+)ai.sign=i; /记录物体的初始编号sort(a,a+n,m); / 对物体按照价值重量比排序xnode = new KNAPNODE; / 建立父亲结点for( i=0; i<n; i+) / 初始化结点xnode->s1i = false;xnode->k = xnode->w = xnode->p = 0;while(xnode->k<n) ynod
23、e = new KNAPNODE; / 建立结点y*ynode = *xnode; /结点x的数据复制到结点yynode->s1ynode->k = true; / 装入第k个物体ynode->w += aynode->k.w; / 背包中物体重量累计ynode->p += aynode->k.p; / 背包中物体价值累计ynode->k +; / 搜索深度+bound(ynode, C, a, n); / 计算结点y的上界y.b = ynode->b;y.p = ynode;insert(heap, y, k); /结点y按上界的值插入堆中zn
24、ode = new KNAPNODE; / 建立结点z*znode = *xnode; /结点x的数据复制到结点zznode->k+; / 搜索深度+bound(znode, C, a, n); /计算节点z的上界z.b = znode->b;z.p = znode;insert(heap, z, k); /结点z按上界的值插入堆中delete xnode;x = del_top(heap, k); /获得堆顶元素作为新的父亲结点xnode = x.p;v = xnode->p;for( i=0; i<n; i+) /取装入背包中物体在排序前的序号if(xnode->s1i)Xai.sign =1 ;elseXai.sign = 0;delete xnode;delete heap;return v; /返回背包中物体的价值/*测试以上算法的主函数*/int main()goods bN;printf("物品种数
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 促进产业技术工人职业发展的策略与实施路径
- 2024年度原材料卸车作业信息管理系统开发合同3篇
- 2024年环境监测与污染治理合同
- 2024年度冷库节能改造财政补贴申请合同2篇
- 内蒙古化工职业学院《证券投资学含证券交易实务模拟实验》2023-2024学年第一学期期末试卷
- 内蒙古鸿德文理学院《中学政治课教师口语训练》2023-2024学年第一学期期末试卷
- 2024年权益让渡与融资担保合同
- 2024版二手房买卖合同范本(含特殊户型改造约定)3篇
- 内蒙古工业大学《化学教材教法实验》2023-2024学年第一学期期末试卷
- 2024版二手房买卖按揭合同附带房屋交易配套设施保障范本3篇
- 改革开放史智慧树知到期末考试答案2024年
- 建筑工程钢结构焊接变形的控制措施
- 文化传媒公司招聘笔试题
- 零售业收银员操作培训
- 信息安全职业生涯规划
- 新能源技术的前沿研究与应用探索
- 装配式结构监理细则
- 工作无纸化总结汇报
- 初中美术八年级上册服装设计(全国一等奖)
- 导医接待与患者情绪管理
- 2024年医院口腔科护理带教计划
评论
0/150
提交评论