結果
| 問題 |
No.2829 GCD Divination
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-08-02 22:16:09 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,044 bytes |
| コンパイル時間 | 362 ms |
| コンパイル使用メモリ | 82,228 KB |
| 実行使用メモリ | 84,896 KB |
| 最終ジャッジ日時 | 2024-08-02 22:16:16 |
| 合計ジャッジ時間 | 6,480 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 4 TLE * 2 -- * 29 |
ソースコード
from functools import lru_cache
import sys, time, random
from collections import deque, Counter, defaultdict
input = lambda: sys.stdin.readline().rstrip()
ii = lambda: int(input())
mi = lambda: map(int, input().split())
li = lambda: list(mi())
inf = 2 ** 61 - 1
mod = 998244353
from math import gcd
n = ii()
ni = 1 / n
@lru_cache(maxsize=None)
def solve(x):
if x == 1:
return 0
ans = 0
nx = x
p = []
for i in range(2, x):
if i * i > nx:
break
if nx % i == 0:
for j in range(i, x, i):
for k in p:
if j % k == 0:
break
else:
ans += solve(gcd(x, j))
while nx % i == 0:
nx //= i
p.append(i)
if nx > 1:
for j in range(nx, x, nx):
for k in p:
if j % k == 0:
break
else:
ans += solve(gcd(x, j))
ans += x
return ans / (x - 1)
print(solve(n))