結果

問題 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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

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