实验二动态规划算法的应用.docVIP

实验二：动态规划算法的应用 ——旅行商问题 院系：信息学院05计算机系 班级： 1班 姓名：阙寿辉 学号：22120051203884 一、问题描述： 某售货员要到若干城市去推销商品，已知各城市之间的路程（或旅费）。他要选定一条从驻地出发，经过每个城市一遍，最后回到驻地的路线，使总的路程（或旅费）最小。 例如：给定4个城市{1，2，3，4}及其各城市之间的路程 请设计出一个动态规划算法，并用程序实现如下功能： 1.任给一个输入实例，能输出最短路程及其路线 2.能用图形演示旅行商的推销路线 二、算法思想： 在无向完全图中，对于任意两个顶点vi和vj，我们可以在多项式时间内找到vi和vj这两个顶点之间的所有路径，选择其中路程最短的一条，令S[i,j]表示vi和vj这两个顶点之间最短距离的那条路径。搜索路径S[i,j]，找到vi到达的在S[i,j]上的第一个顶点，记该顶点为vk，将其记录在数组中R[][]，递归查找vi到vk和vk到vj的最短路径及其相应权值，最后将数组D[]中的顶点和权值之和打印出来即为所求，并用画图函数将行经过程画出。 三、具体实现： 1. int Exp(int k)//计算2的k次方:二进制数左移一位实现*2 { k=1k; return k; } 该函数用以计算2的k次方(二进制数左移一位实现*2计算)。设计思想：用一个二进制数来表示顶点集，其各位数值表示的含义为：当该位为i=1时表示顶点i+1在集合s中，否则若i=0顶点i+1不在集合s中。例如k = 1表示顶点2在集合s中。对给定n个顶点，共需要用n-1位二进制位来记录它的全部顶点。 2. int Readfile()//读入测试文件 { int space,i,j; char s[10]; printf(请输入数据文件名:); gets(s); if(freopen(s,r,stdin)==NULL) { printf(数据文件读取失败!\n); return 0; } scanf(%d,n); printf(经过城市的个数为:n=%d\n,n); G=(int **)malloc(sizeof(int *)*(n+1));//开辟空间，用以存储图的邻接矩阵 for(i=0;i=n;i++) G[i]=(int *)malloc(sizeof(int)*(n+1));//为邻接矩阵分配空间 space=Exp(n-1); D=(int **)malloc(sizeof(int *)*(n+1)); for(i=0;i=n;i++) D[i]=(int *)malloc(sizeof(int)*(space)); R=(int **)malloc(sizeof(int *)*(n+1)); for(i=0;i=n;i++) R[i]=(int *)malloc(sizeof(int)*(space)); for(i=1;i=n;i++) for(j=1;j=n;j++) scanf(%d,G[i][j]); return 1; } 该函数从当前目录下读入测试文件。该函数还可以实现对用户输入的文件名的判断。若输入的文件名在当前目录下存在，那么程序继续运行，否则直接退出以节省程序运行时间。另外，在函数体内实现了对邻接矩阵的空间分配。 3. void FillTable()//采用自底向上填表 { int i,k; unsigned int max,s2; max=Exp(n-1)-1; for(s=1;s=max;s++) { for(i=1;i=n;i++) if((s1(i-2))==0)//如果第i个顶点不属于集合s { D[i][s]=MAXNUM; for(k=2;k=n;k++) { if((s1(k-2))!=0) { s2=s~(1(k-2)); if(G[i][k]+D[k][s2]D[i][s]) { D[i][s]=G[i][k]+D[k][s2]; R[i][s]=k; } } } } } } 在程序运行过程中用该函数来记录添加到集合s中的顶点，本函数采用自底向上来填写。 4. void Tsp()//动态规划求解 { int i; s=s0; for(i=2;i=n;i++) { D[i][s]=G[i][1]; R[i][s]=1; } FillTable(); } 函数void Tsp()实现采用动态规