login
a(n) = denominator of A065642(n)/n.
4

%I #12 May 09 2020 02:41:25

%S 1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,3,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1,

%T 1,3,1,1,1,4,1,1,1,1,3,1,1,8,1,5,1,1,1,3,1,4,1,1,1,2,1,1,3,1,1,1,1,1,

%U 1,1,1,3,1,1,5,1,1,1,1,4,1,1,1,2,1,1,1,1,1,3,1,1,1,1,1,8,1,7,1,5,1,1,1,1,1,1,1,3,1,1,1,4,1,1,1,1,1,1,1,4

%N a(n) = denominator of A065642(n)/n.

%H Antti Karttunen, <a href="/A285337/b285337.txt">Table of n, a(n) for n = 1..10000</a>

%o (Scheme) (define (A285337 n) (denominator (/ (A065642 n) n)))

%o (Python)

%o from sympy import primefactors, prod, Integer

%o def a007947(n): return 1 if n<2 else prod(primefactors(n))

%o def a065642(n):

%o if n==1: return 1

%o r=a007947(n)

%o n += r

%o while a007947(n)!=r:

%o n+=r

%o return n

%o def a(n): return (a065642(n)/Integer(n)).denominator() # _Indranil Ghosh_, Apr 20 2017

%Y Cf. A065642, A285100 (positions of ones), A284342 (positions of terms > 1).

%Y Cf. A285336 for the numerator.

%K nonn,frac

%O 1,12

%A _Antti Karttunen_, Apr 19 2017