(南开大学ACM资料)第6章-分支限界法

1第六章分支限界法2第六章分支限界法本章主要知识点•6.1分支限界法的基本思想•6.2单源最短路径问题•6.3装载问题•6.4布线问题•6.50－1背包问题•6.6最大团问题•6.7旅行售货员问题•6.8电路板排列问题•6.9批处理作业调度36.1分支限界法的基本思想1.分支限界法与回溯法的不同（1）求解目标：回溯法的求解目标是找出解空间树中满足约束条件的所有解，而分支限界法的求解目标则是找出满足约束条件的一个解，或是在满足约束条件的解中找出在某种意义下的最优解。（2）搜索方式的不同：回溯法以深度优先的方式搜索解空间树，而分支限界法则以广度优先或以最小耗费优先的方式搜索解空间树。46.1分支限界法的基本思想2.分支限界法基本思想分支限界法常以广度优先或以最小耗费（最大效益）优先的方式搜索问题的解空间树。在分支限界法中，每一个活结点只有一次机会成为扩展结点。活结点一旦成为扩展结点，就一次性产生其所有儿子结点。在这些儿子结点中，导致不可行解或导致非最优解的儿子结点被舍弃，其余儿子结点被加入活结点表中。此后，从活结点表中取下一结点成为当前扩展结点，并重复上述结点扩展过程。这个过程一直持续到找到所需的解或活结点表为空时为止。56.1分支限界法的基本思想3.常见的两种分支限界法（1）队列式(FIFO)分支限界法按照队列先进先出（FIFO）原则选取下一个节点为扩展节点。（2）优先队列式分支限界法按照优先队列中规定的优先级选取优先级最高的节点成为当前扩展节点。66.2单源最短路径问题1.问题描述下面以一个例子来说明单源最短路径问题：在下图所给的有向图G中，每一边都有一个非负边权。要求图G的从源顶点s到目标顶点t之间的最短路径。76.2单源最短路径问题下图是用优先队列式分支限界法解有向图G的单源最短路径问题产生的解空间树。其中，每一个结点旁边的数字表示该结点所对应的当前路长。86.2单源最短路径问题2.算法思想解单源最短路径问题的优先队列式分支限界法用一极小堆来存储活结点表。其优先级是结点所对应的当前路长。算法从图G的源顶点s和空优先队列开始。结点s被扩展后，它的儿子结点被依次插入堆中。此后，算法从堆中取出具有最小当前路长的结点作为当前扩展结点，并依次检查与当前扩展结点相邻的所有顶点。如果从当前扩展结点i到顶点j有边可达，且从源出发，途经顶点i再到顶点j的所相应的路径的长度小于当前最优路径长度，则将该顶点作为活结点插入到活结点优先队列中。这个结点的扩展过程一直继续到活结点优先队列为空时为止。96.2单源最短路径问题3.剪枝策略在算法扩展结点的过程中，一旦发现一个结点的下界不小于当前找到的最短路长，则算法剪去以该结点为根的子树。在算法中，利用结点间的控制关系进行剪枝。从源顶点s出发，2条不同路径到达图G的同一顶点。由于两条路径的路长不同，因此可以将路长长的路径所对应的树中的结点为根的子树剪去。106.2单源最短路径问题while(true){//搜索问题的解空间for(intj=1;j=n;j++)if(a[enode.i][j]Float.MAX_VALUE&&enode.length+a[enode.i][j]dist[j]){//顶点i到顶点j可达，且满足控制约束dist[j]=enode.length+a[enode.i][j];p[j]=enode.i;HeapNodenode=newHeapNode(j,dist[j]);heap.put(node);//加入活结点优先队列}if(heap.isEmpty())break;elseenode=(HeapNode)heap.removeMin();}顶点I和j间有边，且此路径长小于原先从原点到j的路径长116.3装载问题1.问题描述有一批共个集装箱要装上2艘载重量分别为C1和C2的轮船，其中集装箱i的重量为Wi，且211ccwnii装载问题要求确定是否有一个合理的装载方案可将这个集装箱装上这2艘轮船。如果有，找出一种装载方案。容易证明：如果一个给定装载问题有解，则采用下面的策略可得到最优装载方案。(1)首先将第一艘轮船尽可能装满；(2)将剩余的集装箱装上第二艘轮船。126.3装载问题2.队列式分支限界法在算法的while循环中，首先检测当前扩展结点的左儿子结点是否为可行结点。如果是则将其加入到活结点队列中。然后将其右儿子结点加入到活结点队列中(右儿子结点一定是可行结点)。2个儿子结点都产生后，当前扩展结点被舍弃。活结点队列中的队首元素被取出作为当前扩展结点，由于队列中每一层结点之后都有一个尾部标记-1，故在取队首元素时，活结点队列一定不空。当取出的元素是-1时，再判断当前队列是否为空。如果队列非空，则将尾部标记-1加入活结点队列，算法开始处理下一层的活结点。136.3装载问题2.队列式分支限界法while(true){if(ew+w[i]=c)enQueue(ew+w[i],i);//检查左儿子结点enQueue(ew,i);//右儿子结点总是可行的ew=((Integer)queue.remove()).intValue();//取下一扩展结点if(ew==-1){if(queue.isEmpty())returnbestw;queue.put(newInteger(-1));//同层结点尾部标志ew=((Integer)queue.remove()).intValue();//取下一扩展结点i++;//进入下一层}}146.3装载问题3.算法的改进节点的左子树表示将此集装箱装上船，右子树表示不将此集装箱装上船。设bestw是当前最优解；ew是当前扩展结点所相应的重量；r是剩余集装箱的重量。则当ew+rbestw时，可将其右子树剪去，因为此时若要船装最多集装箱，就应该把此箱装上船。另外，为了确保右子树成功剪枝，应该在算法每一次进入左子树的时候更新bestw的值。156.3装载问题3.算法的改进//检查左儿子结点intwt=ew+w[i];if(wt=c){//可行结点if(wtbestw)bestw=wt;//加入活结点队列if(in)queue.put(newInteger(wt));}提前更新bestw//检查右儿子结点if(ew+rbestw&&in)//可能含最优解queue.put(newInteger(ew));ew=((Integer)queue.remove()).intValue();//取下一扩展结点右儿子剪枝166.3装载问题4.构造最优解为了在算法结束后能方便地构造出与最优值相应的最优解，算法必须存储相应子集树中从活结点到根结点的路径。为此目的，可在每个结点处设置指向其父结点的指针，并设置左、右儿子标志。privatestaticclassQNode{QNodeparent;//父结点booleanleftChild;//左儿子标志intweight;//结点所相应的载重量176.3装载问题找到最优值后，可以根据parent回溯到根节点，找到最优解。//构造当前最优解for(intj=n;j0;j--){bestx[j]=(e.leftChild)?1:0;e=e.parent;}186.3装载问题5.优先队列式分支限界法解装载问题的优先队列式分支限界法用最大优先队列存储活结点表。活结点x在优先队列中的优先级定义为从根结点到结点x的路径所相应的载重量再加上剩余集装箱的重量之和。优先队列中优先级最大的活结点成为下一个扩展结点。以结点x为根的子树中所有结点相应的路径的载重量不超过它的优先级。子集树中叶结点所相应的载重量与其优先级相同。在优先队列式分支限界法中，一旦有一个叶结点成为当前扩展结点，则可以断言该叶结点所相应的解即为最优解。此时可终止算法。196.4布线问题算法的思想解此问题的队列式分支限界法从起始位置a开始将它作为第一个扩展结点。与该扩展结点相邻并且可达的方格成为可行结点被加入到活结点队列中，并且将这些方格标记为1，即从起始方格a到这些方格的距离为1。接着，算法从活结点队列中取出队首结点作为下一个扩展结点，并将与当前扩展结点相邻且未标记过的方格标记为2，并存入活结点队列。这个过程一直继续到算法搜索到目标方格b或活结点队列为空时为止。即加入剪枝的广度优先搜索。206.4布线问题Position[]offset=newPosition[4];offset[0]=newPosition(0,1);//右offset[1]=newPosition(1,0);//下offset[2]=newPosition(0,-1);//左offset[3]=newPosition(-1,0);//上定义移动方向的相对位移for(inti=0;i=size+1;i++){grid[0][i]=grid[size+1][i]=1;//顶部和底部grid[i][0]=grid[i][size+1]=1;//左翼和右翼}设置边界的围墙216.4布线问题for(inti=0;inumOfNbrs;i++){nbr.row=here.row+offset[i].row;nbr.col=here.col+offset[i].col;if(grid[nbr.row][nbr.col]==0){//该方格未标记grid[nbr.row][nbr.col]=grid[here.row][here.col]+1;if((nbr.row==finish.row)&&(nbr.col==finish.col))break;q.put(newPosition(nbr.row,nbr.col));}}找到目标位置后，可以通过回溯方法找到这条最短路径。226.50-1背包问题•算法的思想首先，要对输入数据进行预处理，将各物品依其单位重量价值从大到小进行排列。在下面描述的优先队列分支限界法中，节点的优先级由已装袋的物品价值加上剩下的最大单位重量价值的物品装满剩余容量的价值和。算法首先检查当前扩展结点的左儿子结点的可行性。如果该左儿子结点是可行结点，则将它加入到子集树和活结点优先队列中。当前扩展结点的右儿子结点一定是可行结点，仅当右儿子结点满足上界约束时才将它加入子集树和活结点优先队列。当扩展到叶节点时为问题的最优值。236.50-1背包问题上界函数while(i=n&&w[i]=cleft)//n表示物品总数，cleft为剩余空间{cleft-=w[i];//w[i]表示i所占空间b+=p[i];//p[i]表示i的价值i++;}if(i=n)b+=p[i]/w[i]*cleft;//装填剩余容量装满背包returnb;//b为上界函数246.50-1背包问题while(i!=n+1){//非叶结点doublewt=cw+w[i];if(wt=c){//左儿子结点为可行结点if(cp+p[i]bestp)bestp=cp+p[i];addLiveNode(up,cp+p[i],cw+w[i],i+1,enode,true);}up=bound(i+1);if(up=bestp)//检查右儿子节点addLiveNode(up,cp,cw,i+1,enode,false);//取下一个扩展节点（略）}分支限界搜索过程256.6最大团问题1.问题描述给定无向图G=(V，E)。如果UV，且对任意u，vU有(u，v)E，则称U是G的完全子图。G的完全子图U是G的团当且仅当U不包含在G的更大的完全子图中。G的最大团是指G中所含顶点数最多的团。下图G中，子集{1，2}是G的大小为2的完全子图。这个完全子图不是团，因为它被G的更大的完全子图{1，2，5}包含。{1，2，5}是G的最大团。{1，4，5}和{2，3，5}也是G的最大团。266.6最大团问题2.上界函数用变量cliqueSize表示与该结点相应的团的顶点数；level表示结点在子集空间树中所处的层次；用cliqueSize+n-level+1作为顶点数上界upperSize的值。在此优先队列式分