引入
关于01背包,其实就很简单了,可以看看我的另一篇文章完全背包和01背包对比;那如果把条件换成:求背包刚好装满的时候的最大值,阁下又该如何应对;
题解
我们知道01背包的转移方程如下:
dp[i][j] = max(dp[i - 1][j], dp[i - 1][j - v[i]] + w);
这里有一种很巧妙的处理方式,即把不能刚好填满背包的值赋为-INF(即标记为无效的状态),例如dp[1][0] = 0而不是-INF,因为其确实填满了背包;这样处理过后,我们只需要判断 dp[n][m] 是否大于0即可判断是否能填满背包并且得到能填满背包所有价值的最大值了;
状态的转移
这里有几个需要考虑的问题:
- 是否会出现漏掉某一个数的情况: 每一个 dp[i][j] 代表的是一个只考虑前i个物品,背包只有j个空间的 所有选法 的集合,所以每一个状态都会继承上一个集合,并不会漏掉情况;
- 状态可能由无效的状态转移过来吗: 实际上不会,当 dp[i - 1][j] 和 dp[i - 1][j - v[i]] 中有一个无效的时候的时候max会排除无效的那一个(因为太小了),这样能保证只要这两个有一个是有效的话,那么 dp[i][j] 就一定会有效;
综上所述,当求出来所有的状态之后,若 dp[i][j] 是有效的,说明背包能恰好背装满且有最大值;因此代码如下:
// web:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> PII;
typedef pair<ll, ll> PLL;
const int N = 1e5 + 10;
const int INF = 1e9;
const ll INFL = 1e18;
signed main()
{
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
int n, m;
cin >> n >> m;
vector<int> v(n + 1), w(n + 1);
for (int i = 1; i <= n; i++)
cin >> v[i] >> w[i];
vector<vector<int>> dp(n + 1, vector<int>(m + 1, -INF));
for (int i = 0; i <= n; i++)
dp[i][0] = 0;
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++)
if (j - v[i] >= 0)
dp[i][j] = max(dp[i - 1][j], dp[i - 1][j - v[i]] + w[i]) > 0
? max(dp[i - 1][j], dp[i - 1][j - v[i]] + w[i])
: -INF;
cout << ((dp[n][m] > 0) ? dp[n][m] : -1) << endl;
return 0;
}
另一道例题[P5020 [NOIP2018 提高组] 货币系统]
P5020 [NOIP2018 提高组] 货币系统;这道题其实就可以转化为在 (b, m) 中最少几个数可以凑出 (a, n) 这个很好想到;再转化一下思路,就是 哪几个容量的背包可以被刚好装满,这样就转化成了一个完全背包刚好装满的问题;
就和01背包刚好装满问题类似,关于01背包和完全背包问题的对比可以看看我之前提到的我的另一篇文章;看看码:
// web:https://www.luogu.com.cn/problem/P5020
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> PII;
typedef pair<ll, ll> PLL;
const int N = 1e5 + 10;
const int INF = 1e9;
const ll INFL = 1e18;
signed main()
{
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
int t;
cin >> t;
while (t--)
{
int n;
cin >> n;
vector<int> a(n + 1);
int mxa = -1;
for (int i = 1; i <= n; i++)
cin >> a[i], mxa = max(mxa, a[i]);
vector<int> dp(25000 + 10, -INF);
dp[0] = 0;
for (int i = 1; i <= n; i++)
for (int j = 1; j <= mxa; j++)
if (j - a[i] >= 0)
dp[j] = max(dp[j], dp[j - a[i]] + 1);
int ans = 0;
for (int i = 1; i <= n; i++)
if (dp[a[i]] == 1)
ans++;
cout << ans << endl;
}
return 0;
}