結果

問題 No.2829 GCD Divination
ユーザー hato336hato336
提出日時 2024-08-02 21:54:32
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 226 ms / 2,000 ms
コード長 963 bytes
コンパイル時間 267 ms
コンパイル使用メモリ 82,156 KB
実行使用メモリ 90,368 KB
最終ジャッジ日時 2024-08-02 21:54:45
合計ジャッジ時間 6,300 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 122 ms
85,760 KB
testcase_01 AC 121 ms
85,632 KB
testcase_02 AC 135 ms
85,632 KB
testcase_03 AC 119 ms
85,376 KB
testcase_04 AC 146 ms
89,856 KB
testcase_05 AC 226 ms
90,112 KB
testcase_06 AC 213 ms
90,368 KB
testcase_07 AC 207 ms
90,240 KB
testcase_08 AC 213 ms
89,984 KB
testcase_09 AC 188 ms
90,112 KB
testcase_10 AC 126 ms
86,184 KB
testcase_11 AC 122 ms
85,760 KB
testcase_12 AC 134 ms
89,216 KB
testcase_13 AC 120 ms
85,760 KB
testcase_14 AC 121 ms
85,760 KB
testcase_15 AC 135 ms
89,328 KB
testcase_16 AC 120 ms
85,888 KB
testcase_17 AC 149 ms
89,088 KB
testcase_18 AC 123 ms
85,760 KB
testcase_19 AC 125 ms
85,760 KB
testcase_20 AC 122 ms
85,760 KB
testcase_21 AC 136 ms
89,344 KB
testcase_22 AC 123 ms
85,504 KB
testcase_23 AC 148 ms
89,984 KB
testcase_24 AC 130 ms
85,888 KB
testcase_25 AC 136 ms
85,760 KB
testcase_26 AC 123 ms
85,760 KB
testcase_27 AC 135 ms
89,156 KB
testcase_28 AC 125 ms
85,888 KB
testcase_29 AC 123 ms
85,632 KB
testcase_30 AC 127 ms
85,760 KB
testcase_31 AC 124 ms
85,632 KB
testcase_32 AC 130 ms
85,632 KB
testcase_33 AC 136 ms
89,216 KB
testcase_34 AC 145 ms
85,760 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import collections,sys,math,functools,operator,itertools,bisect,heapq,decimal,string,time,random
#sys.setrecursionlimit(10**9)
#sys.set_int_max_str_digits(0)
#input = sys.stdin.readline
def make_divisors(n):
    lower_divisors , upper_divisors = [], []
    i = 1
    while i*i <= n:
        if n % i == 0:
            lower_divisors.append(i)
            if i != n // i:
                upper_divisors.append(n//i)
        i += 1
    return lower_divisors + upper_divisors[::-1]
n = int(input())
@functools.lru_cache(maxsize=100000000)
def calc(x):
    if x == 1:
        return 0
    res = 0
    div = make_divisors(x)
    c = collections.defaultdict(int)
    for i in reversed(div):
        if i == x:
            c[x] = 1
            continue
        ans = x//i
        for j in div:
            if j % i == 0:
                ans -= c[j]
        c[i] = ans
        res += ans * calc(i)

    res /= x
    res += 1
    
    return res / (1 - 1/x)
print(calc(n))
0