OFFSET
1,2
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
PROG
(Python)
from sympy import primefactors, prod, Integer
def a007947(n): return 1 if n<2 else prod(primefactors(n))
def a065642(n):
if n==1: return 1
r=a007947(n)
n += r
while a007947(n)!=r:
n += r
return n
def a(n): return (a065642(n)/Integer(n)).numerator() # Indranil Ghosh, Apr 20 2017
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Antti Karttunen, Apr 19 2017
STATUS
approved