背包问题(01背包,完全背包,多重背包)(改)

上课没好好听,自己找资料琢磨的,加上了一点点自己的理解。

01背包

简单讲就是有一些物品(每种物品只有一个),告诉你他们的体积以及价值,求怎么你背包中的物品价值和最大。

首先就是打一个表,和dp有点像。

dp(i,j)中i表示第几见物品,j表示背包体积,总的来讲就是,将前i件物品放入体积为j的背包中可以获得的最大价值。

公式:dp[i][j]=max{dp[i-1][j],dp[i-1][j-v[i]]+value[i]},第一种表示第i件物品不取,第二种表示取。

例子如下:(体积,价值)(0,0)(2,3)(3,4)(4,5)(5,6)

d(i/j)012345678
0000000000
1003333333
2003447777
3003457899
40034578910
#include<bits/stdc++.h>
#define rep(i,a,n) for (int i=a;i<=n;i++)
using namespace std;
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
int dp[20][20]={},v[5]={0,2,3,4,5},value[5]={0,3,4,5,6};
int main(int argc, char** argv) {
	int bag=8;
	for(int i=0;i<=4;i++){
		for(int j=1;j<=bag;j++){
			if(j<v[i]){
				dp[i][j]=dp[i-1][j];
			}else{
				dp[i][j]=max(dp[i-1][j],dp[i-1][j-v[i]]+value[i]);
			}
		}
	}
	cout<<dp[4][bag];
	return 0;
}

然后我们发现很多数据是重复的,可以将表压缩成一行(节约空间),注意要从后往前遍历,以免数据覆盖。

优化后代码如下:

#include<bits/stdc++.h>
#define rep(i,a,n) for (int i=a;i<=n;i++)
using namespace std;
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
int dp[20][20]={},v[5]={0,2,3,4,5},value[5]={0,3,4,5,6};
int main(int argc, char** argv) {
	int bag = 8;
	for(int i=0;i<=4;i++){
		for(int j=bag;j>=v[i];j--){
			dp[j]=max(dp[j],dp[j-v[i]]+value[i]);
		}
	}
	cout<<dp[bag];
	return 0;
}

通过上面的方法,我们可以得知背包的最大价值,但不知道此时背包由哪些商品组成。

图中dp(4,8)!=dp(3,8),证明4号物品被选中了,那么我们回溯回上一个就是v(4-1,8-5)=v(3,3)。

由dp[i][j]=dp[i-1][j-v[i]]+value[i]可知,我们下一个开始比较的对象是dp[i-1][j-v[i]]。

要溯源的话就不可以压缩状态了哦

i/j012345678
0000000000
1003333333
2003447777
3003457899
40034578910
#include<bits/stdc++.h>
#define rep(i,a,n) for (int i=a;i<=n;i++)
using namespace std;
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
int v[5] = { 0 , 2 , 3 , 4 , 5 };
int value[5] = { 0 , 3 , 4 , 5 , 6 };
int bag = 8;
int dp[20][20] = {};
int item[5];

void findWhat(int i, int j) {	
	if (i >= 0) {
		if (dp[i][j] == dp[i - 1][j]) {//这件商品没被选中 
			item[i] = 0;
			findWhat(i - 1, j);
		}else if (j - v[i] >= 0 && dp[i][j] == dp[i - 1][j - v[i]] + value[i]) {
			item[i] = 1;
			findWhat(i - 1, j - v[i]);
		}
	}
}

int main(int argc, char** argv) {
	for (int i = 1; i <= 4; i++) {
		for (int j = 1; j <= bag; j++) {
			if (j < v[i]) //写习惯了状压的,不压缩的记得 < v[i] 的部分也要遍历
				dp[i][j]=dp[i-1][j];
			else
				dp[i][j]=max(dp[i-1][j],dp[i-1][j-v[i]]+value[i]);
		}
	}
	findWhat(4, 8);
	cout << dp[4][bag] << endl;
	for (int i = 0; i < 5; i++){
		cout << item[i] << ' ';
	}
	cout << endl;
	return 0;
}

XX背包(正正好装bag个,不多不少)

在背包的基础上初始化为-INF就好啦

#define INF 0x3f3f3f3f //10^9级的
memeset(dp,-INF,sizeof(dp));

完全背包

完全背包就是有一些物品(每种物品有无限个),告诉你他们的体积以及价值,求怎么你背包中的物品价值和最大。

其实和01背包差不多,但是注意要从前往后遍历(这样改变了后面的状态,就可以实现一个物品加入多次)。

01背包中一行表示一种物品,从后往前是为了避免前面改变了影响后面。而完全背包一行表示一种物品,可以有多个,所以前面加入了一个物品后面再加也OK。(救命,我这该死的表达能力)

#include<bits/stdc++.h>
#define rep(i,a,n) for (int i=a;i<=n;i++)
using namespace std;
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
int v[300],value[300],f[300010];
int bag,n;
int main(int argc, char** argv){
    cin>>bag>>n;
    for(int i=1; i<=n; i++){
        cin>>v[i]>>value[i];
    }
    for(int i=1; i<=n; i++){
    	for(int j=v[i];j<=bag; j++){
		    //注意此处,与01背包不同,这里为顺序,01背包为逆序
    	    f[j]=max(f[j],f[j-v[i]]+value[i]);
		}
	}        
    cout<<f[bag];
    return 0;
}

多重背包(朴素版)

又是差差不多呢。就是每种物品的数量改为有限个(题目中会输入哒)。

假设物品 i 有 num[i] 个,可以看作有 num[i] 种物品,每种的体积和质量分别为 k*v[i] , k*value[i] (1<=k<=num[i]),每种物品只能选一次(同01背包啦),再多遍历一维选 k 个物品就好啦。

注意:朴素版多重背包不能状压,因为状压后 f[j]=max(f[j],f[j-k*v[i]]+k*value[i]) 其中 f[j-k*v[i]] 部分可能已经被之前遍历的 k 改变过了。

#include<bits/stdc++.h>
#define rep(i,a,n) for (int i=a;i<=n;i++)
using namespace std;
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
int v[300],value[300],num[300],f[300010];
int bag,n;
int main(int argc, char** argv){
    cin>>bag>>n;
    for(int i=1; i<=n; i++){
        cin>>v[i]>>value[i]>>num[i];
    }
    for(int i=1; i<=n; i++){
    	for(int k=0; k<=num[i]; k++){
    		for(int j=v[i]*k; j<=bag; j++){	
    			f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*value[i]); 
			}
		}
	}        
    cout<<f[n][bag];
    return 0;
}

多重背包(二进制优化

简单讲就是,假设第i号物品有10件,那么把10转化为2进制就是1010,那我们讲i号分为1件,2件,4件,3件的小堆(2^0+2^1+2^2+剩下的3=10)(如果变成1248可能选到15件哦),通过小堆的组合我们可以得到所有(从0到10)的可能选择。

#include<bits/stdc++.h>
#define rep(i,a,n) for (int i=a;i<=n;i++)
using namespace std;
/* run this program using the console pauser or add your own getch, system("pause") or input loop */
int v1,value1,num1,sum,v[300],value[300],num[300],f[300010];
int bag,n,tmp;
int main(int argc, char** argv){
	cin>>bag>>n;
    for(int i=1; i<=n; i++){
        cin>>v1>>value1>>num1;
        int j=1, sum=0;
        while(sum+j<=num1){
        	sum=sum+j; 
        	v[++tmp]=v1*j; value[tmp]=value1*j;
        	j=j*2;
		}
		if(sum!=num1){
			v[++tmp]=v1*(num1-sum);
			value[tmp]=value1*(num1-sum);
		}
    }
    for(int i=1; i<=tmp; i++){
    	for(int j=v[i];j<=bag; j++){
    		// 从 j-v[i] 转移到 j
    	    f[j]=max(f[j],f[j-v[i]]+value[i]);
		}
	}        
    cout<<f[bag];
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值