結果
問題 |
No.2582 Random Average^K
|
ユーザー |
|
提出日時 | 2023-12-10 13:34:20 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,029 bytes |
コンパイル時間 | 259 ms |
コンパイル使用メモリ | 82,408 KB |
実行使用メモリ | 461,128 KB |
最終ジャッジ日時 | 2024-09-27 04:04:06 |
合計ジャッジ時間 | 5,388 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 7 TLE * 1 -- * 7 |
ソースコード
N,K = map(int,input().split()) mod = 998244353 class combination(): def __init__(self,N,p): self.fact = [1, 1] # fact[n] = (n! mod p) self.factinv = [1, 1] # factinv[n] = ((n!)^(-1) mod p) self.inv = [0, 1] # factinv 計算用 self.p = p for i in range(2, N + 1): self.fact.append((self.fact[-1] * i) % p) self.inv.append((-self.inv[p % i] * (p // i)) % p) self.factinv.append((self.factinv[-1] * self.inv[-1]) % p) def cmb(self,n, r): if (r < 0) or (n < r): return 0 r = min(r, n - r) return self.fact[n] * self.factinv[r] * self.factinv[n-r] % self.p C = combination(K+1+N,mod) ans = 0 for i in range(N): if i%2==0: ans = (ans + C.cmb(N,i)*pow(N-i,N+K,mod))%mod else: ans = (ans - C.cmb(N,i)*pow(N-i,N+K,mod))%mod inv = pow(C.fact[N+K]*C.factinv[K],mod-2,mod) ans = (ans*inv)%mod inv = pow(pow(N,K,mod),mod-2,mod) print(ans*inv%mod)