背包问题求方案数(21-21)

52 阅读2分钟

开启掘金成长之旅!这是我参与「掘金日新计划 · 2 月更文挑战」的第 21 天,点击查看活动详情

背包问题求方案数

有 NN 件物品和一个容量是 VV 的背包。每件物品只能使用一次。

第 ii 件物品的体积是 viv_i,价值是 wiw_i

求解将哪些物品装入背包,可使这些物品的总体积不超过背包容量,且总价值最大。

输出 最优选法的方案数。注意答案可能很大,请输出答案模 109+710^9+7 的结果。

输入格式

第一行两个整数,NVN,V, 用空格隔开,分别表示物品数量和背包容积。

接下来有 NN 行,每行两个整数 vi,wiv_i,w_i,用空格隔开,分别表示第 ii 件物品的体积和价值。

输出格式

输出一个整数,表示 方案数 模 109+7109+7 的结果。

数据范围

0<N,V10000<N,V≤1000 0<vi,wi10000<vi,wi≤1000

输入样例

4 5
1 2
2 4
3 4
4 6

输出样例:

2

题目分析

相较于传统01背包,本题多出一个求方案数。

首先定义 f[i]f[i] 用来表示背包体积为 ii 时的最佳方案的总价值,g[i]g[i] 用来表示背包体积为 ii 时总价值为最佳的方案数。

在进行 ff 数组的更新时,我们要以选或不选当前物品的总价值是否改变进行分开讨论。

初始化所有的 g[i]g[i]11,表示价值为 00 时的方案数为 11

如果选择第 ii 件物品时的总价值更大,更新 f[i]f[i]f[jv]+wf[j-v]+w

如果总价值相等,则当前 ii 的最大方案数加上 f[jv]f[j-v] 种。

同时按照题意对 modmod 取模。

时间复杂度为 O(nm)O(nm)

Accept代码 O(nm)

#include <bits/stdc++.h>

using namespace std;

const int N = 1010, mod = 1e9 + 7;

int n, m;
int f[N], g[N];

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

    memset(f, -0x3f, sizeof f);
    f[0] = 0;
    g[0] = 1;

    for (int i = 0; i < n; i ++ )
    {
        int v, w;
        cin >> v >> w;
        for (int j = m; j >= v; j -- )
        {
            int maxv = max(f[j], f[j - v] + w);
            int s = 0;
            if (f[j] == maxv) s = g[j];
            if (f[j - v] + w == maxv) s = (s + g[j - v]) % mod;
            f[j] = maxv, g[j] = s;
        }
    }

    int res = 0;
    for (int i = 1; i <= m; i ++ )
        if (f[i] > f[res])
            res = i;

    int sum = 0;
    for (int i = 0; i <= m; i ++ )
        if (f[i] == f[res])
            sum = (sum + g[i]) % mod;

    cout << sum << endl;

    return 0;
}