結果

問題 No.2645 Sum of Divisors?
ユーザー chineristACchineristAC
提出日時 2024-02-19 22:04:43
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 91 ms / 2,000 ms
コード長 1,345 bytes
コンパイル時間 168 ms
コンパイル使用メモリ 81,700 KB
実行使用メモリ 80,604 KB
最終ジャッジ日時 2024-02-19 22:04:49
合計ジャッジ時間 3,835 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 57 ms
78,092 KB
testcase_01 AC 60 ms
78,092 KB
testcase_02 AC 78 ms
80,604 KB
testcase_03 AC 61 ms
78,092 KB
testcase_04 AC 59 ms
78,092 KB
testcase_05 AC 58 ms
78,092 KB
testcase_06 AC 88 ms
80,476 KB
testcase_07 AC 89 ms
80,476 KB
testcase_08 AC 58 ms
78,092 KB
testcase_09 AC 60 ms
78,092 KB
testcase_10 AC 57 ms
78,092 KB
testcase_11 AC 57 ms
78,092 KB
testcase_12 AC 81 ms
78,092 KB
testcase_13 AC 57 ms
78,092 KB
testcase_14 AC 60 ms
78,092 KB
testcase_15 AC 60 ms
78,092 KB
testcase_16 AC 60 ms
78,092 KB
testcase_17 AC 60 ms
78,092 KB
testcase_18 AC 61 ms
78,348 KB
testcase_19 AC 60 ms
78,348 KB
testcase_20 AC 60 ms
78,348 KB
testcase_21 AC 62 ms
78,348 KB
testcase_22 AC 63 ms
78,348 KB
testcase_23 AC 62 ms
80,476 KB
testcase_24 AC 64 ms
80,604 KB
testcase_25 AC 67 ms
80,604 KB
testcase_26 AC 67 ms
80,596 KB
testcase_27 AC 67 ms
80,604 KB
testcase_28 AC 89 ms
80,604 KB
testcase_29 AC 91 ms
80,604 KB
testcase_30 AC 73 ms
80,604 KB
testcase_31 AC 79 ms
80,592 KB
testcase_32 AC 86 ms
80,604 KB
testcase_33 AC 74 ms
80,604 KB
testcase_34 AC 84 ms
80,604 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
from itertools import permutations
from heapq import heappop,heappush
from collections import deque
import random
import bisect

input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

from math import log
from math import log1p

def isqrt(n):
    t = int(n**.5)
    for i in range(t-2,t+3)[::-1]:
        if i*i <= n:
            return i

"""
d(n) = prod (1+p+p^2+...+p^e)
d(n)/n = prod (1+p^-1+p^-2+...+p^-e)

ij <= N に対する 1/ijの総和
"""

M = 10**6
memo_cumulateive_inv_sum = [0] * (M+1)
for n in range(1,M+1):
    memo_cumulateive_inv_sum[n] = memo_cumulateive_inv_sum[n-1] + 1/n

def approximate(N):
    return log(N+0.5) + 0.57721566490152286060

def calc_cumulative_inv_sum(N):
    M = 10**6
    if N <= M:
        #print(N,memo_cumulateive_inv_sum[N]-approximate(N))
        return memo_cumulateive_inv_sum[N]
    return log(N+1/2) + 0.57721566490152286060

def calc_cumulative_inv_range_sum(l,r):
    return log1p((2*r-2*l+2)/(2*l-1))
    return log((2*r+1)/(2*l-1))

#M = 1000
def solve(N):
    B = isqrt(N)
    res = 0
    for n in range(1,B+1):
        res += 2*calc_cumulative_inv_sum(N//n)/n
    
    #print(res,calc_cumulative_inv_sum(B)**2)
    res -= calc_cumulative_inv_sum(B)**2
    
    
        
    return res

print(solve(int(input())))
0