分支定界问题1

分支定界(branch and bound) s搜索法是一种在问题的解空间树上搜索问题的解的方法。

但与回溯算法不同，分支定界算法采用广度优先或最小耗费优先的方法搜索解空间树，并且，在分支定界算法中，每一个活结点只有一次机会成为扩展结点。

利用分支定界算法对问题的解空间树进行搜索，它的搜索策略是： 1 ．产生当前扩展结点的所有孩子结点；
2 ．在产生的孩子结点中，抛弃那些不可能产生可行解（或最优解）的结点；
3 ．将其余的孩子结点加入活结点表；
4 ．从活结点表中选择下一个活结点作为新的扩展结点。

如此循环，直到找到问题的可行解（最优解）或活结点表为空。

从活结点表中选择下一个活结点作为新的扩展结点，根据选择方式的不同，分支定界算法通常可以分为两种形式： 1 ．FIFO(First In First Out) 分支定界算法：按照先进先出原则选择下一个活结点作为扩展结点，即从活结点表中取出结点的顺序与加入结点的顺序相同。

2 ．最小耗费或最大收益分支定界算法：在这种情况下，每个结点都有一个耗费或收益。

如果要查找一个具有最小耗费的解，那么要选择的下一个扩展结点就是活结点表中具有最小耗费的活结点；如果要查找一个具有最大收益的解，那么要选择的下一个扩展结点就是活结点表中具有最大收益的活结点。

又称分支定界搜索法。

过程系统综合的一类方法。

该法是将原始问题分解，产生一组子问题。

分支是将一组解分为几组子解，定界是建立这些子组解的目标函数的边界。

如果某一子组的解在这些边界之外，就将这一子组舍弃（剪枝）。

分支定界法原为运筹学中求解整数规划(或混合整数规划)问题的一种方法。

用该法寻求整数最优解的效率很高。

将该法原理用于过程系统综合可大大减少需要计算的方案数日。

分支定界法的思想是：首先确定目标值的上下界，边搜索边减掉搜索树的某些支，提高搜索效率。

在竞赛中，我们有时会碰到一些题目，它们既不能通过建立数学模型解决，又没有现成算法可以套用，或者非遍历所有状况才可以得出正确结果。

这时，我们就必须采用搜索算法来解决问题。

搜索算法按搜索的方式分有两类，一类是深度优先搜索，一类是广度优先搜索。

我们知道，深度搜索编程简单，程序简洁易懂，空间需求也比较低，但是这种方法的时间复杂度往往是指数级的，倘若不加优化，其时间效率简直无法忍受；而广度优先搜索虽然时间复杂度比前者低一些，但其庞大的空间需求量又往往让人望而却步。

所以，对程序进行优化，就成为搜索算法编程中最关键的一环。

本文所要讨论的便是搜索算法中优化程序的一种基本方法枣“剪枝”。

什么是剪枝相信刚开始接触搜索算法的人，都做过类似迷宫这样的题目吧。

我们在“走迷宫”的时候，一般回溯法思路是这样的：
1、这个方向有路可走，我没走过
2、往这个方向前进
3、是死胡同，往回走，回到上一个路口
4、重复第一步，直到找着出口这样的思路很好理解，编程起来也比较容易。

但是当迷宫的规模很大时，回溯法的缺点便暴露无遗：搜索耗时极巨，无法忍受。

我们可不可以在向某个方向前进时，先一步判断出这样走会不会走到死胡同里呢？这样一来，搜索的时间不就可以减少了吗？答案是：可以的。

剪枝的概念，其实就跟走迷宫避开死胡同差不多。

若我们把搜索的过程看成是对一棵树的遍历，那么剪枝顾名思义，就是将树中的一些“死胡同”，不能到达我们需要的解的枝条“剪”掉，以减少搜索的时间。

搜索算法，绝大部分需要用到剪枝。

然而，不是所有的枝条都可以剪掉，这就需要通过设计出合理的判断方法，以决定某一分支的取舍。

在设计判断方法的时候，需要遵循一定的原则。

剪枝的原则1、正确性正如上文所述，枝条不是爱剪就能剪的。

如果随便剪枝，把带有最优解的那一分支也剪掉了的话，剪枝也就失去了意义。

所以，剪枝的前提是一定要保证不丢失正确的结果。

2、准确性在保证了正确性的基础上，我们应该根据具体问题具体分析，采用合适的判断手段，使不包含最优解的枝条尽可能多的被剪去，以达到程序“最优化”的目的。

可以说，剪枝的准确性，是衡量一个优化算法好坏的标准。

3、高效性设计优化程序的根本目的，是要减少搜索的次数，使程序运行的时间减少。

但为了使搜索次数尽可能的减少，我们又必须花工夫设计出一个准确性较高的优化算法，而当算法的准确性升高，其判断的次数必定增多，从而又导致耗时的增多，这便引出了矛盾。

因此，如何在优化与效率之间寻找一个平衡点，使得程序的时间复杂度尽可能降低，同样是非常重要的。

倘若一个剪枝的判断效果非常好，但是它却需要耗费大量的时间来判断、比较，结果整个程序运行起来也跟没有优化过的没什么区别，这样就太得不偿失了。

综上所述，我们可以把剪枝优化的主要原则归结为六个字：正确、准确、高效。

剪枝算法按照其判断思路可大致分成两类：可行性剪枝及最优性剪枝。

对于分支定界算法,上界是已求得的可行解的目标函数值中的最小者,分为初始上界和在探测过程中产生的动态上界.
分支定界法在求最优解的迭代过程中, 若某结点估计的下界不小于已知的上界, 则不必从
该节点往下继续搜索. 因此若能产生一个较好的上界, 可以消除许多不必要的列举计算.
分支定界算法的实现在描述分支定界算法步骤之前, 先对算法涉及到的有关术语进行定义如下: p ——分支层数; C*——当前最优目标函数值; P*——相应于
C*的工件顺序; P1——当前节点(现在需要进行分支的节点)所对应的部分序列.
分支定界算法的实施步骤如下: 步骤1初始化: 设置p = 0, P 1 = Á (空集) , C* = ∞.设当前节点总是与P 1 相对应. 此时, 当前节点即根节点. 步骤2计算从当前节点分支得到的各个子节点的下界, 并按下界值由小到大对各子节点排序. 令p ←p + 1.
步骤3如果当前节点被探测尽, 令p ←p - 1, 转步骤6. 否则, 设当前层(第p 层) 各活动子节点中具有最小下界值的节点为Q , 则在P 1末尾加入Q 对应第p 位置上的工件, 此时的当前节点转为Q , 转步骤4. 步骤4因为当前节点是同层同父节点具有最小下界值的节点, 如果当前节点下界值大于或等于C* , 则不必再搜索当前节点及其同层同父的活动节点, 这样, 当前节点的上一层节点(父节点)被探测尽, p ←p - 1, 去掉P 1 中的最后一个工件,转步骤6. 否则, 转步骤5. 步骤5如果p = n, 则得到一个较优顺序.令P* = P 1, C* 是当前节点的下界值, p ←p - 1,去掉P 1 中最后一个工件, 转步骤6; 否则转步骤2. 步骤6若p ≠ 0, 去掉P 1 中最后一个工件,转步骤3; 否则, 算法停止. C* 是最优的目标函数值, P* 是最优顺序.
// 分支限界法之旅行商问题.cpp : Defines the entry point for the console application.//#include "stdafx.h"#include"stdlib.h"#include <iostream.h>//各个城市之间的路径长度int CityVal[4][4]={ {0 ,30,6,4}, {30,0 ,5,10}, {6,5,0 ,20 }, {4,10,20 ,0}, };struct CityNum{ int MinVal[2]; //最短路径长度int PastCity[2][4]; //走过的路径}city[4];int find(int j,int k,int boolnum){ int i; for(i=0;i<4;i++) if(j==city[k].PastCity[boolnum][i]) return 1; return 0;}void copy(int k,int boolnum1,int j,int boolnum2){ int i=0; while(city[k].PastCity[boolnum1][i]!=0)
{ city[j].PastCity[boolnum2][i]=city[k].PastCity[boolnum1][i]; i++; } city[j].PastCity[boolnum2][i]=k; }void main(){ int i,j,k; int boolnum=0; int TempMinVal; int TempCity; int TempVal; for(i=0;i<2;i++) for(j=0;j<4;j++) for(k=0;k<4;k++) city[k].PastCity[i][j]=0; for(i=1;i<4;i++) city[i].MinVal[boolnum]=CityVal[i][0]; for(i=0;i<2;i++) { for(j=1;j<4;j++) { TempMinVal=32767; for(k=1;k<4;k++) { if(j!=k&&!find(j,k,boolnum)) { TempVal=CityVal[j][k]+city[k].MinVal[boolnum]; if(TempMinVal>TempVal) { TempMinVal=TempVal; TempCity=k; } } } city[j].MinVal[(boolnum+1)%2]=TempMinVal;
copy(TempCity,boolnum,j,(boolnum+1)%2); } boolnum=(boolnum+1)%2; } TempMinVal=32767; for(i=1;i<4;i++) { TempVal=CityVal[0][i]+city[i].MinVal[boolnum]; if(TempMinVal>TempVal) { TempMinVal=TempVal; TempCity=i; } } cout<<"最小费用为:"<<TempMinVal<<endl;cout<<"最佳周游路径为:"; cout<<1<<' '; cout<<TempCity+1<<' '; for(i=1;i>=0;i--) cout<<city[TempCity].PastCity[boolnum][i]+1<<' '; cout<<1<<endl;}。