Partition(分割数) (math/fps/partition.hpp)

View this file on GitHub
Last update: 2024-04-24 02:34:04+09:00
Include: #include "math/fps/partition.hpp"

Verified with

test/verify/yosupo-partition-function.test.cpp

Code

/**
 * @brief Partition(分割数)
 */
template< template< typename > class FPS, typename Mint >
FPS< Mint > partition(int N) {
  FPS< Mint > poly(N + 1);
  poly[0] = 1;
  for(int k = 1; k <= N; k++) {
    if(1LL * k * (3 * k + 1) / 2 <= N) poly[k * (3 * k + 1) / 2] += (k % 2 ? -1 : 1);
    if(1LL * k * (3 * k - 1) / 2 <= N) poly[k * (3 * k - 1) / 2] += (k % 2 ? -1 : 1);
  }
  return poly.inv();
}

#line 1 "math/fps/partition.hpp"
/**
 * @brief Partition(分割数)
 */
template< template< typename > class FPS, typename Mint >
FPS< Mint > partition(int N) {
  FPS< Mint > poly(N + 1);
  poly[0] = 1;
  for(int k = 1; k <= N; k++) {
    if(1LL * k * (3 * k + 1) / 2 <= N) poly[k * (3 * k + 1) / 2] += (k % 2 ? -1 : 1);
    if(1LL * k * (3 * k - 1) / 2 <= N) poly[k * (3 * k - 1) / 2] += (k % 2 ? -1 : 1);
  }
  return poly.inv();
}