算法论文旅行商问题的求解方法动态规划法和贪心法讲解Word格式.docx

1、（4） 路径：（5） 路径：（6） 路径：从中，我们可以知道，路径（2）和（4）路径长度最短。我们还应注意到，图1中，有3对不同的路径，对每对路径来说，不同只是路径的方向，因此，可以将这个数量减半，则可能的解有（n-1）！/2个。这是一个非常大的数，随着n的增长，TSP问题的可能解也在迅速增长。如：一个10城市的TSP问题有大约有180,000个可能解。一个20城市的TSP问题有大约有60,000,000,000,000,000个可能解。 一个50城市的TSP问题有大约1062个可能解，而一个行星上也只有1021升水。因此，我们可以知道用蛮力法求解TSP问题，只能解决问题规模很小的实例。2.2

2、动态规划法2.2.1动态规划法的设计思想动态规划法将待求解问题分解成若干个相互重叠的子问题，每个子问题对应决策过程的一个阶段，一般来说，子问题的重叠关系表现在对给定问题求解的递推关系（也就是动态规划函数）中，将子问题的解求解一次并填入表中，当需要再次求解此子问题时，可以通过查表获得该子问题的解而不用再次求解，从而避免了大量重复计算。2.2.2TSP问题的动态规划函数假设从顶点i出发，令表示从顶点i出发经过中各个顶点一次且仅一次，最后回到出发点i的最短路径长度，开始时，于是，TSP问题的动态规划函数为：2.2.3算法讨论（1）for （i=1; iN; i+） /初始化第0列 di0=ci0;

3、（2）for （j=1; j -1; j+） for （i=1;n; i+） /依次进行第i次迭代 if （子集Vj中不包含i） 对Vj中的每个元素k，计算Vm = Vj-k;dij=min（cik+dkm）; （3）对V-1中的每一个元素k，计算Vm = V-1-k;d0 -1=min（c0k+dkm）; （4）输出最短路径长度d0 -1;2.3.4时间复杂性和蛮力法相比，动态规划法求解TSP问题，把原来的时间复杂性是O（n!）的排列问题，转化为组合问题，从而降低了算法的时间复杂性，但它仍需要指数时间。2.3贪心法2.3.1贪心法的设计思想贪心法在解决问题的策略上目光短浅，只根据当前已有的信

4、息就做出选择，而且一旦做出了选择，不管将来有什么结果，这个选择都不会改变。换言之，贪心法并不是从整体最优考虑，它所做出的选择只是在某种意义上的局部最优。这种局部最优选择并不总能获得整体最优解，但通常能获得近似最优解。2.3.2最近邻点策略求解TSP问题贪心法求解TSP问题的贪心策略是显然的，至少有两种贪心策略是合理的：最近邻点策略和最短链接策略。本文仅重点讨论最近邻点策略及其求解过程。最近邻点策略：从任意城市出发，每次在没有到过的城市中选择距离已选择的城市中最近的一个，直到经过了所有的城市，最后回到出发城市。2.3.3算法讨论1P= ; 2V=V-u0; u=u0; /从顶点u0出发3循环直到

5、集合P中包含n-1条边 3.1查找与顶点u邻接的最小代价边（u, v）并且v属于集合V； 3.2 P=P+（u, v）; 3.3 V=V-v; 3.4 u=v; /从顶点v出发继续求解但需注意，用最近邻点贪心策略求解TSP问题所得的结果不一定是最优解。当图中顶点个数较多并且各边的代价值分布比较均匀时，最近邻点策略可以给出较好的近似解，不过，这个近似解以何种程度近似于最优解，却难以保证。2.4分支限界法2.4.1分支限界法的设计思想 假设求解最大化问题，解向量为，其中，的取值范围为某个有穷集合，。在使用分支限界法搜索问题的解空间树时，首先根据限界函数估算目标函数的界down, up，然后从根结点

6、出发，扩展根结点的个孩子结点，从而构成分量的种可能的取值方式。对这个孩子结点分别估算可能取得的目标函数值，其含义是以该孩子结点为根的子树所可能取得的目标函数值不大于，也就是部分解应满足：本文本欲详细讨论该算法，但无奈在编程问题中，尚有问题有待解决，时间所限，不得已放弃。本人编程过程中所用算法思想与老师课上所教略有不同，在寻找下界时，是首先把每个结点所能到达的各个结点及其可能的路径算出来，并添加到PT表中，但最后，不知是何原因，在还有一个城市尚未加入时，PT表的添加出现了问题，思忖良久，仍未解决，时间所限，迫不得已，留待以后有时间再另行研究，本文就只给出动态规划法和贪心法的具体求解过程。3结论本

7、文主要重点讨论了动态规划法和贪心法求解TSP问题算法，并附录给出了相应程序。3.1动态规划法思想动态规划法中对于顶点元素生成的子集本文中用字符串形式存储，然后再用递归方法按照子集中元素个数从小到大开始赋值。因为后面元素个数较多的子集与前面比其元素个数少1的子集间有一定对应关系，所以用递归方式，可以简便很多。个人觉得这算本文的一大特色。另，在计算dij =min（cik+dkj-1）时，获得dkj-1的过程比较困难，运用字符串后，我们就可以首先找到指定字符，然后去掉该字符，返回剩余字符串，在与V逐个比较，找到与其相等的V中元素对应下标，此下标即为j-1；具体求解过程可参考附录源程序，有详细说明。

8、在求解最佳路径所经过城市顺序时，本文是通过边查找dij边记录路径的，这样可以省掉很多麻烦，另，路径也是采用字符串形式的数组，数组规模与存储城市间距离的c数组相同，由于很多元素均不需赋值，这样做可能会浪费内存空间，但是目前还没找到更好地求解方法。3.2贪心法思想贪心法中，由于贪心法相对动态规划法要简单很多，每次在查找最近城市时所得的顶点均为最后该法最佳路径所经过的城市编号，规模相对较小，容易确定，操作相对简单，所以本文用数组V存放最佳路径所经过的城市编号顺序相对来说方便很多。另外，本文用path整型数组存放所经路径的长度，最后相加即可得最短路径。3.3两者比较动态规划法相对贪心法来说虽然要精确些

9、，但代码相对繁杂很多，对时间和空间要求很多，仅适用于城市数量较小的情况。贪心法虽然比较简单，实现起来比较容易，但不是很精确，当图中顶点个数较多并且各边的代价值分布比较均匀时，贪心法可以给出较好的近似解，不过，这个近似解以何种程度近似于最优解，却难以保证。另外，动态规划法有一个明显的缺点，就是出发城市只能是第0个城市（城市从0开始编号），若出发城市改变，则必须以该城市为第0个城市顺序给其他城市编号，输入城市间距离。由于若出发城市任意，编码的难度大大增加，所以最后不得已放弃，但这大大地限制了程序的通用性。而对于贪心法，本文很好地避免了这个问题，一旦城市编号确定，可以从任意城市出发，这也是本文中贪心

10、法优于动态规划法的一点。3.4优点本文程序优点，各个子函数功能分隔很明显，没有大量集中在一个函数里面，而是分成了几个不同功能的小函数，这样程序可阅读性提高。另外，程序中有详细注释，程序中变量取名都是根据变量的性质和所代表的含义命名的，也相应提高了程序的可读性。对于动态规划法，城市个数可以在算法时间允许的范围内任意，于这点来说，通用性较好；对于贪心法，出发城市可以任意，城市个数也可以任意，通用性较好。4建议当城市个数较少时，用动态规划法求出最优解；当城市个数较多并且各边的代价值分布比较均匀时，贪心法可以给出较好的近似解。5参考文献（1）计算机算法分析与设计第二版，王晓东编著，电子工业出版社（2）

11、Java语言与面向对象程序设计（第2版）印旻、王行言编著，清华大学出版社（3）求解TSP算法，周康、强小利、同小军、许进，计算机工程与应用6附录6.1动态规划法6.1.1源代码package exp2;import java.util.Scanner;public class TSPDynamic String V;/顶点生成的子集，这里把每一个子集用一个字符串表示 int c;/顶点间距离 int d;/存放迭代结果 int N; /城市个数 String path;/用于存放每种选择下经过的城市 static int IFINITE = 99999;/无穷大距离 表示城市自己到达自己时，距

12、离无穷大，不作为考虑因素 /构造函数 public TSPDynamic（） initialC（）; initialV1（）; /初始化数组c，即顶点间距离 public void initialC（） Scanner in = new Scanner（ System.in ）; System.out.println（请输入城市个数： （注意根据实际情况城市个数不可小于1！）; N = in.nextInt（）; if（N = 1）不符合要求，请认真核对！ System.exit（0）;/输入错误，结束！请输入城市相邻城市间距离（城市从0开始编号，且出发城市为第0个城市！）： c = new

13、intNN;/为c分配空间 for（int i = 0 ; i N ; i +） for（int j = 0 ; j （int）Math.pow（2,N-1） - 1）return;/如果全部顶点已初始化完成，则返回。 if（m = 0）Vm + = ;/初始化出发顶点，即V0 else int i = m - 1; while（i = 0 & Vi.length（） = len）/找与最后一个初始化的Vm-1子集内元素个数相同的集合，把指针i指向满足条件的集合 i -; i +;/把指针i指向满足条件的第一个集合 while（i m） int ch;/用于表示下一个即将加入子集的数字 if（

14、i = 0）ch = 0;/如果i指向V中第一个元素 else String chStr = + Vi.charAt（Vi.length（） - 1）;/找出Vi中最后一个数字 ch = Integer.parseInt（chStr）;/转换成整型 /比ch大而又比N-1（因为这里顶点是从0开始的）小的数字应该加在子集中 while（ch N - 1） Vm + = Vi + （+ ch）; i +;/对已存在的自己逐个扫描添加 initialV（m,Vm - 1.length（）;/递归调用 /判断自己Vj中是否存在指定元素，即行号i boolean exclude（int i , int

15、j） String str = + i;/把i转换成字符串 if（Vj.contains（str） /System.out.println（i + i return false;/如若存在，则返回false else return true; /获得子集Vj中除指定元素k外的元素，用字符串形式表示 public String getSubString（int k , int j） if（Vj.length（） = 1）return /如果子集中只有一个元素，则返回空串 if（k = 0） return Vj.substring（1,Vj.length（）;/如果k是第一个元素，则返回其后面的元

16、素 else if（k = Vj.length（） - 1） return Vj.substring（0, Vj.length（）-1）;/如果k是最后一个元素，则返回其前面的元素 else return （Vj.substring（0, k） + Vj.substring（k+1, Vj.length（）;/返回除k外的元素 /找出V中与str相同元素的下标号,即找出上一个子集 public int stringEqual（String str） /if（str.equals（）return 0; int i = 0; V.length） if（Vi.equals（str） return i

17、; return -1;/如若没找到，则返回错误符号-1 /求最小距离 public int min（int i , int j） int k = 0;/用于记录Vj中元素个数 String vStr = + Vj.charAt（k）;/铭记Vj.charAt（k）得到的是字符型，转换成整形后是字母对应的ASC码！ int v = Integer.parseInt（vStr）;/把位置k处的字符转换成整形 String str = getSubString（k,j）;/获得Vj中除位置k处外的字符串 /System.out.println（min + str + stringEqual（str

18、） + v）; if（stringEqual（str） = -1）System.exit（0）; int min = civ + dvstringEqual（str）;/先把最小的距离赋值给从Vj中第一个顶点出发的距离 /System.out.println（min）; /stringEqual（str）表示返回与上面获得的字符串相同的V中元素的下标，即找上一个子集 pathij = pathvstringEqual（str） + i; k +; /寻找最小距离 while（k civ + dvstringEqual（str） min = civ + dvstringEqual（str）; p

19、athij = pathvstringEqual（str） + i; k +; /Vj.substring（beginIndex, endIndex） /System.out.println（pathij）; return min;/返回最小值 /处理函数 public void dynamic（） d = new intN（int）Math.pow（2,N-1）;/分配空间 path = new StringN（int）Math.pow（2,N-1）; for（int i = 1 ; i +）/初始化第一列 di0 = ci0; pathi0 = 0/初始化第一个元素，即为出发城市顶点 /S

20、ystem.out.print（di0 + /初始化后面的元素 int j = 1; for（ ; （int）Math.pow（2,N-1） -1 ; j +） for（int i = 1 ; i +） if（exclude（i,j）/判断V子集中是否包含当前顶点，即Vj中是否包含i /System.out.println（done！ + i + + j）; dij = min（i , j）;/寻找最小距离 d0j = min（0,j）;/初始化组后一列 /输出中间结果，各个数组，用于调试程序 public void print（） （int）Math.pow（2,N-1） ; System.

21、out.print（Vi + /for（int i = 0 ; c.length ; ） System.out.println（）; System.out.print（cij + System.out.println（）; System.out.print（dij + /输出最短路径 public void printShortestPath（） /输出所经城市 System.out.print（经过城市： String str = path0（int）Math.pow（2,N-1） - 1; /System.out.println（str）; System.out.print（str.cha

22、rAt（str.length（） - 1）; for（int i = str.length（） - 2 ; i = 0 ; i -） System.out.print（- + str.charAt（i）;会有最短路径最短路径为： + d0（int）Math.pow（2,N-1） - 1）; /主函数 public static void main（String args） TSPDynamic TSP = new TSPDynamic（）; TSP.dynamic（）;/求最短路径 /TSP.print（）; TSP.printShortestPath（）;/输出最短路径/测试数据/*99999 3 6 75 99999 2 36 4 99999 23 7 5 99999 */6.1.2结果（1）（2）（3）（4）6.2贪心法6.2.1源代码public class TSPGreedNode int V;/存放旅行所经过的城市顶点/存放每两座城市间的距离，注意：若路径不存在或同一城市间距离为无穷大 int path;/存放旅行所经过的每两