【算法分析与设计】分支限界法(上)


一、学习要点

  理解分支限界法的剪枝搜索策略。
  掌握分支限界法的算法框架
  (1)队列式(FIFO)分支限界法
  (2)优先队列式分支限界法
  通过应用范例学习 分支限界法的设计策略
  (1)单源最短路径问题
  (2)装载问题;
  (3)布线问题
  (4)0-1背包问题;
  (5)最大团问题;
  (6)旅行售货员问题
  (7)电路板排列问题
  (8)批处理作业调度问题

  三个人和三个鬼要过河,只有一条船能渡,船只有两个位置,人和鬼都会划船。两岸任何一边,如果鬼的个数超过人的个数,鬼就会把人吃掉。问如何在保证所有人安全的情况下,将三人三鬼全部渡过河去
  隐含条件:船不会自己单独行驶,必须有人或鬼驾船才能动。试着用穷举法求解
  回溯法也可求,但是效率低,有无其他解法?


二、分支限界法的基本思想

  分支限界法与回溯法
  (1)求解目标回溯法的求解目标是找出解空间树中满足约束条件的所有解,而分支限界法的求解目标则是找出满足约束条件的一个解,或是在满足约束条件的解中找出在某种意义下的最优解
  (2)搜索方式的不同回溯法以深度优先的方式搜索解空间树,而分支限界法则以广度优先或以最小耗费优先的方式搜索解空间树

  分支限界法常以广度优先或以最小耗费(最大效益)优先的方式搜索问题的解空间树
  在分支限界法中,每一个活结点只有一次机会成为扩展结点。活结点一旦成为扩展结点,就一次性产生其所有儿子结点。在这些儿子结点中,导致不可行解或导致非最优解的儿子结点被舍弃,其余儿子结点被加入活结点表中。
  此后,从活结点表中取下一结点成为当前扩展结点,并重复上述结点扩展过程。这个过程一直持续到找到所需的解或活结点表为空时为止。

  常见的两种分支限界法
  (1)队列式(FIFO)分支限界法
  按照 队列先进先出(FIFO)原则 选取下一个节点为扩展节点。
  (2)优先队列式分支限界法
  按照 优先队列中规定的优先级选取优先级最高的节点 成为当前扩展节点。

  0-1背包问题的队列式分支限界法:
  W=[16,15,15]
  P=[45,25,25]
  C=30
  极大堆
在这里插入图片描述
  旅行售货员问题的优先队列
  可以 用一个限界函数在搜索过程中裁减子树,以减少产生的活结点
  极小堆
在这里插入图片描述
在这里插入图片描述


三、单源最短路径问题

3.1 问题描述

  下面以一个例子来说明单源最短路径问题:在下图所给的有向图G中,每一边都有一个非负边权。要求图G的从源顶点s到目标顶点t之间的最短路径。
在这里插入图片描述
  下图是用优先队列式分支限界法解有向图G的单源最短路径问题产生的解空间树。其中,每一个结点旁边的数字表示该结点所对应的当前路长
在这里插入图片描述
在这里插入图片描述


3.2 算法思想

  解 单源最短路径问题的优先队列式分支限界法 用一极小堆来存储活结点表。其优先级是结点所对应的当前路长
  算法从图G的源顶点s和空优先队列开始。结点s被扩展后,它的儿子结点被依次插入堆中。此后,算法从堆中取出具有最小当前路长的结点作为当前扩展结点,并依次检查与当前扩展结点相邻的所有顶点。如果从当前扩展结点i到顶点j有边可达,且从源出发,途经顶点i再到顶点j的所相应的路径的长度小于当前最优路径长度,则将该顶点作为活结点插入到活结点优先队列中。这个结点的扩展过程一直继续到活结点优先队列为空时为止。


3.3 剪枝策略

  在算法扩展结点的过程中,一旦发现一个结点的下界不小于当前找到的最短路长,则算法剪去以该结点为根的子树

  在算法中,利用结点间的控制关系进行剪枝。从源顶点s出发,2条不同路径到达图G的同一顶点。由于两条路径的路长不同,因此可以将路长长的路径所对应的树中的结点为根的子树剪去。

 while (true) {
     for (int j = 1; j <= n; j++)
       if ((c[E.i][j]<inf)&&(E.length+c[E.i][j]<dist[j])) 
       {//顶点I和j间有边,且此路径长小于原先从原点到j的路径长 
         // 顶点i到顶点j可达,且满足控制约束
         dist[j]=E.length+c[E.i][j];
         prev[j]=E.i;
         // 加入活结点优先队列
         MinHeapNode<Type> N;
         N.i=j;
         N.length=dist[j];
         H.Insert(N);}
     try {H.DeleteMin(E);}         // 取下一扩展结点
     catch (OutOfBounds) {break;}  // 优先队列空
     }
} 

四、装载问题

4.1 问题描述

  有一批共n个集装箱要装上2艘载重量分别为C1和C2的轮船,其中集装箱i的重量为Wi,且在这里插入图片描述。装载问题要求确定是否有一个合理的装载方案可将这个集装箱装上这2艘轮船。如果有,找出一种装载方案。
  容易证明:如果一个给定装载问题有解,则采用下面的策略可得到最优装载方案
  (1)首先将第一艘轮船尽可能装满
  (2)将剩余的集装箱装上第二艘轮船


4.2 队列式分支限界法

  在算法的while循环中,首先检测当前扩展结点的左儿子结点是否为可行结点。如果是则将其加入到活结点队列中。然后将其右儿子结点加入到活结点队列中(右儿子结点一定是可行结点)。2个儿子结点都产生后,当前扩展结点被舍弃。
  活结点队列中的队首元素被取出作为当前扩展结点,由于队列中每一层结点之后都有一个尾部标记-1,故在取队首元素时,活结点队列一定不空。当取出的元素是-1时,再判断当前队列是否为空。如果队列非空,则将尾部标记-1加入活结点队列,算法开始处理下一层的活结点。

while (true) {
      // 检查左儿子结点
      if (Ew + w[i] <= c) // x[i] = 1
         EnQueue(Q, Ew + w[i], bestw, i, n);
      // 右儿子结点总是可行的
      EnQueue(Q, Ew, bestw, i, n); // x[i] = 0
      Q.Delete(Ew);     // 取下一扩展结点
      if (Ew == -1) {      // 同层结点尾部
         if (Q.IsEmpty()) return bestw;
         Q.Add(-1);        // 同层结点尾部标志
         Q.Delete(Ew);  // 取下一扩展结点
         i++;}                 // 进入下一层      }  }

4.3 算法的改进

  结点的左子树表示将此集装箱装上船,右子树表示不将此集装箱装上船。设bestw是当前最优解;ew是当前扩展结点所相应的重量;r是剩余集装箱的重量。则当ew+rbestw时,可将其右子树剪去,因为此时若要船装最多集装箱,就应该把此箱装上船。
  另外,为了确保右子树成功剪枝,应该在算法每一次进入左子树的时候更新bestw的值

// 检查左儿子结点
  Type wt = Ew + w[i];   // 左儿子结点的重量
      if (wt <= c) {     // 可行结点
         if (wt > bestw) bestw = wt;//提前更新bestw 
         // 加入活结点队列
         if (i < n) Q.Add(wt);
}
// 检查右儿子结点
      if (Ew + r > bestw && i < n)//右儿子剪枝 
          Q.Add(Ew);     // 可能含最优解
      Q.Delete(Ew);     // 取下一扩展结点


4.4 构造最优解

  为了在算法结束后能方便地构造出与最优值相应的最优解,算法必须存储相应子集树中从活结点到根结点的路径。为此目的,可在每个结点处设置指向其父结点的指针,并设置左、右儿子标志。

class QNode
 {QNode *parent;  // 指向父结点的指针
      bool LChild;        // 左儿子标志
      Type weight;       // 结点所相应的载重量

  找到最优值后,可以根据parent回溯到根节点,找到最优解。

// 构造当前最优解
for (int j = n - 1; j > 0; j--) {
      bestx[j] = bestE->LChild; 
      bestE = bestE->parent; 
}

4.5 优先队列式分支限界法

  解装载问题的优先队列式分支限界法用最大优先队列存储活结点表。活结点x在优先队列中的优先级定义为从根结点到结点x的路径所相应的载重量再加上剩余集装箱的重量之和。
  优先队列中优先级最大的活结点成为下一个扩展结点。以结点x为根的子树中所有结点相应的路径的载重量不超过它的优先级。子集树中叶结点所相应的载重量与其优先级相同
  在优先队列式分支限界法中,一旦有一个叶结点成为当前扩展结点,则可以断言该叶结点所相应的解即为最优解。此时可终止算法。


五、布线问题

5.1 算法思想

  解此问题的队列式分支限界法从起始位置a开始将它作为第一个扩展结点。与该扩展结点相邻并且可达的方格成为可行结点被加入到活结点队列中,并且将这些方格标记为1,即从起始方格a到这些方格的距离为1。
  接着,算法从活结点队列中取出队首结点作为下一个扩展结点,并将与当前扩展结点相邻且未标记过的方格标记为2,并存入活结点队列。这个过程一直继续到算法搜索到目标方格b或活结点队列为空时为止。即加入剪枝的广度优先搜索

Position offset[4];
//定义移动方向的相对位移
   offset[0].row = 0; offset[0].col = 1; // 右
   offset[1].row = 1; offset[1].col = 0; // 下
   offset[2].row = 0; offset[2].col = -1; // 左
   offset[3].row = -1; offset[3].col = 0; // 上
//设置边界的围墙
 for (int i = 0; i <= m+1; i++)
      grid[0][i] = grid[n+1][i] = 1; // 顶部和底部
 for (int i = 0; i <= n+1; i++)
      grid[i][0] = grid[i][m+1] = 1; // 左翼和右翼
for (int i = 0; i < NumOfNbrs; 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.Add(nbr);} 
         }

  找到目标位置后,可以通过回溯方法找到这条最短路径。


六、0-1背包问题

6.1 算法的思想

  首先,要对输入数据进行预处理,将各物品依其单位重量价值从大到小进行排列
  在下面描述的优先队列分支限界法中,节点的优先级由已装袋的物品价值加上剩下的最大单位重量价值的物品装满剩余容量的价值和确定
  算法首先检查当前扩展结点的左儿子结点的可行性。如果该左儿子结点是可行结点,则将它加入到子集树和活结点优先队列中。当前扩展结点的右儿子结点一定是可行结点,仅当右儿子结点满足上界约束时才将它加入子集树和活结点优先队列。当扩展到叶节点时为问题的最优值。


6.2 上界函数

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;    // 装填剩余容量装满背包
return b;                                       //b为上界函数
 while (i != n+1) {// 非叶结点
      // 检查当前扩展结点的左儿子结点
      Typew wt = cw + w[i];
      //分支限界搜索过程
      if (wt <= c) {// 左儿子结点为可行结点
         if (cp+p[i] > bestp) bestp = cp+p[i];
         AddLiveNode(up, cp+p[i], cw+w[i], true, i+1);}
         up = Bound(i+1);
      // 检查当前扩展结点的右儿子结点
      if (up >= bestp) // 右子树可能含最优解
           AddLiveNode(up, cp, cw, false, i+1);
         //   取下一个扩展节点(略)
}

  当前总重量,当前总价值,当前层数,当前能达到的极限情况下最大解(天花板)
  以当前能达到的极限情况下最大解的值由大到小排序成为优先队列,值最大的先出队列进行扩展,以此类推!
  背包容量为10,数量 4件。物品信息为w=[4,7,5,3];v=[40,42,25,12](排好序了)
在这里插入图片描述

  • 41
    点赞
  • 50
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 76
    评论
#include #include #include #include using namespace std; ifstream infile; ofstream outfile; class Node { friend int func(int*, int, int, int*); public: int ID; double weight;//物品的重量 }; bool comp1(Node a, Node b) //定义比较规则 { return a.weight > b.weight; } class Load; class bbnode; class Current { friend Load; friend struct Comp2; private: int upweight;//重量上界 int weight;//结点相应的重量 int level;//活结点在子集树中所处的层次 bbnode* ptr;//指向活结点在子集树中相应结点的指针 }; struct Comp2 { bool operator () (Current *x, Current *y) { return x->upweightupweight; } }; class Load { friend int func(int*, int, int, int*); public: int Max0(); private: priority_queue<Current*, vector, Comp2>H;//利用优先队列(最大堆)储存 int limit(int i); void AddLiveNode(int up, int cw, bool ch, int level); bbnode *P;//指向扩展结点的指针 int c;//背包的容量 int n;//物品的数目 int *w;//重量数组 int cw;//当前装载量 int *bestx;//最优解方案数组 }; class bbnode { friend Load; friend int func( int*, int, int, int*); bbnode* parent; bool lchild; }; //结点中有双亲指针以及左儿子标志 int Load::limit(int i) //计算结点所相应重量的上界 { int left,a; left= c - cw;//剩余容量 a = cw; //b是重量上界,初始值为已经得到的重量 while (i <= n && w[i] parent = P; b->lchild = ch; Current* N = new Current; N->upweight = up; N->weight = cw; N->level = level; N->ptr = b; H.push(N); } int Load::Max0() { int i = 1; P = 0; cw = 0; int bestw = 0; int up = limit(1); while (i != n + 1) { int wt = cw + w[i]; //检查当前扩展结点的左儿子结点 if (wt bestw) bestw =wt; AddLiveNode(up,wt, true, i + 1); } up = limit(i + 1); //检查当前扩展结点的右儿子结点 if (up >= bestw)//如果右儿子可行 { AddLiveNode(up,cw, false, i + 1); } Current* N = H.top(); //取队头元素 H.pop(); P = N->ptr; cw = N->weight; up = N->upweight; i = N->level; } bestx = new int[n + 1]; for (int j = n; j > 0; --j) { bestx[j] = P->lchild; P = P->parent; } return cw; } int func(int *w, int c, int n, int *bestx) //调用Max0函数对子集树的优先队列式进行分支限界搜索 { int W = 0; //初始化装载的总质量为0 Node* Q = new Node[n]; for (int i = 0; i < n; ++i) { Q[i].ID = i + 1; Q[i].weight = w[i+1]; W += w[i+1]; } if (W <= c)//如果足够装,全部装入 return W; sort(Q, Q + n, comp1); //首先,将各物品按照重量从大到小进行排序; Load K; K.w = new int[n + 1]; for (int j = 0; j < n; j++) K.w[j + 1] = w[Q[j].ID]; K.cw = 0; K.c = c; K.n = n; int bestp = K.Max0(); for (int k = 0; k < n; k++) { bestx[Q[k].ID] = K.bestx[k + 1]; } delete []Q; delete []K.w; delete []K.bestx; return bestp; } int main() { int*w,*Final; int c,n,i,best; infile.open("input.txt",ios::in); if(!infile) { cerr<<"open error"<>c; infile>>n; w=new int[n+1]; for(i=1;i>w[i]; infile.close(); Final = new int[n+1]; best = func( w, c, n, Final); outfile.open("output.txt",ios::out); if(!outfile) { cerr<<"open error"<<endl; exit(1); } outfile << best << endl; for (int i = 1; i <= n; ++i) { outfile<<Final[i]<<" "; } outfile.close(); return 0; }

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 76
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

TJUTCM-策士之九尾

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值