算法设计与分析第9章 近似算法

《算法设计与分析第9章 近似算法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第9章近似算法1第9章近似算法迄今为止，所有的NP完全问题都还没有多项式时间算法。对于这类问题，通常可采取以下几种解题策略。(1)只对问题的特殊实例求解(2)用动态规划法或分支限界法求解(3)用概率算法求解(4)只求近似解(5)用启发式方法求解本章主要讨论解NP完全问题的近似算法。29.1近似算法的性能若一个最优化问题的最优值为c*，求解该问题的一个近似算法求得的近似最优解相应的目标函数值为c，则将该近似算法的性能比定义为=。在通常情况下，该性能比是问题输入规模n的一个函数ρ(n)，即≤ρ(n)。该近似算法的相对误差定义为=。若对问题的输入规模n，有一函数ε(n)使得

2、≤ε(n)，则称ε(n)为该近似算法的相对误差界。近似算法的性能比ρ(n)与相对误差界ε(n)之间显然有如下关系：ε(n)≤ρ(n)-1。39.2顶点覆盖问题的近似算法问题描述：无向图G=(V,E)的顶点覆盖是它的顶点集V的一个子集V’V，使得若(u,v)是G的一条边，则v∈V’或u∈V’。顶点覆盖V’的大小是它所包含的顶点个数

3、V’

4、。VertexSetapproxVertexCover(Graphg){cset=；e1=g.e；while(e1!=){从e1中任取一条边(u,v)；cset=cset∪{u,v}；从e1中删去与u和v相关联的所有边；}return

5、c}Cset用来存储顶点覆盖中的各顶点。初始为空，不断从边集e1中选取一边(u,v)，将边的端点加入cset中，并将e1中已被u和v覆盖的边删去，直至cset已覆盖所有边。即e1为空。49.2顶点覆盖问题的近似算法图(a)～(e)说明了算法的运行过程及结果。(e)表示算法产生的近似最优顶点覆盖cset，它由顶点b,c,d,e,f,g所组成。(f)是图G的一个最小顶点覆盖，它只含有3个顶点：b,d和e。算法approxVertexCover的性能比为2。59.3旅行售货员问题近似算法问题描述：给定一个完全无向图G=(V,E)，其每一边(u,v)∈E有一非负整数费用c(u,v

6、)。要找出G的最小费用哈密顿回路。比如，费用函数c往往具有三角不等式性质，即对任意的3个顶点u,v,w∈V，有：c(u,w)≤c(u,v)+c(v,w)。当图G中的顶点就是平面上的点，任意2顶点间的费用就是这2点间的欧氏距离时，费用函数c就具有三角不等式性质。旅行售货员问题的一些特殊性质：69.3.1具有三角不等式性质的旅行售货员问题对于给定的无向图G，可以利用找图G的最小生成树的算法设计找近似最优的旅行售货员回路的算法。voidapproxTSP(Graphg){(1)选择g的任一顶点r；(2)用Prim算法找出带权图g的一棵以r为根的最小生成树T；(3)前序遍历树T得

7、到的顶点表L；(4)将r加到表L的末尾，按表L中顶点次序组成回路H，作为计算结果返回；}当费用函数满足三角不等式时，算法找出的旅行售货员回路的费用不会超过最优旅行售货员回路费用的2倍。79.3.1具有三角不等式性质的旅行售货员问题举例(b)表示找到的最小生成树T；(c)表示对T作前序遍历的次序；(d)表示L产生的哈密顿回路H；(e)是G的一个最小费用旅行售货员回路。89.3.2一般的旅行售货员问题在费用函数不一定满足三角不等式的一般情况下，不存在具有常数性能比的解TSP问题的多项式时间近似算法，除非P=NP。换句话说，若P≠NP，则对任意常数ρ>1，不存在性能比为ρ的解旅

8、行售货员问题的多项式时间近似算法。99.4集合覆盖问题的近似算法问题描述：给定一个完全无向图G=(V,E)，其每一边(u,v)∈E有一非负整数费用c(u,v)。要找出G的最小费用哈密顿回路。集合覆盖问题的一个实例〈X,F〉由一个有限集X及X的一个子集族F组成。子集族F覆盖了有限集X。也就是说X中每一元素至少属于F中的一个子集，即X=。对于F中的一个子集CF，若C中的X的子集覆盖了X，即X=，则称C覆盖了X。集合覆盖问题就是要找出F中覆盖X的最小子集C*，使得

9、C*

10、=min{

11、C

12、

13、CF且C覆盖X}109.4集合覆盖问题的近似算法集合覆盖问题举例：用12个黑点表示集合

14、X。F={S1,S2,S3,S4,S5,S6,}，如图所示。容易看出，对于这个例子，最小集合覆盖为：C={S3,S4,S5,}。119.4集合覆盖问题的近似算法集合覆盖问题近似算法——贪心算法算法的循环体最多执行min{

15、X

16、，

17、F

18、}次。而循环体内的计算显然可在O(

19、X

20、

21、F

22、)时间内完成。因此，算法的计算时间为O(

23、X

24、

25、F

26、min{

27、X

28、，

29、F

30、})。由此即知，该算法是一个多项式时间算法。SetgreedySetCover(X,F){U=X；C=；while(U!=){选择F中使

31、S∩U

32、最大的子集S；U=U-S；C