結果
問題 |
No.766 金魚すくい
|
ユーザー |
|
提出日時 | 2022-12-29 01:06:13 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 129 ms / 1,500 ms |
コード長 | 852 bytes |
コンパイル時間 | 402 ms |
コンパイル使用メモリ | 82,176 KB |
実行使用メモリ | 98,840 KB |
最終ジャッジ日時 | 2024-11-24 03:33:49 |
合計ジャッジ時間 | 5,698 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 42 |
ソースコード
MOD = 10 ** 9 + 7 n, m, p = map(int, input().split()) l = max(n + m, 100) fact = [1 for _ in range(l + 1)] inv = [1 for _ in range(l + 1)] ifact = [1 for _ in range(l + 1)] for x in range(2, l + 1): fact[x] = (fact[x - 1] * x) % MOD inv[x] = (- (MOD // x) * inv[MOD % x]) % MOD ifact[x] = (ifact[x - 1] * inv[x]) % MOD v = list(map(int, input().split())) v.sort(reverse=True) s = [0 for _ in range(n + 1)] for i in range(n): s[i + 1] = (s[i] + v[i]) % MOD pm = pow(p * inv[100], m, MOD) prob = [0 for _ in range(n + 2)] v = (100 - p) * inv[100] % MOD cur = 1 for i in range(1, n + 1): prob[i] = fact[i + m - 2] * ifact[i - 1] * ifact[m - 1] % MOD prob[i] = (prob[i] * pm * cur) % MOD cur = (cur * v) % MOD prob[n + 1] = (1 - sum(prob) % MOD) % MOD ans = 0 for i in range(1, n + 2): ans = (ans + prob[i] * s[i - 1]) % MOD print(ans)