結果
問題 | No.2135 C5 |
ユーザー |
![]() |
提出日時 | 2022-11-25 23:29:58 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 158 ms / 2,000 ms |
コード長 | 1,378 bytes |
コンパイル時間 | 1,892 ms |
コンパイル使用メモリ | 172,940 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-10-02 05:59:39 |
合計ジャッジ時間 | 4,176 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 48 |
ソースコード
#include <bits/stdc++.h> using namespace std; const long long HALF = 499122177; const long long MOD = 998244353; int main(){ int N, M; cin >> N >> M; M = N * (N - 1) / 2 - M; if (M > N){ cout << 0 << endl; } else { vector<long long> fact(N + 1); fact[0] = 1; for (int i = 1; i <= N; i++){ fact[i] = fact[i - 1] * i % MOD; } vector<vector<long long>> binom(N + 1, vector<long long>(N + 1, 0)); for (int i = 0; i <= N; i++){ binom[i][0] = 1; for (int j = 1; j < i; j++){ binom[i][j] = (binom[i - 1][j - 1] + binom[i - 1][j]) % MOD; } binom[i][i] = 1; } vector<vector<long long>> dp(N + 1, vector<long long>(M + 1, 0)); dp[0][0] = 1; for (int i = 0; i < N; i++){ for (int j = 0; j <= M; j++){ for (int k = 5; k <= N - i; k++){ if (j + k <= M){ dp[i + k][j + k] += dp[i][j] * binom[N - i - 1][k - 1] % MOD * fact[k - 1] % MOD * HALF % MOD; dp[i + k][j + k] %= MOD; } } for (int k = 2; k <= N - i; k++){ if (j + k - 1 <= M){ dp[i + k][j + k - 1] += dp[i][j] * binom[N - i - 1][k - 1] % MOD * fact[k] % MOD * HALF % MOD; dp[i + k][j + k - 1] %= MOD; } } dp[i + 1][j] += dp[i][j]; dp[i + 1][j] %= MOD; } } cout << dp[N][M] << endl; } }