#include "bits/stdc++.h" #pragma region cscpptemp // require at least C++17 to compile using llint = long long int; template using Heap = std::priority_queue; template using RevHeap = std::priority_queue, std::greater>; // ex. for (const int i : range(a, b)) ... loop [a, b) template >> class range { struct range_iterator { T itr; constexpr range_iterator(const int pos) noexcept : itr(pos) {} constexpr void operator++() noexcept { ++itr; } constexpr bool operator!=(const range_iterator& other) const noexcept { return itr != other.itr; } constexpr T operator*() const noexcept { return itr; } }; public: constexpr range(const T first_, const T last_) noexcept : first(first_), last(last_) {} constexpr range_iterator begin() const noexcept { return first; } constexpr range_iterator end() const noexcept { return last; } private: const range_iterator first, last; }; // ex. for (const int i : range(a, b)) ... revloop [a, b] template >> class revrange { struct revrange_iterator { T itr; constexpr revrange_iterator(const int pos) noexcept : itr(pos) {} constexpr void operator++() noexcept { --itr; } constexpr bool operator!=(const revrange_iterator& other) const noexcept { return itr != other.itr; } constexpr T operator*() const noexcept { return itr; } }; public: constexpr revrange(const T first_, const T last_) noexcept : first(first_), last(last_ - 1) {} constexpr revrange_iterator begin() const noexcept { return first; } constexpr revrange_iterator end() const noexcept { return last; } private: const revrange_iterator first, last; }; // ex. rec_lambda([](auto&& f, int n) -> int {return n < 1 ? 1 : f(--n); }) template class rec_lambda { public: constexpr rec_lambda(F&& f_) : f(std::forward(f_)) {} template constexpr auto operator()(Args&&... args) const { return f(*this, std::forward(args)...); } private: F f; }; // ex. const int N = read() struct read { template operator T() { T ret; std::cin >> ret; return ret; } }; // ex. chmax(a, b) ... a = max(a, b) template constexpr void chmax(T& a, const T b) noexcept { if (a < b) a = b; } template constexpr void chmin(T& a, const T b) noexcept { if (a > b) a = b; } // ex. ceildiv(a, b) ... ceil(a / b) template >> constexpr T ceildiv(const T dividend, const T divisor) { return (dividend + divisor - 1) / divisor; } // ex. mk_vec(a, b, c, v) ... 3D vector size: [a][b][c] init value: v template auto mk_vec(const int n, const T value) { return std::vector(n, value); } template auto mk_vec(const int n, Args... args) { return std::vector(n, mk_vec(args...)); } #pragma endregion int main() { std::ios::sync_with_stdio(false); std::cin.tie(nullptr); const int N = read(); const int M = read(); // M is prime // 0sum <= 1sum (all i) std::vector DP(N + 1, std::vector(N + 1, std::vector(N * N + 1, 0))); /* (0, 1) = (a, b) * invnum = k * DP[a + 1][b] */ DP[0][0][0] = 1; for (const int a : range(0, N + 1)) for (const int b : range(0, a + 1)) for (const int k : range(0, N * N + 1)) { if (a != N and k + b <= N * N) { DP[a + 1][b][k + b] += DP[a][b][k]; if (DP[a + 1][b][k + b] >= M) DP[a + 1][b][k + b] -= M; } if (a > b) { DP[a][b + 1][k] += DP[a][b][k]; if (DP[a][b + 1][k] >= M) DP[a][b + 1][k] -= M; } } for (const int i : range(0, N * N + 1)) std::cout << DP[N][N][i] << '\n'; // ............. return 0; }