結果

問題 No.454 逆2乗和
ユーザー ukikagiukikagi
提出日時 2016-12-05 04:10:11
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 33 ms / 2,000 ms
コード長 519 bytes
コンパイル時間 175 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 10,752 KB
最終ジャッジ日時 2024-11-30 02:43:08
合計ジャッジ時間 2,275 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,624 KB
testcase_01 AC 31 ms
10,752 KB
testcase_02 AC 31 ms
10,624 KB
testcase_03 AC 31 ms
10,624 KB
testcase_04 AC 30 ms
10,752 KB
testcase_05 AC 31 ms
10,624 KB
testcase_06 AC 31 ms
10,752 KB
testcase_07 AC 31 ms
10,624 KB
testcase_08 AC 31 ms
10,624 KB
testcase_09 AC 31 ms
10,624 KB
testcase_10 AC 31 ms
10,624 KB
testcase_11 AC 31 ms
10,624 KB
testcase_12 AC 31 ms
10,624 KB
testcase_13 AC 31 ms
10,624 KB
testcase_14 AC 30 ms
10,752 KB
testcase_15 AC 31 ms
10,624 KB
testcase_16 AC 31 ms
10,752 KB
testcase_17 AC 30 ms
10,752 KB
testcase_18 AC 30 ms
10,624 KB
testcase_19 AC 31 ms
10,624 KB
testcase_20 AC 32 ms
10,624 KB
testcase_21 AC 31 ms
10,752 KB
testcase_22 AC 31 ms
10,752 KB
testcase_23 AC 32 ms
10,624 KB
testcase_24 AC 32 ms
10,624 KB
testcase_25 AC 33 ms
10,624 KB
testcase_26 AC 31 ms
10,752 KB
testcase_27 AC 30 ms
10,624 KB
testcase_28 AC 31 ms
10,624 KB
testcase_29 AC 31 ms
10,624 KB
testcase_30 AC 31 ms
10,624 KB
testcase_31 AC 31 ms
10,752 KB
testcase_32 AC 31 ms
10,752 KB
testcase_33 AC 31 ms
10,624 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from functools import lru_cache

@lru_cache(maxsize = None)
def comb(n, k):
  if k == 0 or k == n:
    return 1
  else:
    return comb(n-1, k-1) + comb(n-1, k)

@lru_cache(maxsize = None)
def B(n):
  if n == 0:
    return 1
  else:
    return -sum(comb(n+1, k) * B(k) for k in range(0, n)) / (n+1)

def solve(x, N = 100, K = 100):
  S  = sum((x+n)**(-2) for n in range(1, N+1))
  S += 1/(x+N) - 0.5/(x+N)**2
  S += sum(B(2*k)*(x+N)**(-2*k-1) for k in range(1,K+1))
  return S

x = float(input())
print(solve(x, 8, 4))
0