結果
問題 | No.1173 Endangered Species |
ユーザー | anagohirame |
提出日時 | 2020-08-14 23:18:55 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 462 bytes |
コンパイル時間 | 88 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 74,576 KB |
最終ジャッジ日時 | 2024-10-10 16:44:41 |
合計ジャッジ時間 | 6,738 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | TLE | - |
testcase_01 | AC | 1,380 ms
67,992 KB |
testcase_02 | RE | - |
testcase_03 | AC | 1,379 ms
67,620 KB |
testcase_04 | RE | - |
testcase_05 | AC | 1,375 ms
68,004 KB |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
testcase_10 | RE | - |
testcase_11 | RE | - |
testcase_12 | RE | - |
testcase_13 | RE | - |
testcase_14 | RE | - |
testcase_15 | RE | - |
testcase_16 | TLE | - |
testcase_17 | -- | - |
testcase_18 | -- | - |
testcase_19 | -- | - |
testcase_20 | -- | - |
testcase_21 | -- | - |
ソースコード
from math import log from scipy import optimize n = int(input()) p = list(map(float, input().split())) q = list(map(float, input().split())) a = list(map(int, input().split())) def f(x): return sum((1-q[i]) / (1-q[i]*x) for i in range(n)) - x def f_prime(x): return sum((1/q[i] - 1) / ((1/q[i] - x)**2) for i in range(n)) - 1 r = optimize.newton(f, 0, fprime=f_prime) #print(r) ans = 0 for i in range(n): ans += a[i] * log((1-q[i]) / (1-q[i]*r)) print(ans)