結果
問題 |
No.3153 probability max K
|
ユーザー |
|
提出日時 | 2025-05-21 21:01:19 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 162 ms / 2,000 ms |
コード長 | 839 bytes |
コンパイル時間 | 1,957 ms |
コンパイル使用メモリ | 203,380 KB |
実行使用メモリ | 15,752 KB |
最終ジャッジ日時 | 2025-05-21 21:01:27 |
合計ジャッジ時間 | 6,762 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include<bits/stdc++.h> #include<atcoder/modint> #define int long long using namespace std; using mint=atcoder::modint998244353; #define double long double signed main(){ int n,k;cin>>n>>k; vector<int>a(n); for(int i=0;i<n;i++)cin>>a[i]; vector<vector<mint>>dp(n+1,vector<mint>(2,0)); dp[0][0]=1; for(int i=0;i<n;i++){ if(a[i]>=k){ dp[i+1][0]+=dp[i][0]*((mint)(k-1)/a[i]); dp[i+1][1]+=dp[i][0]*(mint(1)/a[i]); dp[i+1][1]+=dp[i][1]*((mint)k/a[i]); // dp[i+1][0]+=dp[i][1]*((a[i]-k)/a[i]); }else{ dp[i+1][1]+=dp[i][1]; dp[i+1][0]+=dp[i][0]; } } // for(int i=0;i<2;i++){ // for(int j=0;j<=n;j++)cout<<fixed<<setprecision(18)<<dp[j][i]<<" "; // cout<<endl; // } cout<<dp[n][1].val()<<endl; }