結果
問題 | No.3004 ヤング図形 |
ユーザー |
![]() |
提出日時 | 2025-01-17 22:44:00 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
MLE
|
実行時間 | - |
コード長 | 1,361 bytes |
コンパイル時間 | 3,816 ms |
コンパイル使用メモリ | 276,448 KB |
実行使用メモリ | 784,900 KB |
最終ジャッジ日時 | 2025-01-17 22:45:06 |
合計ジャッジ時間 | 45,593 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | MLE * 3 |
other | MLE * 25 |
ソースコード
#include <bits/stdc++.h>using namespace std;using ll = long long;const int INF = 1e9 + 10;const ll INFL = 4e18;const int mod = 998244353;const int mx = 1e8 + 3;int fact[mx + 5], factinv[mx + 5];int modpow(ll a, int b) {ll ret = 1;while (b) {if (b & 1) (ret *= a) %= mod;(a *= a) %= mod;b /= 2;}return ret;}int main() {int K;cin >> K;vector<int> L(K), M(K);int N = 0;for (int i = 0; i < K; i++) {cin >> L[i] >> M[i];N += L[i] * M[i];}ll tmp;fact[0] = 1;for (int i = 0; i < mx; i++) {tmp = fact[i];tmp *= i + 1;tmp %= mod;fact[i + 1] = tmp;}factinv[mx] = modpow(fact[mx], mod - 2);for (int i = mx; i > 0; i--) {tmp = factinv[i];tmp *= i;tmp %= mod;factinv[i - 1] = tmp;}ll ans = fact[N];for (int i = 0; i < K; i++) {ans *= factinv[L[i] * M[i]];ans %= mod;}for (int i = 0; i < K; i++) {int base = L[i];for (int t = 0; t < M[i] - 1; t++) {int add = L[i] - 1;base += add;(ans *= fact[base]) %= mod;(ans *= factinv[add]) %= mod;(ans *= factinv[base - add]) %= mod;base++;}}cout << ans << endl;}