%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