結果
問題 | No.2582 Random Average^K |
ユーザー | Navier_Boltzmann |
提出日時 | 2023-12-10 13:40:34 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,678 ms / 2,000 ms |
コード長 | 1,124 bytes |
コンパイル時間 | 272 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 75,776 KB |
最終ジャッジ日時 | 2025-01-17 02:07:14 |
合計ジャッジ時間 | 11,306 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 15 |
ソースコード
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 C = 1 for i in range(N): if i%2==0: ans = (ans + C*pow(N-i,N+K,mod))%mod else: ans = (ans - C*pow(N-i,N+K,mod))%mod C = (N-i)*pow(i+1,mod-2,mod)*C%mod inv = 1 for i in range(K+1,N+K+1): inv = (inv*i)%mod inv = pow(inv,mod-2,mod) ans = (ans*inv)%mod inv = pow(pow(N,K,mod),mod-2,mod) print(ans*inv%mod)