《第9章近似算法.ppt》由会员分享，可在线阅读，更多相关《第9章近似算法.ppt（16页珍藏版）》请在三一办公上搜索。
1、1,第9章 近似算法,2,第9章 近似算法,迄今为止，所有的NP完全问题都还没有多项式时间算法。对于这类问题，通常可采取以下几种解题策略。(1)只对问题的特殊实例求解(2)用动态规划法或分支限界法求解(3)用概率算法求解(4)只求近似解(5)用启发式方法求解 本章主要讨论解NP完全问题的近似算法。,3,9.1 近似算法的性能,若一个最优化问题的最优值为c*，求解该问题的一个近似算法求得的近似最优解相应的目标函数值为c，则将该近似算法的性能比定义为=。在通常情况下，该性能比是问题输入规模n的一个函数(n)，即(n)。,该近似算法的相对误差定义为=。若对问题的输入规模n，有一函数(n)使得(n)，
2、则称(n)为该近似算法的相对误差界。近似算法的性能比(n)与相对误差界(n)之间显然有如下关系：(n)(n)-1。,4,9.2 顶点覆盖问题的近似算法,问题描述：无向图G=(V,E)的顶点覆盖是它的顶点集V的一个子集VV，使得若(u,v)是G的一条边，则vV或uV。顶点覆盖V的大小是它所包含的顶点个数|V|。,VertexSet approxVertexCover(Graph g)cset=；e1=g.e；while(e1!=)从e1中任取一条边(u,v)；cset=csetu,v；从e1中删去与u和v相关联的所有边；return c,Cset用来存储顶点覆盖中的各顶点。初始为空，不断从边集e
3、1中选取一边(u,v)，将边的端点加入cset中，并将e1中已被u和v覆盖的边删去，直至cset已覆盖所有边。即e1为空。,5,9.2 顶点覆盖问题的近似算法,图(a)(e)说明了算法的运行过程及结果。(e)表示算法产生的近似最优顶点覆盖cset，它由顶点b,c,d,e,f,g所组成。(f)是图G的一个最小顶点覆盖，它只含有3个顶点：b,d和e。,算法approxVertexCover的性能比为2。,6,9.3 旅行售货员问题近似算法,问题描述：给定一个完全无向图G=(V,E)，其每一边(u,v)E有一非负整数费用c(u,v)。要找出G的最小费用哈密顿回路。,比如，费用函数c往往具有三角不等式
4、性质，即对任意的3个顶点u,v,wV，有：c(u,w)c(u,v)+c(v,w)。当图G中的顶点就是平面上的点，任意2顶点间的费用就是这2点间的欧氏距离时，费用函数c就具有三角不等式性质。,旅行售货员问题的一些特殊性质：,7,9.3.1 具有三角不等式性质的旅行售货员问题,对于给定的无向图G，可以利用找图G的最小生成树的算法设计找近似最优的旅行售货员回路的算法。,void approxTSP(Graph g)(1)选择g的任一顶点r；(2)用Prim算法找出带权图g的一棵以r为根的最小生成树T；(3)前序遍历树T得到的顶点表L；(4)将r加到表L的末尾，按表L中顶点次序组成回路H，作为计算结果
5、返回；,当费用函数满足三角不等式时，算法找出的旅行售货员回路的费用不会超过最优旅行售货员回路费用的2倍。,8,9.3.1 具有三角不等式性质的旅行售货员问题举例,(b)表示找到的最小生成树T；(c)表示对T作前序遍历的次序；(d)表示L产生的哈密顿回路H；(e)是G的一个最小费用旅行售货员回路。,9,9.3.2 一般的旅行售货员问题,在费用函数不一定满足三角不等式的一般情况下，不存在具有常数性能比的解TSP问题的多项式时间近似算法，除非P=NP。换句话说，若PNP，则对任意常数1，不存在性能比为的解旅行售货员问题的多项式时间近似算法。,10,9.4 集合覆盖问题的近似算法,集合覆盖问题的一个实
6、例X,F由一个有限集X及X的一个子集族F组成。子集族F覆盖了有限集X。也就是说X中每一元素至少属于F中的一个子集，即X=。对于F中的一个子集CF，若C中的X的子集覆盖了X，即X=，则称C覆盖了X。集合覆盖问题就是要找出F中覆盖X的最小子集C*，使得|C*|=min|C|CF且C覆盖X,11,9.4 集合覆盖问题的近似算法,集合覆盖问题举例：,用12个黑点表示集合X。F=S1,S2,S3,S4,S5,S6,，如图所示。容易看出，对于这个例子，最小集合覆盖为：C=S3,S4,S5,。,12,9.4 集合覆盖问题的近似算法,集合覆盖问题近似算法贪心算法,算法的循环体最多执行min|X|，|F|次。而
7、循环体内的计算显然可在O(|X|F|)时间内完成。因此，算法的计算时间为O(|X|F|min|X|，|F|)。由此即知，该算法是一个多项式时间算法。,Set greedySetCover(X,F)U=X；C=；while(U!=)选择F中使|SU|最大的子集S；U=U-S；C=CS；return C；,13,9.5 子集合问题的近似算法,问题描述：设子集和问题的一个实例为S,t。其中，S=x1，x2，xn是一个正整数的集合，t是一个正整数。子集和问题判定是否存在S的一个子集S1，使得。,14,9.5.1 子集合问题的指数时间算法,int exactSubsetSum(S,t)int n=|S|
8、；L0=0；for(int i=1；i=n；i+)Li=mergeLists(Li-1,Li-1+Si)；删去Li中超过t的元素；return max(Ln)；,算法以集合S=x1，x2，xn和目标值t作为输入。算法中用到将2个有序表L1和L2合并成为一个新的有序表的算法mergeLists(L1,L2)。,15,9.5.2 子集合问题的完全多项式时间近似格式,基于算法exactSubsetSum，通过对表Li作适当的修整建立一个子集和问题的完全多项式时间近似格式。,在对表Li进行修整时，用到一个修整参数，01。用参数修整一个表L是指从L中删去尽可能多的元素，使得每一个从L中删去的元素y，都有
9、一个修整后的表L1中的元素z满足(1-)yzy。可以将z看作是被删去元素y在修整后的新表L1中的代表。,举例：若=0.1，且L=10,11,12,15,20,21,22,23,24,29，则用对L进行修整后得到L1=10，12，15，20，23，29。其中被删去的数11由10来代表，21和22由20来代表，24由23来代表。,16,9.5.2 子集合问题的完全多项式时间近似格式,对有序表L修整算法,List trim(L,)int m=|L|；L1=L1；int last=L1；for(int i=2；i=m；i+)if(last(1-)*Li)将Li加入表L1的尾部；last=Li；return L1；,子集和问题近似格式,int approxSubsetSum(S,t,)n=|S|；L0=0；for(int i=1；i=n；i+)Li=Merge-Lists(Li-1,Li-1+Si)；Li=Trim(Li,/n)；删去Li中超过t的元素；return max(Ln)；,
