結果
| 問題 |
No.2105 Avoid MeX
|
| コンテスト | |
| ユーザー |
Nachia
|
| 提出日時 | 2022-10-21 22:10:05 |
| 言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 27 ms / 2,000 ms |
| コード長 | 1,055 bytes |
| コンパイル時間 | 2,453 ms |
| コンパイル使用メモリ | 112,864 KB |
| 最終ジャッジ日時 | 2025-02-08 10:08:20 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 6 |
| other | AC * 19 |
ソースコード
#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <utility>
#include <atcoder/modint>
using namespace std;
using i32 = int32_t;
using u32 = uint32_t;
using i64 = int64_t;
using u64 = uint64_t;
#define rep(i,n) for(int i=0; i<(int)(n); i++)
const i64 INF = 1001001001001001001;
using Modint = atcoder::static_modint<998244353>;
int main(){
int C, X; cin >> C >> X;
if(X == 0){
Modint ans = Modint(C+1).inv();
cout << ans.val() << '\n';
return 0;
}
vector<Modint> dp0(X+1);
dp0[X] = 1;
for(int i=C+1; i<=X; i++){
Modint per = Modint(i).inv();
for(int j=0; j<X; j++){
auto mv = dp0[j+1] * (per * (i-(X-j-1)));
dp0[j] += mv;
dp0[j+1] -= mv;
}
}
Modint ans = 1;
rep(i,X+1) ans -= dp0[i] / (i+1);
cout << ans.val() << '\n';
return 0;
}
struct ios_do_not_sync{
ios_do_not_sync(){
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
}
} ios_do_not_sync_instance;
Nachia