結果
| 問題 |
No.2105 Avoid MeX
|
| コンテスト | |
| ユーザー |
risujiroh
|
| 提出日時 | 2022-10-21 23:47:09 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 21 ms / 2,000 ms |
| コード長 | 867 bytes |
| コンパイル時間 | 3,136 ms |
| コンパイル使用メモリ | 250,980 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-07-01 07:52:44 |
| 合計ジャッジ時間 | 4,126 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 6 |
| other | AC * 19 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/modint>
using Fp = atcoder::modint998244353;
void solve(std::istream& is, std::ostream& os) {
int C, X;
is >> C >> X;
if (X == 0) {
os << Fp(C + 1).inv().val() << '\n';
} else if (X <= C) {
os << (1 - Fp(X + 1).inv()).val() << '\n';
} else {
std::vector<Fp> f{1};
for (int i = C + 1; i <= X; ++i) {
std::vector<Fp> nf(f.size() + 1);
Fp r = Fp(i).inv();
for (int j = 0; j < int(f.size()); ++j) {
nf[j] += f[j] * j * r;
nf[j + 1] += f[j] * (i - j) * r;
}
f = std::move(nf);
}
Fp ans = 0;
for (int i = 1; i < int(f.size()); ++i) {
ans += f[i] * Fp(X - i + 1).inv();
}
ans = 1 - ans;
os << ans.val() << '\n';
}
}
int main() {
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
solve(std::cin, std::cout);
}
risujiroh