結果
問題 |
No.2582 Random Average^K
|
ユーザー |
![]() |
提出日時 | 2023-06-11 01:14:37 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 995 bytes |
コンパイル時間 | 7,536 ms |
コンパイル使用メモリ | 346,892 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-09-27 03:49:23 |
合計ジャッジ時間 | 9,551 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 5 RE * 10 |
ソースコード
#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("Ofast,unroll-loops") #include <bits/stdc++.h> #include <atcoder/all> #define rep(i, n) for (int i = 0; i < n; i++) #define ALL(a) a.begin(), a.end() #define ll long long using namespace std; using namespace atcoder; constexpr int mod = 998244353; using mint = modint998244353; mint fac[5050]; mint inv[5050]; void solve() { fac[0] = 1; for (int i = 1; i < 5050; i++) fac[i] = i * fac[i - 1]; inv[5049] = 1 / fac[5049]; for (int i = 5049; i > 0; i--) inv[i - 1] = i * inv[i]; int n, m; cin >> n >> m; mint ans = 0; for (int k = 1; k <= n; k++) { mint mul = fac[n] * inv[k] * inv[n - k]; if ((n - k) & 1) mul = -mul; mul *= mint(k).pow(n + m); ans += mul; } mint ff = 1; for (int i = 1; i <= n; i++) ff *= i + m; ans /= ff * mint(n).pow(m); cout << ans.val() << '\n'; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(13); solve(); return 0; }