分支限界法——装载问题

上一节我们了解了什么是分支限界法,按照惯例我们要根据例题来总结分支限界法的思想和模式。

简单来说我们学习到分支限界法实际上是对解空间树的广度优先遍历,而其核心可以拆解为分支和限界两个词语,分支代表按照广度展开的可能分支,而限界代表我们所规定的剪枝的上界和下界条件。将可行的结点加入活结点表,不可行的则剪枝,直至我们找到符合输出条件的分支。


让我们回顾一下装载问题, 

 这是一个子集树问题,还记得我们用回溯法来解决的时候,就是依次遍历全部的解空间并按照限定条件来剪枝。现在我们要用分支限界法来解决,那么对于这个题目来说我们要搞清楚分支,从根节点开始就是c1c2c3,如果选择了c1那么分支就是c2c3,选择c1c2分支就是c3。当然可以不选。限界就是规定总重量不超过c并且要求记录最大的bestw。

 (注意MaxLoading不是只执行一次,而是需要遍历整个子集树)

 改进了之后增加了一个剪枝条件,对于当前重量Ew+剩余重量r<=bestw的进行剪枝,因为全加上也无法超过bestw,没必要再遍历。


以上代码是队列式的分支限界法,还有一种优先队列式的分支限界法。

找了网上的一张图

 简单来说,对于一个货物,我们的选择只有选中或者不选中。想要尽可能装重一点的货物,那么其上界就要尽可能大。在图中所给的例子中我们发现,up=所有选中的货物重量+可供选择的货物重量,它所代表的是最大重量的上界,那么如果我们要寻找最优解,肯定要从上界大的子树里寻找。这个up的值当且仅当我们明确某货物不会被选中,才会减去它的重量,比如第二层发现8+6=14>c1=12。在超重剪枝的情况下我们就会直接走到右子树。而在第三层,虽然左子树上界较大,但是我们发现最终up较大值是在右子树当up=11时,这也是为什么我们不再遍历A到C这条路径,因为A到C这里up也是11,不大于我们记录的bestw=11(bestw仅当为叶子节点才更新),所以剪枝了。

 


从这两题我们总结一下分支限界法的一般规律:

分支限界法的概念是分支和限界,分支即我们如何选取拓展节点,限界即我们需要限定的上界和下界。判断是否得到了解的条件较多:

比如如果我们的解达到了上界,说明已经是最优了,那么就是我们所求的解。

如果某条分支的上界已经低于总体下界,那么说明已经不在解范围内,我们就需要剪枝。

而优先队列实际上表达的是对界限的动态更新,每层的上界或者下界都是实时更新着的,我们每次所遍历的最优解都可以作为新的上界或下界,这样做的好处是将一部分的路径可以直接剪枝,使得计算量更少。

#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; }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值