結果
問題 | No.1554 array_and_me |
ユーザー |
![]() |
提出日時 | 2021-06-18 23:33:27 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 955 bytes |
コンパイル時間 | 183 ms |
コンパイル使用メモリ | 81,972 KB |
実行使用メモリ | 124,244 KB |
最終ジャッジ日時 | 2024-06-22 21:33:20 |
合計ジャッジ時間 | 32,915 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 19 WA * 22 |
ソースコード
from fractions import Fractionimport heapqimport sysinput = sys.stdin.buffer.readlinesys.setrecursionlimit(10 ** 7)mod = 998244353U = 10 ** 5 + 10fact = [1] * (U+1)inv = [1] * (U+1)for i in range(1, U+1):fact[i] = fact[i-1] * i % modinv[U] = pow(fact[i], mod-2, mod)for i in reversed(range(U)):inv[i] = inv[i+1] * (i + 1) % modT = int(input())for _ in range(T):N, K = map(int, input().split())A = list(map(int, input().split()))S = sum(A) % modS_inv = pow(S, mod-2, mod)que = [(-Fraction(a), i) for i, a in enumerate(A)]heapq.heapify(que)cnt = [0] * Nfor _ in range(K):f, idx = heapq.heappop(que)f = -fnum, div = f.as_integer_ratio()cnt[idx] += 1heapq.heappush(que, (-Fraction(num, div + 1), idx))res = fact[K]for i in range(N):res = res * pow(A[i] * S_inv, cnt[i], mod)res = res * inv[cnt[i]] % modprint(res)