結果
| 問題 | 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;
}
            
            
            
        