0027算法笔记回溯法回溯法与装载问题.docx

1、0027算法笔记回溯法回溯法与装载问题1、回溯法 (1)描述:回溯法是一种选优搜索法，按选优条件向前搜索，以达到目标。但当探索到某一步时，发现原先选择并不优或达不到目标，就退回一步重新选择，这种走不通就退回再走的技术为回溯法。 (2)原理:回溯法在问题的解空间树中，按深度优先策略，从根结点出发搜索解空间树。算法搜索至解空间树的任意一点时，先判断该结点是否包含问题的解。如果肯定不包含，则跳过对该结点为根的子树的搜索，逐层向其祖先结点回溯；否则，进入该子树，继续按深度优先策略搜索。 回溯法的基本做法是搜索，或是一种组织得井井有条的，能避免不必要搜索的穷举式搜索法。这种方法适用于解一些组合数相当大的

2、问题。有许多问题，当需要找出它的解集或者要求回答什么解是满足某些约束条件的最佳解时，往往要使用回溯法。 (3)问题的解空间 问题的解向量：回溯法希望一个问题的解能够表示成一个n元式(x1,x2,xn)的形式。 显约束：对分量xi的取值限定。 隐约束：为满足问题的解而对不同分量之间施加的约束。 解空间：对于问题的一个实例，解向量满足显式约束条件的所有多元组，构成了该实例的一个解空间。 注意：同一个问题可以有多种表示，有些表示方法更简单，所需表示的状态空间更小（存储量少，搜索方法简单）。 例1：n=3的01 背包问题的回溯法搜索过程。W=16,15,15 p=45,25,25 C=30 例2：旅行

3、售货员问题。某售货员要到若干城市去推销商品，已知各城市之间的路程（旅费），他要选定一条从驻地出发，经过每个城市一遍，最后回到驻地的路线，使总的路程（总旅费）最小。 (4)生成问题状态的基本方法 扩展结点:一个正在产生儿子的结点称为扩展结点。 活结点:一个自身已生成但其儿子还没有全部生成的节点称做活结点。 死结点:一个所有儿子已经产生的结点称做死结点。 深度优先的问题状态生成法：如果对一个扩展结点R，一旦产生了它的一个儿子C，就把C当做新的扩展结点。在完成对子树C（以C为根的子树）的穷尽搜索之后，将R重新变成扩展结点，继续生成R的下一个儿子（如果存在）。 宽度优先的问题状态生成法：在一个扩展结点

4、变成死结点之前，它一直是扩展结点。 回溯法：为了避免生成那些不可能产生最佳解的问题状态，要不断地利用限界函数(bounding function)来处死那些实际上不可能产生所需解的活结点，以减少问题的计算量。具有限界函数的深度优先生成法称为回溯法。 (5)回溯法的基本思想 基本思想： 用回溯法解题的一个显著特征是在搜索过程中动态产生问题的解空间。在任何时刻，算法只保存从根结点到当前扩展结点的路径。如果解空间树中从根结点到叶结点的最长路径的长度为h(n)，则回溯法所需的计算空间通常为O(h(n)。而显式地存储整个解空间则需要O(2h(n)或O(h(n)!)内存空间。 解题步骤： 1)针对所给问题

5、，定义问题的解空间； 2)确定易于搜索的解空间结构； 3)以深度优先方式搜索解空间，并在搜索过程中用剪枝函数避免无效搜索。 常用剪枝函数：用约束函数在扩展结点处剪去不满足约束的子树；用限界函数剪去得不到最优解的子树。 递归回溯： 回溯法对解空间作深度优先搜索，因此，在一般情况下用递归方法实现回溯法。cppview plaincopy1. voidbacktrack(intt)2. 3. if(tn)4. output(x);/已到叶子结点，输出结果5. else6. for(inti=f(n,t);i0)5. if(f(n,t)=g(n,t)6. for(inti=f(n,t);in)outp

6、ut(x);4. else5. for(inti=0;in)output(x);4. else5. for(inti=t;ic1时，以结点z为根的子树中所有结点都不满足约束条件，因而该子树中的解均为不可行解，故可将该子树剪去。（该约束函数去除不可行解，得到所有可行解）。 可以引入一个上界函数，用于剪去不含最优解的子树，从而改进算法在平均情况下的运行效率。设z是解空间树第i层上的当前扩展结点。cw是当前载重量；bestw是当前最优载重量；r是剩余集装箱的重量，即r=。定义上界函数为cw+r。在以z为根的子树中任一叶结点所相应的载重量均不超过cw+r。因此，当cw+r=bestw时，可将z的右子树

7、剪去。 递归回溯具体代码如下：cppview plaincopy1. #includestdafx.h2. #include3. usingnamespacestd;4. 5. template6. classLoading7. 8. /friendTypeMaxLoading(Type,Type,int,int);9. /private:10. public:11. voidBacktrack(inti);12. intn,/集装箱数13. *x,/当前解14. *bestx;/当前最优解15. Type*w,/集装箱重量数组16. c,/第一艘轮船的载重量17. cw,/当前载重量18.

8、bestw,/当前最优载重量19. r;/剩余集装箱重量20. ;21. 22. template23. voidLoading:Backtrack(inti);24. 25. template26. TypeMaxLoading(Typew,Typec,intn,intbestx);27. 28. intmain()29. 30. intn=3,m;31. intc=50,c2=50;32. 33. intw4=0,10,40,40;34. intbestx4;35. 36. m=MaxLoading(w,c,n,bestx);37. 38. cout轮船的载重量分别为：endl;39. c

9、outc(1)=c,c(2)=c2endl;40. 41. cout待装集装箱重量分别为：endl;42. coutw(i)=;43. for(inti=1;i=n;i+)44. 45. coutwi;46. 47. coutendl;48. 49. cout回溯选择结果为：endl;50. coutm(1)=mendl;51. coutx(i)=;52. 53. for(inti=1;i=n;i+)54. 55. coutbestxi;56. 57. coutendl;58. 59. intm2=0;60. for(intj=1;j=n;j+)61. 62. m2=m2+wj*(1-best

10、xj);63. 64. coutm(2)=m2c2)67. 68. cout因为m(2)大于c(2),所以原问题无解！endl;69. 70. return0;71. 72. 73. template74. voidLoading:Backtrack(inti)/搜索第i层结点75. 76. if(in)/到达叶结点77. 78. if(cwbestw)79. 80. for(intj=1;j=n;j+)81. 82. bestxj=xj;/更新最优解83. bestw=cw;84. 85. 86. return;87. 88. 89. r-=wi;90. if(cw+wibestw)99.

11、100. xi=0;/搜索右子树101. Backtrack(i+1);102. 103. r+=wi;104. 105. 106. template107. TypeMaxLoading(Typew,Typec,intn,intbestx)/返回最优载重量108. 109. LoadingX;110. /初始化X111. X.x=newintn+1;112. X.w=w;113. X.c=c;114. X.n=n;115. X.bestx=bestx;116. X.bestw=0;117. X.cw=0;118. /初始化r119. X.r=0;120. 121. for(inti=1;i=

12、n;i+)122. 123. X.r+=wi;124. 125. 126. X.Backtrack(1);127. deleteX.x;128. returnX.bestw;129. 迭代回溯具体代码如下：cppview plaincopy1. #includestdafx.h2. #include3. usingnamespacestd;4. 5. template6. TypeMaxLoading(Typew,Typec,intn,intbestx);7. 8. intmain()9. 10. intn=3,m;11. intc=50,c2=50;12. intw4=0,10,40,40;

13、13. intbestx4;14. 15. m=MaxLoading(w,c,n,bestx);16. 17. cout轮船的载重量分别为：endl;18. coutc(1)=c,c(2)=c2endl;19. 20. cout待装集装箱重量分别为：endl;21. coutw(i)=;22. for(inti=1;i=n;i+)23. 24. coutwi;25. 26. coutendl;27. 28. cout回溯选择结果为：endl;29. coutm(1)=mendl;30. coutx(i)=;31. 32. for(inti=1;i=n;i+)33. 34. coutbestxi

14、;35. 36. coutendl;37. 38. intm2=0;39. for(intj=1;j=n;j+)40. 41. m2=m2+wj*(1-bestxj);42. 43. coutm(2)=m2c2)46. 47. cout因为m(2)大于c(2),所以原问题无解！endl;48. 49. return0;50. 51. 52. 53. template54. TypeMaxLoading(Typew,Typec,intn,intbestx)/迭代回溯法，返回最优载重量及其相应解，初始化根结点55. 56. inti=1;/当前层，x1:i-1为当前路径57. int*x=newi

15、ntn+1;58. 59. Typebestw=0,/当前最优载重量60. cw=0,/当前载重量61. r=0;/剩余集装箱重量62. 63. for(intj=1;j=n;j+)64. 65. r+=wj;66. 67. 68. while(true)/搜索子树69. 70. while(i=n&cw+win)/到达叶结点79. 80. for(intj=1;j=n;j+)81. 82. bestxj=xj;83. 84. bestw=cw;85. 86. else/进入右子树87. 88. r-=wi;89. xi=0;i+;90. 91. while(cw+r0&!xi)95. 96. r+=wi;97. i-;98. 99. /从右子树返回100. if(i=0)101. 102. deletex;103. returnbestw;104. 105. xi=0;106. cw-=wi;