0-1背包问题(回溯法解决)

 

0-1背包问题(回溯法解决)

给定一个物品集合s={1,2,3,…,n},物品i的重量是wi,其价值是vi,背包的容量为W,即最大载重量不超过W。在限定的总重量W内,我们如何选择物品,才能使得物品的总价值最大。

输入
第一个数据是背包的容量为c(1≤c≤1500),第二个数据是物品的数量为n(1≤n≤50)。接下来n行是物品i的重量是wi,其价值为vi。所有的数据全部为整数,且保证输入数据中物品的总重量大于背包的容量。

输出:对每组测试数据,输出装入背包中物品的最大价值。

令cw(i)表示目前搜索到第i层已经装入背包的物品总重量,即部分解(x1, x2 , …, xi)的重量:

在这里插入图片描述
对于左子树, xi =1 ,其约束函数为:
在这里插入图片描述

若constraint(i)>W,则停止搜索左子树,否则继续搜索。

对于右子树,为了提高搜索效率,采用限界函数Bound(i)剪枝。
令cv(i)表示目前到第i层结点已经装入背包的物品价值:

在这里插入图片描述
令r(i)表示剩余物品的总价值:
在这里插入图片描述

则限界函数Bound(i)为:

在这里插入图片描述

假设当前最优值为bestv,若Bound(i)<bestv,则停止搜索第i层结点及其子树,否则继续搜索。
显然r(i)越小, Bound(i)越小,剪掉的分支就越多(在高层剪枝,剪掉的分支越多。从而能加快搜索速度)。
为了构造更小的r(i) ,将物品以单位重量价值比di=vi/wi递减的顺序进行排列(贪心策略):
d1≥d2≥… ≥dn
对于第i层,背包的剩余容量为W-cw(i),采用贪心算法把剩余的物品放进背包。

在这里插入图片描述

为什么没有按照重量进行排序呢?

0-1背包问题的目标是:
在不超重的前提下计算背包内物品的最大价值
因此,贪心策略是:
按照单位价值由大到小进行处理
物品的价值会对限界函数产生影响
物品的重量会对约束函数产生影响
由于根据物品价值制定贪心策略
所以,通过单位价值排序,加速剪枝

#include<iostream>
#include<algorithm>
using namespace std;
#define NUM 100
int c;			//背包的容量		
int n;			//物品的数量
int cw;			//当前背包内物品的重量
int cv;			//当前背包内物品的总价值
int bestv;		//当前最优价值



//描述每个物品的数据结构
struct Object
{
public:
	int w;		//物品的重量
	int v;		//物品的价值
	double d;	//物品的单位价值
public:
	double getd()
	{
		return d;
	}
};		//物品数组

Object Q[NUM];
void backtrack(int i);
int Bound(int i);
bool cmp(Object, Object);


//物品的单位价值重量比是在输入数据时计算的


int main()
{
	cin >> c >> n;


	for (int i = 0; i < n; i++)
	{
		//物品的单位价值重量比是在输入数据时计算的
		scanf_s("%d%d", &Q[i].w, &Q[i].v);
		Q[i].d = 1.0 * Q[i].v / Q[i].w;

	}

	sort(Q, Q + n, cmp);

	backtrack(1);

	cout << bestv;

}








//以物品单位价值重量比递减排序的因子:
bool cmp(Object a, Object b)
{
	if (a.d>= b.d)
		return true;
	else
		return false;
	
}


void backtrack(int i)
{
	//到达叶子节点时更新最优值
	if (i + 1 > n)
	{
		bestv = cv;
		return;
	}

	//进入左子树搜索
	if (cw + Q[i].w <= c)
	{
		cw += Q[i].w;
		cv += Q[i].v;
		backtrack(i + 1);
		cw -= Q[i].w;
		cv -= Q[i].v;

	}

	//进入右子树搜索
	if (Bound(i + 1) > bestv)
		backtrack(i + 1);
}

int Bound(int i)
{
	int cleft = c - cw;				//背包剩余的容量
	int b = cv;						//上界
	//尽量装满背包
	while (i < n && Q[i].w <= cleft)
	{
		cleft -= Q[i].w;
		b += Q[i].v;
		i++;
	}

	//剩余的部分空间也装满
	if (i < n)
		b += 1.0 * cleft * Q[i].v / Q[i].w;
	return b;
}
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页