北航研究生 算法设计与分析 Assignment_2.docVIP

用分支定界算法求以下问题： 某公司于乙城市的销售点急需一批成品，该公司成品生产基地在甲城市。甲城市与乙城市之间共有 n 座城市，互相以公路连通。甲城市、乙城市以及其它各城市之间的公路连通情况及每段公路的长度由矩阵M1 给出。每段公路均由地方政府收取不同额度的养路费等费用，具体数额由矩阵M2 给出。 请给出在需付养路费总额不超过 1500 的情况下，该公司货车运送其产品从甲城市到乙城市的最短运送路线。 具体数据参见文件： m1.txt: 各城市之间的公路连通情况及每段公路的长度矩阵(有向图); 甲城市为城市Num.1，乙城市为城市Num.50。 m2.txt: 每段公路收取的费用矩阵（非对称）。 思想：利用Floyd算法的基本方法求解。 程序实现流程说明： 1.将m1.txt和m2.txt的数据读入两个50×50的数组。 2.用Floyd算法求出所有点对之间的最短路径长度和最小费用。 3.建立一个堆栈，初始化该堆栈。 4.取出栈顶的结点，检查它的相邻的所有结点，确定下一个当前最优路径上的结点，被扩展的结点依次加入堆栈中。在检查的过程中，如果发现超出最短路径长度或者最小费用，则进行”剪枝”，然后回溯。 5.找到一个解后，保存改解，然后重复步骤4。 6.重复步骤4、5，直到堆栈为空，当前保存的解即为最优解。 时间复杂度分析： Floyd算法的时间复杂度为，N为所有城市的个数。 该算法的时间复杂度等于DFS的时间复杂度，即O(N+E)。其中，E为所有城市构成的有向连通图的边的总数。但是因为采用了剪枝，会使实际运行情况的比较次数远小于E。 求解结果： 算法所得结果: 甲乙之间最短路线长度是:464 最短路线收取的费用是:1448 最短路径是: 1 3 8 11 15 21 23 26 32 37 39 45 47 50 C源代码（注意把m1.txt与m2.txt放到与源代码相同的目录下， 下面代码可直接复制运行）： #includestdlib.h #includestdio.h #includetime.h #includestring.h #define N 50 #define MAX 52 void input(int a[N][N],int b[N][N]); void Floyd(int d[N][N]); void fenzhi(int m1[N][N],int m2[N][N],int mindist[N][N],int mincost[N][N]); int visited[N],bestPath[N]; void main() { clock_t start,finish; double duration; int i,j,mindist[N][N],mincost[N][N],m1[N][N],m2[N][N]; /* m1[N][N]和m2[N][N]分别代表题目所给的距离矩阵和代价矩阵 */ // int visited[N],bestPath[N]; FILE *fp,*fw; // system(cls); time_t ttime; time(ttime); printf(%s,ctime(ttime)); start=clock(); for(i=0;iN;i++) { visited[i]=0; bestPath[i]=0; } fp=fopen(m1.txt,r); /* 把文件中的距离矩阵m1读入数组mindist[N][N] */ if(fp==NULL) { printf(can not open file\n); return; } for(i=0;iN;i++) for(j=0;jN;j++) fscanf(fp,%d,mindist[i][j]); fclose(fp); /* 距离矩阵m1读入完毕 */ fp=fopen(m2.txt,r); /* 把文件中的代价矩阵m2读入数组mincost[N][N] */ if(fp==NULL) { printf(can not open file\n); return; } for(i=0;iN;i++) for(j=0;jN;j++) fscanf(fp,%d,mincost[i][j]); fclose(fp); /* 代