結果

問題 No.3004 ヤング図形
ユーザー Today03Today03
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0