結果
問題 | No.2141 Enumeratest |
ユーザー | 👑 rin204 |
提出日時 | 2022-12-02 21:35:10 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 23 ms / 2,000 ms |
コード長 | 943 bytes |
コンパイル時間 | 2,228 ms |
コンパイル使用メモリ | 200,940 KB |
実行使用メモリ | 19,200 KB |
最終ジャッジ日時 | 2024-10-09 22:41:19 |
合計ジャッジ時間 | 3,852 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 34 |
ソースコード
#include<bits/stdc++.h> using namespace std; using ll = long long; const ll MOD = 998244353; const int N = 1000100; ll fac[N], invfac[N]; ll modpow(ll a, ll b, ll MOD){ ll ret = 1; while (b > 0){ if(b & 1){ ret *= a; ret %= MOD; } a *= a; a %= MOD; b >>= 1; } return ret; } void init(){ fac[0] = 1; for(int i = 1; i < N; i++) fac[i] = fac[i - 1] * i % MOD; invfac[N - 1] = modpow(fac[N - 1], MOD - 2, MOD); for(int i = N - 2; i >= 0; i--) invfac[i] = invfac[i + 1] * (i + 1) % MOD; } ll nCk(int n, int k){ if(k < 0 || k > n) return 0; else return (fac[n] * invfac[k] % MOD) * invfac[n - k] % MOD; } void solve(){ init(); ll n, m; cin >> n >> m; ll ans = 1; int mi = m / n; int ma = mi + 1; int mac = m % n; int mic = n - m % n; while(mic--){ ans *= nCk(m, mi); ans %= MOD; m -= mi; } while(mac--){ ans *= nCk(m, ma); ans %= MOD; m -= ma; } cout << ans << endl; } int main(){ solve(); }