设计 郑宗汉郑明 第8章 分支与限界.pptVIP

if (ynode-k==2) { //矩阵只剩2阶 if (ynode-c[0][0]==0) {//登记最后两边 ynode-ad[ynode-row_init[0]]= ynode-col_init[0]; ynode-ad[ynode-row_init[1]]= ynode-col_init[1]; } else { ynode-ad[ynode-row_init[0]]= ynode-col_init[1]; ynode-ad[ynode-row_init[1]]= ynode-col_init[0]; } ynode-k = 0; } insert(heap,n_heap,y); //y插入堆中 delete xnode; //释放没用的x缓冲区 x = delete_min(heap,n_heap); xnode = x.p; } w = xnode-w //保存最短路费用 for (i=0;in;i++) //保存顶点邻接表 ad[i] = xnode-ad[i]; delete xnode; //释放x结点缓冲区 for (i=1;i=n_heap;i++) //释放堆缓 delete heap[i].p; delete heap; return w; } 算法的复杂性分析: 算法的时间复杂度: O(n22n) 算法所需存储空间: 　　每个结点需O(n2)存放费用矩阵, 　　共2n个结点, 算法的空间复杂度: O(n22n) * 令S是归约矩阵中cij=0的元素集合, 　dkl是S中使dij最大的元素, 即 　　　dkl=maxS{dij} 则边(vk,vl)为所选择的分支方向。 注: dij为第i行除cij外的最小元素与第j列除cij外的最小元素之和。 考虑上述5顶点TSP问题, 其归约矩阵如右所示: ? 0 16 3 3 0 ? 13 22 21 19 15 ? 2 0 4 45 19 ? 0 16 2 0 0 ? 2 1 0 4 3 4 3 2 1 0 dij最大的元素为d10=17, 由此确定所选择的方向为v1v0, 建立分支结点Y和Y?: w(Y?)= w(X)+dkl= 48+17= 65 c01=c10=c24=c34=c42 　 =c43=0 d01=3+2=5, d10=13+4=17 d24=2+0=2, d34=4+0=4 d42=0+13=13, d43=0+2=2 使用分支限界法的求解过程中, 将动态生成很多结点, 用结点表存放动态生成的结点信息。 由于必须按费用的下界来确定搜索方向, 故可用优先队列或堆来维护结点表。在此使用最小堆维护结点表。 8.2.3 TSP问题的求解过程 分支限界法求解TSP的求解过程: 分配堆缓冲区, 初始化为空堆; 建立父结点X。X的费用矩阵X.c初始化为费用矩阵c, X的费用矩阵的阶X.k初始化为n; 归约X.c得归约常数h, 令结点X的下界X.w=h; 初始化回路的顶点邻接表X.ad; 由归约后的X.c中所有cij=0的元素计算dij: (4)选取dij最大的元素dkl, 边vkvl为分支方向; (5)建立子结点Y?。把X的费用矩阵X.c复制到Y?.c; 把Y?.c中元素ckl置为?, 得到Y?.c; 按w(Y?)=w(X)+dkl计算Y?的下界Y?.w; 把X的回路顶点邻接表X.ad复制到Y?.ad, X.k复制到Y?.k; 把结点Y?按Y?.w插入最小堆; (6)建立子结点Y。把X的费用矩阵X.c复制到Y.c, 把X的回路顶点邻接表X.ad复制到Y.ad, X.k复制到Y.k; 按回路顶点邻接表Y.ad把费用矩阵Y.c的有关元素置为?; (7)删去Y.c的第k行第l列元素, 使Y.k减1, 即费用矩阵Y.c的阶数减1; 归约降阶后的费用矩阵Y.c, 按w(Y)=w(X)+h计算结点Y的下界Y.w; (8)若Y.k=2, 直接判断最矩回路的两条边, 并登记到邻接表Y.ad, 使Y.k=0; (9)把结点Y按Y.w插入最小堆; (10)取堆顶元素作为结点X, 若X.k=0且X.w为结点表中的最大(小)值, 则算法结束, 否则转(3)计算d