結果
問題 | No.2105 Avoid MeX |
ユーザー |
![]() |
提出日時 | 2022-10-21 22:45:08 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 23 ms / 2,000 ms |
コード長 | 1,398 bytes |
コンパイル時間 | 2,143 ms |
コンパイル使用メモリ | 218,424 KB |
最終ジャッジ日時 | 2025-02-08 10:19:14 |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 6 |
other | AC * 19 |
ソースコード
#include<bits/stdc++.h>namespace {#pragma GCC diagnostic push#pragma GCC diagnostic ignored "-Wunused-function"#include<atcoder/all>#pragma GCC diagnostic popusing namespace std;using namespace atcoder;#define rep(i,n) for(int i = 0; i < (int)(n); i++)#define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--)#define all(x) begin(x), end(x)#define rall(x) rbegin(x), rend(x)template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; }template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; }using ll = long long;using P = pair<int,int>;using VI = vector<int>;using VVI = vector<VI>;using VL = vector<ll>;using VVL = vector<VL>;using mint = modint998244353;} int main() {ios::sync_with_stdio(false);cin.tie(0);int c, x;cin >> c >> x;if (x == 0) {cout << mint(c + 1).inv().val() << '\n';return 0;}vector<mint> dp(x + 1);dp[0] = 1;for(int i = 1; c + i <= x; i++) {int cnt = c + i;mint cinv = mint(cnt).inv();vector<mint> ndp(x + 1);rep(j, x) {mint v = dp[j];if (!v.val()) continue;ndp[j] += j * cinv * v;if (j + 1 <= x) {ndp[j+1] += (cnt - j) * cinv * v;}}swap(dp, ndp);}mint ans;rep(i, x + 1) ans += dp[i] / (x + 1 - i);ans = 1 - ans;cout << ans.val() << '\n';}