結果

問題 No.2528 pop_(backfront or not)
ユーザー Carpenters-Cat
提出日時 2023-11-03 22:03:10
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 430 ms / 2,000 ms
コード長 946 bytes
コンパイル時間 2,085 ms
コンパイル使用メモリ 192,728 KB
最終ジャッジ日時 2025-02-17 18:18:30
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
ll const m = 998244353;
int main () {
	int N;
	cin >> N;
	ll dp[4020][4020];
	ll cb[4020][4020];
	for (int n = 0; n <= N * 2 + 1; n ++) {
		cb[n][0] = cb[n][n] = 1;
		for (int r = 1; r < n; r ++) {
			cb[n][r] = (cb[n - 1][r] + cb[n - 1][r - 1]) % m;
		}
	}
	auto comb = [&](int n, int r) -> ll {
		if (r < 0 || n < r) return 0;
		return cb[n][r];
	};
	for (int l = 0; l <= N * 2; l ++) {
		for (int r = 0; r <= N * 2; r ++) {
			if (l + r == 0) {
				dp[l][r] = 1;
			} else if (l * r == 0) {
				dp[l][r] = 0;
			} else {
				dp[l][r] = (1ll + (l - 1) * (r - 1)) * dp[l - 1][r - 1];
				dp[l][r] %= m;
				if (l > 2) {
					dp[l][r] += comb(l - 1, 2) * dp[l - 2][r];
					dp[l][r] %= m;
				}
				if (r > 2) {
					dp[l][r] += comb(r - 1, 2) * dp[l][r - 2];
					dp[l][r] %= m;
				}
			}
		}
	}
	for (int i = 0; i < N * 2 + 1; i ++) {
		cout << dp[i][N * 2 - i] << endl;
	}
}
0