【算法导论08】分支限界法-最大装载问题

08分支限界法-最大装载问题

问题描述:
已知有一个两个容量为C1,C2的轮船,现有集装箱n个,集装箱体积分别为W1-n,试问两个轮船能否装下所有集装箱,如果能装下,请设计算法计算最佳装载方案。

问题分析:
此问题为固定容量大小最大装载问题,如果想要判断最大装载能否装下所有集装箱,则首先考虑在C1最大装载的情况下,剩余集装箱能否全部装进C2,因此此问题转化为C1最大装载问题。由于每个集装箱都有装和不装两个选项,所以此问题的解答是一个解空间树,将设计用分支限界法来解决。

算法介绍:
分支限界法的核心思想是通过设计合适的限界函数,对全排列空间树的结点进行剪枝,减少不必要的解数量,以此提高算法的效率。限界函数的作用是判断该结点是否满足继续执行的必要,如果该条路径的解偏离最优方案,则不满足。例如在01背包问题中,如果加入某个商品结点,其最大价值相较于已经得出的价值变小,则进行剪枝。这里的最大价值便是此结点选择下的答案的上界(若是最小解则是下界)。因此,只要能计算出每个结点的上界,便可以对该结点进行剪枝判断。而某个结点的上界,往往需要将该结点往下所有分支全部计算后,才能得出,因此,在分支限界法中,通常使用估计上界(估计代价)来代替上界。估计上界是假设选择此结点后,下面所有的执行步骤都按照最优解来计算,如在01背包中,第i个商品装入后的估计上界为已经装入的总价值+后面所有商品在不考虑不可分情况下能装入的最大价值(此前已经对商品基于单位价值进行排序),即理论上在未来能装入的最大价值,将此作为上界函数。而剪枝则发生在活结点的上界<当前函数值时,表明不会再产生比当前函数值更大的数值。
分支限界法用队列实现,结点入队顺序与层次遍历类似,在扩展结点处,将其所有孩子结点入队,再从中选择一个作为新的扩展结点。传统的FIFBOO为按照先进先出的顺序选择扩展结点,而LCBB则利用优先队列,通过对活节点中每个结点的上界进行比较,优先选择最优上界作为新的扩展结点。在分支限界法中,每个结点通常有两个参数,一个为当前解(下界),另一个为上界,当上界小于当前解(下界)时,进行剪枝。在FIFBOO中,随着结点的选择,下界不断向上修正,直到上界=下界(当前解为最大解)或者上界<下界,进行剪枝。而在LCBB中,由于每次选择的均为最大上界,因此只要发生剪枝操作,算法结束运行,得出最优解。

笔记:
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

在这里插入图片描述
在这里插入图片描述

最大装载问题代码:

//FLFOBB
int MaxLoading(int w[],int c,int n){//物品重量,最大容量,物品个数
	Queue Q;
	Q.enqueue(-1);
	int i = 1;
	int ew = 0,r = 0,bestw = 0;//截止拓展结点处载重量,剩余载重量,最优载重量
	for(int j=2;j<=n;j++){
		r+ = w[j];
	}	
	while(true){
		int wt = ew+w[i];//左子树结点重量,即选了w[i]的载重量
		if(wt<=c){//判断左子树可行性
			if(wt>bestw)//该解是否大于最优解
				bestw = wt;
			if(i<n)
				Q.enqueue(wt);	
		}
		if(ew+r>bestw&&i<n)//判断右子树可行性
			Q.equeue(ew)
		Q.dequeue(ew);//取下一个扩展结点			
		if(ew = -1){
			if(Q.isempty())
				return bestw;
			Q.enqueue(-1);
			Q.dequeue(ew);
			i++;
			r -= w[i];
		}
	}
}
#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、付费专栏及课程。

余额充值