Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A285336
a(n) = numerator of A065642(n)/n.
2
1, 2, 3, 2, 5, 2, 7, 2, 3, 2, 11, 3, 13, 2, 3, 2, 17, 4, 19, 2, 3, 2, 23, 3, 5, 2, 3, 2, 29, 2, 31, 2, 3, 2, 5, 4, 37, 2, 3, 5, 41, 2, 43, 2, 5, 2, 47, 9, 7, 8, 3, 2, 53, 4, 5, 7, 3, 2, 59, 3, 61, 2, 7, 2, 5, 2, 67, 2, 3, 2, 71, 4, 73, 2, 9, 2, 7, 2, 79, 5, 3, 2, 83, 3, 5, 2, 3, 2, 89, 4, 7, 2, 3, 2, 5, 9, 97, 8, 3, 8, 101, 2, 103, 2, 3, 2, 107, 4, 109, 2, 3
OFFSET
1,2
LINKS
PROG
(Scheme) (define (A285336 n) (numerator (/ (A065642 n) n)))
(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
Cf. A065642.
Cf. A285337 for the denominator.
Sequence in context: A071889 A091963 A067695 * A273282 A079866 A134332
KEYWORD
nonn,frac
AUTHOR
Antti Karttunen, Apr 19 2017
STATUS
approved