結果

問題 No.454 逆2乗和
ユーザー はむ吉🐹はむ吉🐹
提出日時 2017-01-01 21:50:31
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 664 bytes
コンパイル時間 154 ms
コンパイル使用メモリ 10,896 KB
実行使用メモリ 9,896 KB
最終ジャッジ日時 2023-08-22 09:21:40
合計ジャッジ時間 3,143 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 29 ms
9,692 KB
testcase_01 AC 29 ms
9,708 KB
testcase_02 AC 29 ms
9,708 KB
testcase_03 AC 29 ms
9,704 KB
testcase_04 AC 29 ms
9,756 KB
testcase_05 AC 29 ms
9,736 KB
testcase_06 AC 30 ms
9,676 KB
testcase_07 AC 30 ms
9,788 KB
testcase_08 AC 29 ms
9,832 KB
testcase_09 AC 28 ms
9,732 KB
testcase_10 AC 29 ms
9,876 KB
testcase_11 AC 29 ms
9,676 KB
testcase_12 AC 29 ms
9,796 KB
testcase_13 AC 30 ms
9,872 KB
testcase_14 AC 29 ms
9,676 KB
testcase_15 AC 29 ms
9,704 KB
testcase_16 AC 29 ms
9,700 KB
testcase_17 AC 29 ms
9,672 KB
testcase_18 AC 29 ms
9,688 KB
testcase_19 AC 29 ms
9,708 KB
testcase_20 AC 29 ms
9,896 KB
testcase_21 AC 29 ms
9,592 KB
testcase_22 WA -
testcase_23 AC 29 ms
9,764 KB
testcase_24 AC 30 ms
9,608 KB
testcase_25 AC 29 ms
9,596 KB
testcase_26 AC 30 ms
9,752 KB
testcase_27 AC 29 ms
9,708 KB
testcase_28 AC 29 ms
9,604 KB
testcase_29 AC 29 ms
9,704 KB
testcase_30 AC 30 ms
9,716 KB
testcase_31 AC 29 ms
9,596 KB
testcase_32 AC 29 ms
9,800 KB
testcase_33 AC 28 ms
9,836 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#!/usr/bin/env python3

import fractions


def inv(x):
    return fractions.Fraction(1, x)


# http://www.cis.nagasaki-u.ac.jp/~masada/Digamma2010042801.pdf
def trigamma(x):
    ans = inv(x)
    ans += inv(2 * x ** 2)
    ans += inv(6 * x ** 3)
    ans -= inv(30 * x ** 5)
    ans += inv(43 * x ** 7)
    ans -= inv(30 * x ** 9)
    ans += 5 * inv(66 * x ** 11)
    ans -= 691 * inv(2730 * x ** 13)
    ans += 7 * inv(6 * x ** 15)
    return ans


def main():
    x = fractions.Fraction(input())
    if x == 0:
        res = 1.6449340668482426
    else:
        res = trigamma(x + 1)
    print("{:.16f}".format(float(res)))


if __name__ == '__main__':
    main()
0