分支限界---->装载问题

装载问题

一、问题描述

有两艘船,n个货箱。第一艘船的载重量是c1,第二艘船的载重量是c2,wi是货箱i 的重量,且w1+w2+……+wn≤c1+c2。确定是否有一种可将所有n个货箱全部装船的方法。若有的话,找出该方法。

例子:

有两艘船,3个货箱。第一艘船的载重量是50,第二艘船的载重量是50。

         当w =[10,40,40]时,有解;

         当w =[20,40,40]时,无解。

二、问题是否可解的判断

当第一艘船的可行最大装载为bestw时,若w1+w2+……+wn-bestw≤c2,则问题有解。即

总的重量 - 第一艘船实际能够装载的最大重量  < 第二艘船的载重量 时

因此,此问题可以转化为第一艘船的装载问题。

三、分析

  • “装载上界”,就是假设装载当前物品以后的所有物品。
  • 若当前分支的“装载上界”,比现有的最大装载小,则该分支就无需继续搜索。这样就缩小搜索范围,提高了搜索效率。
  • 优先队列中结点优先级常规定为一个与该结点相关的数值p,它一般表示其接近最优解的程度,本例就以当前结点所在分支的装载上界为优先值

四、例子

W={10,30,50},C1=60, 所构成的子集树如下图所表示:

4.1、FIFO限界搜索过程

1)初始队列中只有结点A;

2)结点A变为E-结点扩充B入队,bestw=10;

  结点C的装载上界为30+50=80> bestw,也入队;

3)结点B变为E-结点扩充D入队,bestw=40;

  结点E的装载上界为60> bestw,也入队;

4)结点C变为E-结点扩充F入队,bestw仍为40;

  结点G的装载上界为50> bestw,也入队;

5)结点D变为E-结点,叶结点H超过容量,

  叶结点I的装载为40,bestw仍为40;

6)结点E变为E-结点,叶结点J装载量为60,bestw为60;

  叶结点K被剪掉;

7)结点F变为E-结点,叶结点L超过容量,bestw为60;

  叶结点M被剪掉;

8)结点G变为E-结点,叶结点N、O都被剪掉;

  此时队列空算法结束。

4.2、LC-搜索的过程

1) 初始队列中只有结点A;

2) 结点A变为E-结点扩充B入堆,bestw=10;

结点C的装载上界30+50=80>bestw,入堆;堆中B上界为90在优先队列首。

3) 结点B变为E-结点扩充D入堆,bestw=40;

结点E的装载上界60>bestw,入堆;堆中D上界为90为优先队列首。

4) 结点D变为E-结点,叶结点H超过重量,叶结点I的装载为40,入堆,bestw仍为40;此时堆中C上界为80为优先队列首。

5) 结点C变为E-结点扩充F入堆,bestw仍为40;

结点G的装载上界50> bestw,入堆;此时堆中E上界为60为优先队首

6) 结点E变为E-结点,叶结点J装载量为60入堆,bestw变为60;

叶结点K上界10<bestw被剪掉;此时堆中J上界为60为优先队列首。

7) 结点J变为E-结点,扩展的层次为4算法结束。

虽然此时堆并不空,但可以确定已找到了最优解。

优先队列限界搜索解空间的过程是:A-B-D-C-E-J

  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
#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、付费专栏及课程。

余额充值