分支限界算法

分支限界算法

基本思想

◇ 概念

类似于回溯法，也是一种在问题的解空间树T上搜索问题解的算法。但在一般情况下，分支限界法与回溯法的求解目标不同。回溯法的求解目标是找出T中满足约束条件的所有解，而分支限界法的求解目标则是找出满足约束条件的一个解，或是在满足约束条件的解中找出使某一目标函数值达到极大或极小的解，即在某种意义下的最优解。
⑴ 分支搜索算法
所谓“分支”就是采用广度优先的策略，依次搜索E-结点的所有分支，也就是所有相邻结点，抛弃不满足约束条件的结点，其余结点加入活结点表。然后从表中选择一个结点作为下一个E-结点，继续搜索。
选择下一个E-结点的方式不同，则会有几种不同的分支搜索方式。
① FIFO搜索
② LIFO搜索
③ 优先队列式搜索
⑵ 分支限界搜索算法

◇ 分支限界法的一般过程

由于求解目标不同，导致分支限界法与回溯法在解空间树T上的搜索方式也不相同。回溯法以深度优先的方式搜索解空间树T，而分支限界法则以广度优先或以最小耗费优先的方式搜索解空间树T。
分支限界法的搜索策略是：在扩展结点处，先生成其所有的儿子结点（分支），然后再从当前的活结点表中选择下一个扩展对点。为了有效地选择下一扩展结点，以加速搜索的进程，在每一活结点处，计算一个函数值（限界），并根据这些已计算出的函数值，从当前活结点表中选择一个最有利的结点作为扩展结点，使搜索朝着解空间树上有最优解的分支推进，以便尽快地找出一个最优解。
分支限界法常以广度优先或以最小耗费（最大效益）优先的方式搜索问题的解空间树。问题的解空间树是表示问题解空间的一棵有序树，常见的有子集树和排列树。在搜索问题的解空间树时，分支限界法与回溯法对当前扩展结点所使用的扩展方式不同。在分支限界法中，每一个活结点只有一次机会成为扩展结点。活结点一旦成为扩展结点，就一次性产生其所有儿子结点。在这些儿子结点中，那些导致不可行解或导致非最优解的儿子结点被舍弃，其余儿子结点被子加入活结点表中。此后，从活结点表中取下一结点成为当前扩展结点，并重复上述结点扩展过程。这个过程一直持续到找到所求的解或活结点表为空时为止。

◇ 回溯法和分支限界法的一些区别

有一些问题其实无论用回溯法还是分支限界法都可以得到很好的解决，但是另外一些则不然。也许我们需要具体一些的分析——到底何时使用分支限界而何时使用回溯呢？
回溯法和分支限界法的一些区别：
方法对解空间树的搜索方式
存储结点的常用数据结构
结点存储特性常用应用
回溯法深度优先搜索堆栈活结点的所有可行子结点被遍历后才被从栈中弹出找出满足约束条件的所有解
分支限界法广度优先或最小消耗优先搜索队列、优先队列每个结点只有一次成为活结点的机会找出满足约束条件的一个解或特定意义下的最优解

例题

最小重量机器设计问题

设某一机器由n个部件组成，每一种部件都可以从m个不同的供应商处购得。设 wij 是从供应商j处购得的部件 i 的重量， cij 是相应的价格。试设计一个算法，给出总价格不超过 c 的最小重量机器设计。
主要实现代码如下：

import static org.junit.Assert.*;
 
import java.util.LinkedList;
import java.util.Queue;
import java.util.Scanner;
 
import javax.management.Query;
 
import org.junit.Test;
 
class Node{
    private int curw, curv;
    private int units_i;
    private int[] route;
    Node(){
        route = new int[110];
    }
    public Node(int curw, int curv, int units_i, int[] route) {
        super();
        this.curw = curw;
        this.curv = curv;
        this.units_i = units_i;
        this.route = route;
    }
    public int getCurw() {
        return curw;
    }
    public void setCurw(int curw) {
        this.curw = curw;
    }
    public int getCurv() {
        return curv;
    }
    public void setCurv(int curv) {
        this.curv = curv;
    }
    public int getUnits_i() {
        return units_i;
    }
    public void setUnits_i(int units_i) {
        this.units_i = units_i;
    }
    public int[] getRoute() {
        return route;
    }
    public void setRoute(int[] route) {
        this.route = route;
    }
    
}
 
public class 最小重量设计_分支限界法 {
    
    /*
     * 3 3
     * 3 2  1 4  5 6
     * 1 4  3 2  5 6
     * 5 6  3 2  1 4
     * 10
     * 
     * 答案是5
     */
    private static int w[][];
    private static int c[][];
    private static int n, m, cc;
    private static Scanner cin;
    private static Queue<Node>q;
    static{
        cin = new Scanner(System.in);
        q = new LinkedList<Node>();
    }
 
 
    public static void main(String[] args) {
        System.out.println("请输入部件数目以及供货商数量n和m：");
        n = cin.nextInt();
        m = cin.nextInt();
        w = new int[n][m];
        c = new int[n][m];
        
        System.out.println("n行代表n个部件，每行输入每个供货商供应此部件的重量以及价格：");
        for(int i = 0; i < n; i++){
            for(int j = 0; j < m; j++){
                w[i][j] = cin.nextInt();
                c[i][j] = cin.nextInt();
            }
        }
        System.out.println("请输入不超过的价格p：");
        cc = cin.nextInt();
        q.clear();
        int[] route = new int[110];
        Node e = new Node(0, 0, -1, route);
        q.add(e);
        int ans = Integer.MAX_VALUE;
        while(!q.isEmpty()){
            e = q.poll();
            for(int j = 0; j < m; j++){
                int i = e.getUnits_i() + 1;
                if(i >= n){
                    if(e.getCurv() <= cc){
                        if(ans > e.getCurw()){
                            ans = Math.min(ans, e.getCurw());
                            route = e.getRoute();
                        }
                    }
                    continue;
                }
                int curv = e.getCurv() + c[i][j];
                int curw = e.getCurw() + w[i][j];
                int[] lastroute = e.getRoute();
                int[] curroute = new int[110];
                for(int k = 0; k <= e.getUnits_i(); k++){
                    curroute[k] = lastroute[k];
                }
                curroute[i] = j;
                
                q.add(new Node(curw, curv, i, curroute));
            }
        }
        if(ans == Integer.MAX_VALUE){
            System.out.println("不能找出总价格不超过 c的最小重量机器的方案");
        }else{
            System.out.println("满足方案的最小重量是：" + ans);
            System.out.println("方案是：");
            for(int j = 0; j < n; j++){
                System.out.print(route[j] + " ");
            }
            
        }
    }
}