#pragma region Macros // #include #include using namespace std; // using namespace __gnu_pbds; // using namespace __gnu_cxx; // #include // using namespace atcoder; // #include // namespace mp = boost::multiprecision; // using Bint = mp::cpp_int; #define TO_STRING(var) # var #define pb emplace_back #define int ll #define endl '\n' using ll = long long; using ld = long double; const ld PI = acos(-1); const ld EPS = 1e-10; const ll INFL = 1LL << 61; // const int MOD = 998244353; const int MOD = 1000000007; __attribute__((constructor)) void constructor() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(15); } int modpow(int x, int N, int mod) { int ret = 1; while (N > 0) { if (N % 2 == 1) ret = ret * x % mod; x = x * x % mod; N /= 2; } return ret; } int ceil(int x, int y) { return (x > 0 ? (x + y - 1) / y : x / y); } #pragma endregion signed main() { int N, M; cin >> N >> M; vector> dp(N + 1,vector(N * N + 1)); dp[0][0] = 1; for (int i = 0; i < 2 * N; i++) { // 長さiの01列であって vector> ep(N + 1,vector(N * N + 1)); swap(dp, ep); for (int j = 0; j <= (i + 1) / 2; j++) { // 1の個数がj個で for (int k = 0; k < N * N; k++) { // 転倒数がk // 末尾に1を付け足しても転倒数は変わらない if (j != (i + 1) / 2) { dp[j + 1][k] += ep[j][k]; dp[j + 1][k] %= M; } if (k + j >= N * N) continue; // 末尾に0を付け足すと左にある1の個数分だけ転倒数が増える dp[j][k + j] += ep[j][k]; dp[j][k + j] %= M; } } } for (int i = 0; i <= N * N; i++) { cout << dp[N][i] << endl; } }