新《计算机算法设计与分析》课程设计用分治法解决快速排序问题及用回溯法解决0-1背包问题一、课程设计目的:《计算机算法设计与分析》这门课程是一门实践性非常强的课程,要求我们能够将所学的算法应用到实际中,灵活解决实际问题。
通过这次课程设计,能够培养我们独立思考、综合分析与动手的能力,并能加深对课堂所学理论和概念的理解,可以训练我们算法设计的思维和培养算法的分析能力。
二、课程设计内容:1、分治法:(2)快速排序;2、回溯法:(2)图的着色。
三、概要设计:分治法—快速排序:分治法的基本思想是将一个规模为n 的问题分解为k个规模较小的子问题,这些子问题互相独立且与原问题相同。
递归地解这些子问题,然后将各个子问题的解合并得data[i]=data[j]; data[j]=temp;}data[p]=data[j];data[j]=n;return j;}void quick_sort(int data[],int p,int r){ if(p>=r)return;int q=partition(data,p,r);quick_sort(data,p,q-1); //对左半段排序quick_sort(data,q+1,r); //对右半段排序}int main(){int i,n,data[size];printf("请输入要排列的数目(<=20):");scanf("%d",&n);printf("请输入要排列的数列:\n");for(i=0;i<n;++i)scanf("%d",&data[i]);quick_sort(data,0,n-1);printf("排列后的数列为:\n");for(i=0;i<n;++i)printf( "%d ",data[i]);printf("\n");return 0;}运行结果如下:图1图5●回溯法—0-1背包问题●回溯法是一个既带有系统性又带有跳跃性的的搜索算法。
它在包含问题的所有解的解空间树中,按照深度优先的策略,从根结点出发搜索解空间树。
算法搜索至解空间树的任一结点时,总是先判断该结点是否肯定不包含问题的解。
如果肯定不包含,则跳过对以该结点为根的子树的系统搜索,逐层向其祖先结点回溯。
否则,进入该子树,继续按深度优先的策略进行搜索。
回溯法在用来求问题的所有解时,要回溯到根,且根结点的所有子树都已被搜索遍才结束。
而回溯法在用来求问题的任一解时,只要搜索到问题的一个解就可以结束。
这种以深度优先的方式系统地搜索问题的解的算法称为回溯法,它适用于解一些组合数较大的问题。
二、算法框架:1、问题的解空间:应用回溯法解问题时,首先应明确定义问题的解空间。
问题的解空间应到少包含问题的一个(最优)解。
2、回溯法的基本思想:确定了解空间的组织结构后,回溯法就从开始结点(根结点)出发,以深度优先的方式搜索整个解空间。
这个开始结点就成为一个活结点,同时也成为当前的扩展结点。
在当前的扩展结点处,搜索向纵深方向移至一个新结点。
这个新结点就成为一个新的活结点,并成为当前扩展结点。
如果在当前的扩展结点处不能再向纵深方向移动,则当前扩展结点就成为死结点。
换句话说,这个结点不再是一个活结点。
此时,应往回移动(回溯)至最近的一个活结点处,并使这个活结点成为当前的扩展结点。
回溯法即以这种工作方式递归地在解空间中搜索,直至找到所要求的解或解空间中已没有活结点时为止。
运用回溯法解题通常包含以下三个步骤:(1)针对所给问题,定义问题的解空间;(2)确定易于搜索的解空间结构;(3)以深度优先的方式搜索解空间,并且在搜索过程中用剪枝函数避免无效搜索;3、递归回溯:由于回溯法是对解空间的深度优先搜索,因此在一般情况下可用递归函数来实现回溯法。
1 问题描述对于0-1背包问题回溯法的一个实例,n=4,c=7,p=[9,10,7,4],w=[3,5,2,1].这4个物品的单位重量价值分别为[3,2,3,5,4].以物品为单位价值的递减序装入物品。
先装入物品4,然后装入物品3和1.装入这3个物品后,剩余的背包容量为1,只能装入0.2的物品2.由此可得到一个解为x=[1,0,2,1,1],其相应的价值为22.尽管这不是一个可行解,但可以证明其价值是最有大的上界。
因此,对于这个实例,最优值不超过22. 回溯法的基本思想是按深度优先策略,从根节点出发搜索解空间树,算法搜索至解空间的任一点时,先判断该结点是否包含问题的解,如果肯定不包含,则跳过以该结点为根的子树的搜索,逐层向其祖先结点回溯,否则,进入该子树,继续按深度优先进行搜索。
2 问题分析利用回溯法的设计思想来解决背包问题。
首先是将可供选择的物品的个数输入程序,将物品排成一列,计算总物品的体积s,然后输入背包的实际体积V,如果背包的体积小于0或者大于物品的总体积s,则判断输入的背包体积错误,否则开始顺序选取物品装入背包,假设已选取了前i 件物品之后背包还没有装满,则继续选取第i+1件物品,若该件物品"太大"不能装入,则弃之而继续选取下一件,直至背包装满为止。
但如果在剩余的物品中找不到合适的物品以填满背包,则说明"刚刚"装入背包的那件物品"不合适",应将它取出"弃之一边",继续再从"它之后"的物品中选取,如此重复,直至求得满足条件的解。
因为回溯求解的规则是"后进先出",所以要用到栈来存储符合条件的解,在存储过程中,利用数组来存储各个物品的体积,然后用深度优先的搜索方式求解,将符合条件的数组元素的下标存入栈里,最后得到符合条件的解并且实现输出。
3 建立数学模型0-l背包问题是子集选取问题。
一般情况下,0-1背包问题是NP难题。
0-1背包问题的解空间可用子集树表示。
解0-1背包问题的回溯法与装载问题的回溯法十分类似。
在搜索解空间树时,只要其左儿子结点是一个可行结点,搜索就进入其左子树。
当右子树有可能包含最优解时才进入右子树搜索。
否则将右子树剪去。
设r是当前剩余物品价值总和;cp是当前价值;bestp是当前最优价值。
当cp+r≤bestp时,可剪去右子树。
计算右子树中解的上界的更好方法是将剩余物品依其单位重量价值排序,然后依次装入物品,直至装不下时,再装入该物品的一部分而装满背包。
由此得到的价值是右子树中解的上界。
算法复杂性是算法运行所需要的计算机资源的量,需要时间资源的量称为时间复杂性,需要的空间资源的量称为空间复杂性。
这个量应该只依赖于算法要解的问题的规模、算法的输入和算法本身的函数。
如果分别用N 、I 和A 表示算法要解问题的规模、算法的输入和算法本身,而且用C 表示复杂性,那么,应该有C=F(N,I,A)。
一般把时间复杂性和空间复杂性分开,并分别用T 和S 来表示,则有: T=T(N,I)和S=S(N,I) 。
最坏情况下的时间复杂性:最好情况下的时间复杂性:平均情况下的时间复杂性:其中DN 是规模为N 的合法输入的集合;I*是DN 中使T(N, I*)达到Tmax(N)的合法输入; 是中使T(N, )达到Tmin(N)的合法输入;而P(I)是在算法的应用中出现输入I 的概率。
),(max max I N T (N)T N D I ∈=),(max 1I N e t k i i i D I N ∑=∈=),(*1I N e t k i i i ∑==),(*I N T =),(min min I N T (N)T N D I ∈=),(m in 1I N e t k i i i D I N ∑=∈=)~,(1I N e t k i i i ∑==)~,(I N T =∑∈=N D I I N T I P (N)T ),()(avg ∑∑∈==N D I k i i i I N e t I P ),()(14 算法设计使用栈作为该程序的数据结构,利用栈进行语法检查,以深度优先的搜索方式解空间,实现递归过程和函数的调用,在设计时还使用C语言的数组及其循环语言来实现程序。
运用回溯法解题,在搜索解空间树时,只要其左儿子节点是一个可行结点,搜索就进入左子树,在右子树中有可能包含最优解是才进入右子树搜索。
否则将右子树剪去。
具体步骤如下:1.针对所给问题,定义问题的解空间;2.确定易于搜索的解空间结3.以深度优先的方式搜索解空间,并且在搜索过程中用剪枝函数避免无效搜索;5 算法实现运用Pc机,C/C++程序学习系统软件编程C语言程序源代码核心源代码#include <iostream.h>struct bag{double w;double cw,cp,bestp;};struct goods{double wc,pc;int v;};void paixu(goods g[],int n){int i,j;double temp;for(i=0;i<n-1;i++){for(j=n-1;j>i;j--){if((g[j].pc/g[j].wc)>(g[j-1].pc/g[j-1].wc)){temp=g[j].pc;g[j].pc=g[j-1].pc;g[j-1].pc=temp;temp=g[j].wc;g[j].wc=g[j-1].wc;g[j-1].wc=temp;}}}}void backtruck(int i,int n,goods g[],bag *b) {if(i>=n&&b->cw<=b->w){if(b->bestp<b->cp)b->bestp=b->cp;}if(i<n&&b->cw<=b->w){b->cw=b->cw+g[i].wc;b->cp=b->cp+g[i].pc;backtruck(i+1,n,g,b);b->cw=b->cw-g[i].wc;b->cp=b->cp-g[i].pc;backtruck(i+1,n,g,b);}}void main(){goods g[100];bag b;int n,i;cout<<"输入背包容量:";cin>>b.w;cout<<"物品数量:";cin>>n;cout<<"输入物品的质量:";for(i=0;i<n;i++){g[i].v=1;cin>>g[i].wc;}cout<<"输入物品价值:";for(i=0;i<n;i++){g[i].v=1;cin>>g[i].pc;}paixu(g,n);b.bestp=0;b.cw=0;b.cp=0;backtruck(0,n,g,&b);cout<<"最优价值:";cout<<b.bestp<<endl;}6 测试分析1.当输入的背包体积V小于物品的总体积s 时:2.当输入的背包体积V大于物品的总体积s 时:当输入的背包体积小于等于可供选择的物品的总体积时,可以正确进行选择装入的物品,并且统计符合条件的解的方法;当输入的背包体积小于0或者大于可供选择的物品的总体积时,程序会自动提醒输入背包体积错误,这样就方便了程序的正确运行与解决实际问题的重要作用。