結果

問題 No.1173 Endangered Species
ユーザー anagohirameanagohirame
提出日時 2020-08-14 23:35:21
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 476 bytes
コンパイル時間 427 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 145,808 KB
最終ジャッジ日時 2024-10-10 16:52:02
合計ジャッジ時間 7,404 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 AC 1,381 ms
67,772 KB
testcase_02 AC 1,388 ms
67,904 KB
testcase_03 AC 1,405 ms
68,016 KB
testcase_04 AC 1,399 ms
68,032 KB
testcase_05 AC 1,408 ms
68,160 KB
testcase_06 AC 1,407 ms
68,040 KB
testcase_07 AC 1,391 ms
67,896 KB
testcase_08 AC 1,398 ms
68,020 KB
testcase_09 AC 1,387 ms
68,164 KB
testcase_10 AC 1,404 ms
68,036 KB
testcase_11 AC 1,384 ms
67,776 KB
testcase_12 AC 1,386 ms
68,016 KB
testcase_13 AC 1,418 ms
67,632 KB
testcase_14 AC 1,463 ms
68,508 KB
testcase_15 AC 1,446 ms
67,904 KB
testcase_16 AC 1,889 ms
75,752 KB
testcase_17 TLE -
testcase_18 TLE -
testcase_19 TLE -
testcase_20 TLE -
testcase_21 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

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(p[i] * (1-q[i]) / (1-q[i]*x) for i in range(n)) - x

def f_prime(x):
	return sum(p[i] * (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)
0