結果

問題 No.2105 Avoid MeX
ユーザー hitonanodehitonanode
提出日時 2022-09-05 10:11:51
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 312 ms / 2,000 ms
コード長 1,173 bytes
コンパイル時間 7,813 ms
コンパイル使用メモリ 271,432 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-01 05:47:23
合計ジャッジ時間 10,380 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 6
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <vector>
using namespace std;

#include "testlib.h"

#include <atcoder/modint>
using mint = atcoder::modint998244353;

int main(int argc, char *argv[]) {
    registerValidation(argc, argv);
    int C = inf.readInt(0, 2000);
    inf.readSpace();
    int X = inf.readInt(0, 2000);
    inf.readEoln();
    inf.readEof();

    if (X == 0) {
        cout << mint(C + 1).inv().val() << '\n';
        return 0;
    }

    vector<mint> dp(X + 1);
    dp.at(0) = 1;
    for (int ub = C; ub < X; ++ub) {
        vector<mint> dpnxt(dp.size());
        for (int i = 0; i < int(dp.size()); ++i) {
            if (dp.at(i).val()) {
                mint c = i * mint(ub + 1).inv();
                dpnxt.at(i) += dp.at(i) * c;
                dpnxt.at(i + 1) += dp.at(i) * (1 - c);
            }
        }
        dp = dpnxt;
    }
    dp.pop_back();

    mint ret = 0;
    for (int i = 0; i < int(dp.size());  ++i) {
        // 0, ..., X - 1 のうち i 個がすでに埋まっている
        // 残りの X - i 個を埋めるより先に X を引けば良い
        ret += dp.at(i) * (1 - mint(X - i + 1).inv());
    }
    cout << ret.val() << endl;
}
0