結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 65 ms
78,488 KB
testcase_01 AC 65 ms
78,124 KB
testcase_02 AC 83 ms
80,572 KB
testcase_03 AC 63 ms
77,984 KB
testcase_04 AC 63 ms
78,344 KB
testcase_05 AC 67 ms
78,352 KB
testcase_06 AC 95 ms
79,204 KB
testcase_07 AC 96 ms
79,644 KB
testcase_08 AC 63 ms
77,700 KB
testcase_09 AC 64 ms
79,096 KB
testcase_10 AC 65 ms
78,108 KB
testcase_11 AC 65 ms
78,816 KB
testcase_12 AC 64 ms
77,868 KB
testcase_13 AC 63 ms
78,040 KB
testcase_14 AC 64 ms
78,172 KB
testcase_15 AC 63 ms
78,068 KB
testcase_16 AC 65 ms
78,940 KB
testcase_17 AC 63 ms
78,064 KB
testcase_18 AC 66 ms
78,988 KB
testcase_19 AC 65 ms
80,352 KB
testcase_20 AC 67 ms
79,072 KB
testcase_21 AC 66 ms
79,164 KB
testcase_22 AC 69 ms
79,472 KB
testcase_23 AC 67 ms
79,256 KB
testcase_24 AC 72 ms
80,752 KB
testcase_25 AC 74 ms
81,856 KB
testcase_26 AC 73 ms
80,908 KB
testcase_27 AC 72 ms
80,508 KB
testcase_28 AC 78 ms
80,748 KB
testcase_29 AC 96 ms
80,088 KB
testcase_30 AC 81 ms
80,960 KB
testcase_31 AC 85 ms
81,444 KB
testcase_32 AC 89 ms
80,280 KB
testcase_33 AC 79 ms
81,364 KB
testcase_34 AC 89 ms
81,416 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