POJ3634-动态规划01背包

Charm Bracelet

Time Limit: 1000MS		Memory Limit: 65536K
Total Submissions: 14893		Accepted: 6809

Description

Bessie has gone to the mall's jewelry store and spies a charm bracelet. Of course, she'd like to fill it with the best charms possible from the N (1 ≤ N ≤ 3,402) available charms. Each charm i in the supplied list has a weightW_i (1 ≤ W_i ≤ 400), a 'desirability' factor D_i (1 ≤ D_i ≤ 100), and can be used at most once. Bessie can only support a charm bracelet whose weight is no more than M (1 ≤ M ≤ 12,880).

Given that weight limit as a constraint and a list of the charms with their weights and desirability rating, deduce the maximum possible sum of ratings.

Input

* Line 1: Two space-separated integers: N and M
* Lines 2..N+1: Line i+1 describes charm i with two space-separated integers: W_i and D_i

Output

* Line 1: A single integer that is the greatest sum of charm desirabilities that can be achieved given the weight constraints

Sample Input

4 6
1 4
2 6
3 12
2 7

Sample Output

23

Source

USACO 2007 December Silver

题目大意：说Bessie要选N个珠宝，总重量不能超过M，每个珠宝都有自己的重量和价值，求最大的总价值。

解题思路：最开始我用f[i][j]来代表前i个珠宝放入载重量为j的背包的最大价值,w[i]表示第i个珠宝的重量，用v[i]表示第i个珠宝的价值。考虑第i个珠宝，只有放和不放两种可能，那么如果不放第i个珠宝，价值为f[i-1][j], 如果放第i个珠宝，价值为f[i-1][j-w[i]]+v[i]

要使总价值最大，那么有f[i][j]=max(f[i-1][j],f[i-1][j-w[i]]+v[i])，思考出这个方程后，接着就是要填充它。

刚开始做dp，我最开始对这个方程还是挺迷茫的，不知道怎么写循环。后来想了一下，这个方程的重要思想，或者说dp的重要思想就是借助状态之间的关联性，用初始状态推出末状态，也就是我们所需的状态来。因此，应该按照这个方程的思路，依次来填充这个方程，首先是珠宝数，应该是从1->N，然后是载重量，应该是从0->M,因此循环如下：

	for(int i=1;i<=N;++i)
			for(int j=0;j<=M;++j)
				if(j>=weight[i])
					f[i][j]=max(f[i-1][j],f[i-1][j-weight[i]]+value[i]);
				else f[i][j]=f[i-1][j];

抱着侥幸的心理我提交了下，结果超内存了。算算消耗的空间为 12,880*3,402=43817760，我又开的int，每个占4字节，因此消耗的总内存约为17wkb，题目要求在6w kb以下，看来还得改进。

搜索了一下资料，自己又想了好久，想到了其实f[i][j]中的[i]是可以去掉的。因为观察f[j]=max(f[j],f[j-w[i]]+v[i])，只需要等号右边的f[j]保持在之前那个状态，这样算出来的等号左边的f[j]就是正确的,而且j循环也需要逆向循环。

AC代码如下

11325798

TSERROF

3624

Accepted

324K

329MS

C++

391B

2013-03-08 17:14:54

using namespace std;
#define MAXN 3403
int N,M, weight[MAXN], value[MAXN], f[12881];
int main()
{
	while(cin>>N>>M)
	{
		for(int i=1;i<=N;++i)
			cin>>weight[i]>>value[i];

		memset(f,0,sizeof(f));

		for(int i=1;i<=N;++i)
			for(int j=M;j>=0;--j)
				if(j>=weight[i])
					f[j]=max(f[j],f[j-weight[i]]+value[i]);

		cout<<f[M]<<endl;
	}
	return 0;
}