435-分支限界算法-集装箱装载问题

集装箱装载问题

有一批共n个集装箱要装上2艘载重量分别为c1,c2的轮船,其中集装箱的重量为wi,且要求确定是否有一个合理的装载方案可将这n个集装箱装上这2艘轮船。

思路:
我们要尽最大可能把一个轮船塞得越满越好,让c1的轮船最大化装载。然后把剩下的集装箱装上另一个轮船。

我们定义

int w[] = { 12, 8, 15 };//集装箱的重量
const int n = sizeof(w) / sizeof(w[0]);//集装箱的数量
int c = 27;//轮船的容量
int cw = 0;//已选择物品的重量
int bestw = 0;//记录最优的装载量

我们的左孩子表示选择(1),右孩子表示不选择(0),我们先处理左孩子。

因为是层序遍历,一层一层的遍历,我们先处理第一层:
每一层是从左向右处理的
第一层(从左向右)
在这里插入图片描述
处理第二层(从左向右)
在这里插入图片描述
处理最后一层(从左向右)
在这里插入图片描述
我们看出,最优解是这个红色框起来的叶子节点
在这里插入图片描述
在这里插入图片描述
我们是广度优先遍历,找到最优解。但是不能获取到哪些集装箱进了轮船。
所以我们在扩展孩子节点的时候要记录一下父节点是谁,通过parent找见。还得知道是否选择该节点,所以还要增加个标识判断!!!

分支限界解题代码

#include <iostream>
#include <queue>
using namespace std;


int w[] = { 12, 8, 15 };//集装箱的重量
const int n = sizeof(w) / sizeof(w[0]);//集装箱的数量
int c = 20;//轮船的容量
int cw = 0;//已选择物品的重量
int bestw = 0;//记录最优的装载量

//描述节点类型
struct Node 
{
	Node(int w, int l, Node *p, bool left) 
	{
		weight = w;
		level = l;
		parent = p;
		isleft = left;
	}
	int weight;//从根节点到当前节点所选择物品的总重量
	int level;//当前节点的层数
	Node *parent;//记录当前节点的父节点
	bool isleft;//记录当前节点是否被选择,true就是被往左选择
};

//广度优先遍历子集树的FIFO队列(先进先出队列)
queue<Node*> que;
int i = 0;//表示起始的层数,从根节点开始
Node *bestnode = nullptr;
int r = 0;//记录未被处理的物品的总重量

//添加活结点到队列当中
void addLiveNode(int w, int level, Node *parent, bool isleft) 
{
	Node *node = new Node(w, level, parent, isleft);
	que.push(node);

	//在最后一层,记录最优值节点
	if (level == n && w == bestw) 
	{
		bestnode = node;
	}
}
//节点的上界函数, 返回未被处理的物品的总重量
int maxBound(int level) 
{
	int s = 0;
	for (int j = level + 1; j < n; ++j)//当前level是正在处理的节点 
	{
		s += w[j];
	}
	return s;
}
int main()
{
	Node *node = nullptr; //初始状态,根节点无父节点
	while (i < n) 
	{   //当前节点是第i层 0   开始扩展node节点
		//处理左孩子,表示选择i节点
		int wt = cw + w[i];//已选择的物品加上即将选择的物品 
		if (wt <= c)//选择物品i以后,其总重量不能超过轮船的容量 
		{ 
			if (wt > bestw)//更新一下 bestw
			{
				bestw = wt;
			}

			//活结点孩子入队列
			//que.push(Node(cw + w[i], i + 1));
			addLiveNode(cw+w[i], i+1, node, true);
		}

		//处理右孩子,表示不选择i节点
		//que.push(Node(cw, i + 1));
		r = maxBound(i);//求第i个节点的重量值上界  
		if (cw + r >= bestw) 
		{ //>=这里的=操作不能少,否则无法选择到叶子节点上 c=20
			addLiveNode(cw, i + 1, node, false);
		}
		
		//处理完i节点后,它成为死节点,然后出队
		node = que.front();
		que.pop();
		//恢复cw和i的值,表示从i节点跳到广度遍历的下一个节点了
		//因为是层序遍历,所以要恢复一下
		cw = node->weight; 
		i = node->level;
	}

	cout << bestw << endl;
	int bestx[n] = { 0 };
	for (int j = n - 1; j >= 0; --j) 
	{
		bestx[j] = bestnode->isleft ? 1 : 0;
		bestnode = bestnode->parent;
	}
	for (int v : bestx) 
	{
		cout << v << " ";
	}
	cout << endl;

	return 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
发出的红包

打赏作者

林林林ZEYU

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

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

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

打赏作者

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

抵扣说明:

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

余额充值