結果
問題 | No.2105 Avoid MeX |
ユーザー | 沙耶花 |
提出日時 | 2022-10-21 22:56:35 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 244 ms / 2,000 ms |
コード長 | 1,755 bytes |
コンパイル時間 | 4,581 ms |
コンパイル使用メモリ | 269,468 KB |
実行使用メモリ | 4,700 KB |
最終ジャッジ日時 | 2023-09-13 23:07:51 |
合計ジャッジ時間 | 7,000 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 3 ms
4,380 KB |
testcase_01 | AC | 4 ms
4,472 KB |
testcase_02 | AC | 3 ms
4,376 KB |
testcase_03 | AC | 3 ms
4,380 KB |
testcase_04 | AC | 119 ms
4,424 KB |
testcase_05 | AC | 3 ms
4,376 KB |
testcase_06 | AC | 51 ms
4,384 KB |
testcase_07 | AC | 86 ms
4,444 KB |
testcase_08 | AC | 24 ms
4,380 KB |
testcase_09 | AC | 141 ms
4,404 KB |
testcase_10 | AC | 3 ms
4,476 KB |
testcase_11 | AC | 4 ms
4,384 KB |
testcase_12 | AC | 25 ms
4,508 KB |
testcase_13 | AC | 3 ms
4,380 KB |
testcase_14 | AC | 68 ms
4,560 KB |
testcase_15 | AC | 89 ms
4,400 KB |
testcase_16 | AC | 4 ms
4,380 KB |
testcase_17 | AC | 4 ms
4,380 KB |
testcase_18 | AC | 4 ms
4,700 KB |
testcase_19 | AC | 237 ms
4,408 KB |
testcase_20 | AC | 240 ms
4,692 KB |
testcase_21 | AC | 244 ms
4,576 KB |
testcase_22 | AC | 244 ms
4,400 KB |
testcase_23 | AC | 4 ms
4,380 KB |
testcase_24 | AC | 4 ms
4,376 KB |
ソースコード
#include <stdio.h> #include <bits/stdc++.h> #include <atcoder/all> using namespace atcoder; using mint = modint998244353; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf32 1000000001 #define Inf64 3000000000000000 struct combi{ deque<mint> kaijou; deque<mint> kaijou_; combi(int n){ kaijou.push_back(1); for(int i=1;i<=n;i++){ kaijou.push_back(kaijou[i-1]*i); } mint b=kaijou[n].inv(); kaijou_.push_front(b); for(int i=1;i<=n;i++){ int k=n+1-i; kaijou_.push_front(kaijou_[0]*k); } } mint combination(int n,int r){ if(r>n)return 0; mint a = kaijou[n]*kaijou_[r]; a *= kaijou_[n-r]; return a; } mint junretsu(int a,int b){ mint x = kaijou_[a]*kaijou_[b]; x *= kaijou[a+b]; return x; } mint catalan(int n){ return combination(2*n,n)/(n+1); } }; int main(){ combi c(100000); int C,X; cin>>C>>X; vector dp(X+1,vector<mint>(2,0)); dp[0][0] =1; for(int i=1;i<=100000;i++){ int l = 0; int r = C+i-1; mint len = r-l+1; len = len.inv(); if(r>X)break; vector ndp(X+1,vector<mint>(2,0)); rep(j,X+1){ { mint p = j; p *= len; ndp[j][0] += dp[j][0] * p; } if(j!=X){ mint p = min(r-l+1,X)-j; p *= len; ndp[j+1][0] += dp[j][0] * p; } if(r==X){ mint p = 1; p *= len; ndp[j][1] += dp[j][0]*p; } } swap(dp,ndp); } mint ans = 0; rep(i,X){ rep(j,2){ //cout<<(dp[i][j]*6).val()<<','; if(j==0){ int r = X-i; mint p = c.kaijou[r]; p *= c.kaijou_[r+1]; ans += dp[i][j] * (1-p); } else{ ans += dp[i][j]; } } //cout<<endl; } if(X==0){ int l = 0; int r = C+1-1; ans = 1; ans /= r-l+1; //ans = 1-ans; } cout<<ans.val()<<endl; return 0; }