結果

問題 No.3072 Sum of sqrt(x)
ユーザー AtamaokaCAtamaokaC
提出日時 2021-10-24 23:02:23
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 629 bytes
コンパイル時間 513 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 89,216 KB
最終ジャッジ日時 2024-10-02 05:26:41
合計ジャッジ時間 12,242 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 132 ms
85,376 KB
testcase_01 AC 159 ms
80,256 KB
testcase_02 AC 130 ms
80,000 KB
testcase_03 AC 131 ms
80,256 KB
testcase_04 AC 129 ms
80,256 KB
testcase_05 AC 130 ms
80,128 KB
testcase_06 AC 133 ms
80,128 KB
testcase_07 TLE -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

def sqrt_floor(x):
    ceil = x
    floor = 1
    while ceil > floor + 1:
        m = (ceil+floor)//2
        if m*m <= x:
            floor = m
        else:
            ceil = m
    return floor

import decimal
from decimal import Decimal
decimal.getcontext().prec=20
inv_error = 10**20
N = int(input())
ans = Decimal(0)
for _ in range(N):
    x = int(input())
    p = sqrt_floor(x)
    if p*p == x:
        ans += p
    else:
        R = 1
        U, V = 2, 2*p
        a, b = 2*p, x-p*p
        while R < inv_error:
            R *= 2*p
            U, V = V, a*V + b*U
        ans += Decimal(V)/Decimal(U) - p
    print(ans)
0