Question

2902. 和带限制的子多重集合的数目

English Version

题目描述

给你一个下标从 0 开始的非负整数数组 nums 和两个整数 l 和 r 。

请你返回 nums 中子多重集合的和在闭区间 [l, r] 之间的 子多重集合的数目 。

由于答案可能很大，请你将答案对 109 + 7 取余后返回。

子多重集合 指的是从数组中选出一些元素构成的 无序 集合，每个元素 x 出现的次数可以是 0, 1, ..., occ[x] 次，其中 occ[x] 是元素 x 在数组中的出现次数。

注意：

• 如果两个子多重集合中的元素排序后一模一样，那么它们两个是相同的 子多重集合 。
• 空 集合的和是 0 。

 

示例 1：

输入：nums = [1,2,2,3], l = 6, r = 6
输出：1
解释：唯一和为 6 的子集合是 {1, 2, 3} 。

示例 2：

输入：nums = [2,1,4,2,7], l = 1, r = 5
输出：7
解释：和在闭区间 [1, 5] 之间的子多重集合为 {1} ，{2} ，{4} ，{2, 2} ，{1, 2} ，{1, 4} 和 {1, 2, 2} 。

示例 3：

输入：nums = [1,2,1,3,5,2], l = 3, r = 5
输出：9
解释：和在闭区间 [3, 5] 之间的子多重集合为 {3} ，{5} ，{1, 2} ，{1, 3} ，{2, 2} ，{2, 3} ，{1, 1, 2} ，{1, 1, 3} 和 {1, 2, 2} 。

 

提示：

• 1 <= nums.length <= 2 * 104
• 0 <= nums[i] <= 2 * 104
• nums 的和不超过 2 * 104 。
• 0 <= l <= r <= 2 * 104

解法

方法一

class Solution:
  def countSubMultisets(self, nums: List[int], l: int, r: int) -> int:
    kMod = 1_000_000_007
    # dp[i] := # of submultisets of nums with sum i
    dp = [1] + [0] * r
    count = collections.Counter(nums)
    zeros = count.pop(0, 0)

    for num, freq in count.items():
      # stride[i] := dp[i] + dp[i - num] + dp[i - 2 * num] + ...
      stride = dp.copy()
      for i in range(num, r + 1):
        stride[i] += stride[i - num]
      for i in range(r, 0, -1):
        if i >= num * (freq + 1):
          # dp[i] + dp[i - num] + dp[i - freq * num]
          dp[i] = stride[i] - stride[i - num * (freq + 1)]
        else:
          dp[i] = stride[i]

    return (zeros + 1) * sum(dp[l : r + 1]) % kMod

class Solution {
  static final int MOD = 1_000_000_007;
  public int countSubMultisets(List<Integer> nums, int l, int r) {
    Map<Integer, Integer> count = new HashMap<>();
    int total = 0;
    for (int num : nums) {
      total += num;
      if (num <= r) {
        count.merge(num, 1, Integer::sum);
      }
    }
    if (total < l) {
      return 0;
    }
    r = Math.min(r, total);
    int[] dp = new int[r + 1];
    dp[0] = count.getOrDefault(0, 0) + 1;
    count.remove(Integer.valueOf(0));
    int sum = 0;
    for (Map.Entry<Integer, Integer> e : count.entrySet()) {
      int num = e.getKey();
      int c = e.getValue();
      sum = Math.min(sum + c * num, r);
      // prefix part
      // dp[i] = dp[i] + dp[i - num] + ... + dp[i - c*num] + dp[i-(c+1)*num] + ... + dp[i %
      // num]
      for (int i = num; i <= sum; i++) {
        dp[i] = (dp[i] + dp[i - num]) % MOD;
      }
      int temp = (c + 1) * num;
      // correction part
      // subtract dp[i - (freq + 1) * num] to the end part.
      // leves dp[i] = dp[i] + dp[i-num] +...+ dp[i - c*num];
      for (int i = sum; i >= temp; i--) {
        dp[i] = (dp[i] - dp[i - temp] + MOD) % MOD;
      }
    }
    int ans = 0;
    for (int i = l; i <= r; i++) {
      ans += dp[i];
      ans %= MOD;
    }
    return ans;
  }
}

class Solution {
public:
  int countSubMultisets(const vector<int> &nums, int l, int r) {
    int cnt[20001] = {};
    int memo[20001] = {};
    const int mod = 1000000007;
    for (int n : nums) {
      ++cnt[n];
    }
    fill_n(memo, cnt[1] + 1, 1);
    for (int n = 2, total = cnt[1]; n <= r; ++n) {
      if (!cnt[n]) {
        continue;
      }
      int top = (cnt[n] + 1) * n;
      total += n * cnt[n];
      for (int i = n, ii = min(total, r); i <= ii; ++i) {
        memo[i] = (memo[i] + memo[i - n]) % mod;
      }
      for (int i = min(total, r); i >= top; --i) {
        memo[i] = (mod + memo[i] - memo[i - top]) % mod;
      }
    }
    return accumulate(memo + l, memo + r + 1, 0LL) * (cnt[0] + 1) % mod;
  }
};

func countSubMultisets(nums []int, l int, r int) int {
  multiset := make(map[int]int)
  for _, num := range nums {
    multiset[num]++
  }
  mem := make([]int, r+1)
  mem[0] = 1
  prefix := make([]int, len(mem))
  for num, occ := range multiset {
    copy(prefix, mem)
    for sum := num; sum <= r; sum++ {
      prefix[sum] = (prefix[sum] + prefix[sum-num]) % mod
    }
    for sum := r; sum >= 0; sum-- {
      if num > 0 {
        mem[sum] = prefix[sum]
        if sum >= num*(occ+1) {
          mem[sum] = (mem[sum] - prefix[sum-num*(occ+1)] + mod) % mod
        }
      } else {
        mem[sum] = (mem[sum] * (occ + 1)) % mod
      }
    }
  }
  var result int
  for sum := l; sum <= r; sum++ {
    result = (result + mem[sum]) % mod
  }
  return result
}
var mod int = 1e9 + 7

function countSubMultisets(nums: number[], l: number, r: number): number {
  const cnt: number[] = Array(20001).fill(0);
  const memo: number[] = Array(20001).fill(0);
  const mod: number = 1000000007;
  for (const n of nums) {
    cnt[n]++;
  }
  memo.fill(1, 0, cnt[1] + 1);
  let total: number = cnt[1];
  for (let n = 2; n <= r; ++n) {
    if (!cnt[n]) {
      continue;
    }
    const top: number = (cnt[n] + 1) * n;
    total += n * cnt[n];
    for (let i = n, ii = Math.min(total, r); i <= ii; ++i) {
      memo[i] = (memo[i] + memo[i - n]) % mod;
    }
    for (let i = Math.min(total, r); i >= top; --i) {
      memo[i] = (mod + memo[i] - memo[i - top]) % mod;
    }
  }
  let result: number = 0;
  for (let i = l; i <= r; i++) {
    result = (result + memo[i]) % mod;
  }
  return (result * (cnt[0] + 1)) % mod;
}