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