結果
問題 | 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)