結果

問題 No.454 逆2乗和
ユーザー ukikagiukikagi
提出日時 2016-12-05 03:50:12
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
RE  
実行時間 -
コード長 511 bytes
コンパイル時間 275 ms
コンパイル使用メモリ 12,416 KB
実行使用メモリ 13,656 KB
最終ジャッジ日時 2024-05-05 20:12:59
合計ジャッジ時間 2,389 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 53 ms
13,396 KB
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
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 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 AC 47 ms
13,528 KB
testcase_22 RE -
testcase_23 AC 48 ms
13,272 KB
testcase_24 AC 49 ms
13,276 KB
testcase_25 AC 51 ms
13,396 KB
testcase_26 AC 52 ms
13,400 KB
testcase_27 AC 51 ms
13,528 KB
testcase_28 AC 48 ms
13,524 KB
testcase_29 AC 51 ms
13,272 KB
testcase_30 AC 50 ms
13,396 KB
testcase_31 AC 47 ms
13,400 KB
testcase_32 AC 49 ms
13,272 KB
testcase_33 AC 48 ms
13,656 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 = int(input())
print(solve(x))
0