結果

問題 No.1967 Sugoroku Optimization
ユーザー asaringoasaringo
提出日時 2022-06-03 23:09:39
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 644 ms / 2,000 ms
コード長 1,275 bytes
コンパイル時間 1,877 ms
コンパイル使用メモリ 192,360 KB
最終ジャッジ日時 2025-01-29 17:59:13
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std ;
#define fast_input_output ios::sync_with_stdio(false); cin.tie(nullptr);
typedef long long ll ;
typedef long double ld ;
typedef pair<ll,ll> P ;
typedef tuple<ll,ll,ll> TP ;
#define chmin(a,b) a = min(a,b)
#define chmax(a,b) a = max(a,b)
#define bit_count(x) __builtin_popcountll(x)
#define gcd(a,b) __gcd(a,b)
#define lcm(a,b) a / gcd(a,b) * b
#define rep(i,n) for(int i = 0 ; i < n ; i++)
#define rrep(i,a,b) for(int i = a ; i < b ; i++)
#define endl "\n"

const int mod = 998244353 ;
ll powmod(ll x , ll n){
    ll res = 1 ;
    while(n > 0){
        if(n & 1) (res *= x) %= mod ;
        (x *= x) %= mod ;
        n >>= 1 ;
    }
    return res ;
}

int n , k ;
ll dp[2020][2020] ;

int main(){
    fast_input_output
    cin >> n >> k ;
    dp[0][0] = 1 ;
    rep(i,k){
        rep(j,n+1) {
            ll rem = n - j ;
            ll val = powmod(rem,mod-2) ;
            dp[i+1][j+1] += dp[i][j] * val % mod ;
            dp[i+1][n+1] -= dp[i][j] * val % mod ;
            dp[i+1][j+1] %= mod ;
            dp[i+1][n+1] += mod ;
            dp[i+1][n+1] %= mod ;
        }
        rep(j,n+2) (dp[i+1][j+1] += dp[i+1][j]) %= mod ;
    }
    ll res = 0 ;
    rep(i,k+1) (res += dp[i][n]) %= mod ;
    cout << res << endl ;
}
0