結果
| 問題 |
No.1554 array_and_me
|
| コンテスト | |
| ユーザー |
Kiri8128
|
| 提出日時 | 2021-06-18 23:07:43 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 366 ms / 2,000 ms |
| コード長 | 884 bytes |
| コンパイル時間 | 316 ms |
| コンパイル使用メモリ | 82,944 KB |
| 実行使用メモリ | 95,076 KB |
| 最終ジャッジ日時 | 2024-06-22 21:20:28 |
| 合計ジャッジ時間 | 11,938 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 41 |
ソースコード
from heapq import heappush, heappop
P = 998244353
nn = 202020
fa = [1] * (nn+1)
fainv = [1] * (nn+1)
for i in range(nn):
fa[i+1] = fa[i] * (i+1) % P
fainv[-1] = pow(fa[-1], P-2, P)
for i in range(nn)[::-1]:
fainv[i] = fainv[i+1] * (i+1) % P
C = lambda a, b: fa[a] * fainv[b] % P * fainv[a-b] % P if 0 <= b <= a else 0
T = int(input())
for _ in range(T):
N, K = map(int, input().split())
A = [int(a) for a in input().split()]
su = sum(A)
B = [max(K * a // su - 1, 0) for a in A]
su2 = sum(B)
H = []
for i in range(N):
heappush(H, (-A[i] / (B[i] + 1), i))
for k in range(K - su2):
a, i = heappop(H)
B[i] += 1
heappush(H, (-A[i] / (B[i] + 1), i))
ans = fa[K]
for i, b in enumerate(B):
ans = ans * fainv[b] % P * pow(A[i], b, P) % P
ans = ans * pow(su, (P - 2) * K, P) % P
print(ans)
Kiri8128