poj - 3624（01背包的常数优化）

poj-3624

 

01背包：

（1）二维数组表示dp[i][j]表示对第i个物品进行操作后，j大小的背包体积所能得到的最大价值。

所以每次对物品i操作时的状态转移方程为：dp[i][j] = MAX(dp[i-1][j],dp[i-1][j-a[i]]+b[i]);

表示选与不选物品i。

伪代码：

for(i = 1~n)

    for（j = ai~m）

      dp[i][j] = MAX(dp[i-1][j],dp[i-1][j-a[i]]+b[i]);

 

（2）一维数组的优化：

因为对j的更新有先后顺序，所以我们可以先更新m，就是j = （m~ai）的情况下进行更新，

这样是先更新后面的dp[j-a[i]]+b[i] , dp[j]等同于dp[i-1][j-a[i]]+b[i],dp[i-1][j],所以可以将化简为以为数组。

 

（3）常数的优化

将j的范围更改为（m->MAX(ai,m-∑ak（k是i~n）)）

含义是如果取完剩下的i到n个物品，剩下的空间还是比ai大（就是还有空间留给下一个将要装入背包的物体），

那就不用更新更左边的体积的最大价值（因为接下来的更新用不到那么左边的价值），

不过只对背包求最大价值这题成立，其他情况下可能不成立（eg：求多少种价值）。

 

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn = 1e5+10;
int a[maxn],b[maxn],dp[maxn]={0},m,n;
int MAX(int x,int y){
	return x>y?x:y;
}
int main(void){
	scanf("%d%d",&n,&m);
	int sum = 0;
	for(int i=1;i<=n;i++) scanf("%d%d",&a[i],&b[i]),sum+=a[i];
	for(int i=1;i<=n;i++){
		int tp = MAX(a[i],m-sum);
		for(int j=m;j>=tp;j--)
		if(dp[j-a[i]]+b[i]>dp[j]){
			dp[j] = dp[j-a[i]]+b[i];
		}
		sum-=a[i];
	}
	int ans = dp[1];
	for(int i=1;i<=m;i++) ans = MAX(ans,dp[i]);
	printf("%d\n",ans);
	return 0;
}