#925669 (C++17) No.2528 pop_(backfront or not)

• 問題: No.2528 pop_(backfront or not)
• ユーザー: Carpenters-Cat
• 提出日時: 2023-11-03 22:03:10
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 229 ms / 2,000 ms
• コード長: 946 bytes
• コンパイル時間: 2,747 ms
• コンパイル使用メモリ: 201,700 KB
• 実行使用メモリ: 206,720 KB
• 最終ジャッジ日時: 2024-09-25 20:22:27
• 合計ジャッジ時間: 4,995 ms
• ジャッジサーバーID （参考情報）: judge3 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	3 ms; 5,376 KB
testcase_04	AC	2 ms; 5,376 KB
testcase_05	AC	4 ms; 5,504 KB
testcase_06	AC	5 ms; 6,528 KB
testcase_07	AC	9 ms; 9,984 KB
testcase_08	AC	11 ms; 12,416 KB
testcase_09	AC	30 ms; 26,112 KB
testcase_10	AC	72 ms; 65,280 KB
testcase_11	AC	75 ms; 69,504 KB
testcase_12	AC	102 ms; 95,104 KB
testcase_13	AC	110 ms; 103,296 KB
testcase_14	AC	131 ms; 124,160 KB
testcase_15	AC	191 ms; 181,888 KB
testcase_16	AC	229 ms; 206,464 KB
testcase_17	AC	229 ms; 206,592 KB
testcase_18	AC	227 ms; 206,720 KB

ソースコード

#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;
	}
}