結果

問題 No.2530 Yellow Cards
ユーザー syndro_6syndro_6
提出日時 2023-11-03 22:35:59
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 318 ms / 2,000 ms
コード長 4,452 bytes
コンパイル時間 1,946 ms
コンパイル使用メモリ 208,036 KB
実行使用メモリ 198,912 KB
最終ジャッジ日時 2024-09-25 20:59:04
合計ジャッジ時間 5,675 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 3 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 307 ms
198,784 KB
testcase_08 AC 310 ms
198,784 KB
testcase_09 AC 307 ms
198,912 KB
testcase_10 AC 312 ms
198,784 KB
testcase_11 AC 318 ms
198,272 KB
testcase_12 AC 313 ms
198,528 KB
testcase_13 AC 310 ms
198,272 KB
testcase_14 AC 219 ms
140,032 KB
testcase_15 AC 150 ms
90,240 KB
testcase_16 AC 144 ms
92,544 KB
testcase_17 AC 238 ms
152,320 KB
testcase_18 AC 89 ms
60,928 KB
testcase_19 AC 23 ms
16,768 KB
testcase_20 AC 35 ms
23,424 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// (⁠◕⁠ᴗ⁠◕⁠✿⁠)
 
#include <bits/stdc++.h>
#define rep(i, n) for (ll i = 0; i < (n); i++)
#define srep(i, s, n) for (ll i = s; i < (n); i++)
#define len(x) ((int)(x).size())
#define all(x) (x).begin(), (x).end()
using namespace std;
template<typename T> using vc = vector<T>;
template<typename T> using vv = vc<vc<T>>;
template<typename T> using vvv = vv<vc<T>>;
using vi = vc<int>;using vvi = vv<int>; using vvvi = vv<vi>;
using ll = long long;using vl = vc<ll>;using vvl = vv<ll>; using vvvl = vv<vl>;
using ld = long double; using vld = vc<ld>; using vvld = vc<vld>; using vvvld = vc<vvld>;
using uint = unsigned int;
using ull = unsigned long long;
const ld pi = 3.141592653589793;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
// const ll mod = 1000000007;
const ll mod = 998244353;
inline bool inside(long long y, long long x, long long H, long long W) {return 0 <= (y) and (y) < (H) and 0 <= (x) and (x) < (W); }

#define debug(var)  do{std::cout << #var << " : \n";view(var);}while(0)
template<typename T> void view(T e){cout << e << endl;}
template<typename T> void view(const vc<T>& v){for(const auto& e : v){ cout << e << " "; } cout << endl;}
template<typename T> void view(const vv<T>& vv){ for(const auto& v : vv){ view(v); } }

struct modint {
    ll val;
    using mint = modint;
    modint() {};
    modint(ll v) {
        if (0 <= v && v < mod) val = v;
        else{
            v %= mod;
            if (v < 0) v += mod;
            val = v;
        }
    }
    mint operator++(int){
        val++;
        if (val == mod) val = 0;
        return *this;
    }
    mint operator--(int){
        if (val == 0) val = mod;
        val--;
        return *this;
    }
    mint& operator+=(const mint& v){
        val += v.val;
        if (val > mod) val -= mod;
        return *this;
    }
    friend mint operator+(const mint& u, const mint& v){
        return mint(u) += v;
    }
    mint& operator-=(const mint& v){
        val -= v.val;
        if (val < 0) val += mod;
        return *this;
    }
    friend mint operator-(const mint& u, const mint& v){
        return mint(u) -= v;
    }
    mint& operator*=(const mint& v){
        val *= v.val;
        if (val > mod) val %= mod;
        return *this;
    }
    friend mint operator*(const mint& u, const mint& v){
        return mint(u) *= v;
    }
    mint pow(ll n){
        mint a = *this, ret = 1;
        while (n > 0){
            if (n & 1) ret *= a;
            a *= a;
            n >>= 1;
        }
        return ret;
    }
    mint operator/=(const mint& v){
        val *= mint(v).pow(mod - 2).val;
        if (val > mod) val %= mod;
        return *this;
    }
    friend mint operator/(const mint& u, const mint& v){
        return mint(u) /= v;
    }
    friend bool operator==(const mint& u, const mint& v){
        return u.val == v.val;
    }
    friend bool operator!=(const mint& u, const mint& v){
        return u.val != v.val;
    }
    friend bool operator>(const mint& u, const mint& v){
        return u.val > v.val;
    }
    friend bool operator<(const mint& u, const mint& v){
        return u.val < v.val;
    }
};

using mint = modint;

struct Binoms {
    public : 
        Binoms(int n){
            fact = vc<mint>(1, 1);
            inv_fact = vc<mint>(1, 1);
            update(n);
        };

        void update(int n){
            if (sz >= n) return;
            fact.resize(n + 1, 1);
            inv_fact.resize(n + 1, 1);
            for (int i = sz; i < n; i++) fact[i + 1] = fact[i] * (i + 1);
            inv_fact[n] = 1 / fact[n];
            for (int i = n - 1; i >= sz + 1; i--) inv_fact[i] = inv_fact[i + 1] * (i + 1);
            sz = n;
            return;
        }

        mint C(int n, int k){
            if (sz < n) update(n);
            return fact[n] * inv_fact[n - k] * inv_fact[k];
        }

        mint fac(int n){
            if (sz < n) update(n);
            return fact[n];
        }

    private :
        int sz = 0;
        vc<mint> fact, inv_fact;
};

int main(){
    int N, K; cin >> N >> K;
    vv<mint> dp(K + 1, vc<mint>(N + 1, 0));
    mint inv_N = mint(N).pow(mod - 2);
    dp[0][N] = 1;
    rep(i, K) rep(j, N + 1){
        if (j > 0) dp[i + 1][j - 1] += dp[i][j] * j * inv_N;
        if (j < N) dp[i + 1][j + 1] += dp[i][j] * (N - j) * inv_N;
    }
    mint ret = N;
    rep(j, N + 1) ret += dp[K][j] * max(0ll, (K - (N - j)) / 2);
    cout << ret.val << endl;
}
0