Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #14 Aug 28 2023 12:45:13
%S 0,2,1,3,2,4,3,2,3,5,3,3,3,2,1,4,2,5,1,3,2,6,2,4,2,1,3,5,3,6,1,2,3,2,
%T 3,10,4,4,5,5,8,7,7,2,4,7,3,2,4,3,2,5,3,4,2,8,3,4,4,5,3,3,7,2,5,10,4,
%U 2,6,8,3,6,6,4,3,6,4,7,4,4,3,4,8,5,7,4
%N Number of prime factors (counted with multiplicity) in factorization of A007408(n).
%H Amiram Eldar, <a href="/A124876/b124876.txt">Table of n, a(n) for n = 1..106</a>
%F a(n) = A001222(A007408(n)). - _R. J. Mathar_, May 18 2007
%e a(1) = 0 since A007408(1) = 1 contains no prime factor,
%e a(2) = 2 since A007408(2) = 9 = 3 * 3,
%e a(3) = 1 since A007408(3) = 251 is prime,
%e a(6) = 4 since A007408(6) = 7 * 7 * 11 * 53.
%p seq( add(op(2,j),j=op(2,(ifactors@A007408)(n))), n=1..28 );
%p A001222 := proc(n) numtheory[bigomega](n) ; end: b := fscanf("b007408.txt","%d %d") : while b <> [] do printf("%d, ",A001222(op(2,b))) ; b := fscanf("b007408.txt","%d %d") : od : # _R. J. Mathar_, May 18 2007
%t Table[PrimeOmega[Numerator[Sum[1/k^3, {k, 1, n}]]], {n, 1, 50}] (* _Amiram Eldar_, Feb 09 2020 *)
%t PrimeOmega[Numerator[Accumulate[1/Range[50]^3]]] (* _Harvey P. Dale_, Aug 28 2023 *)
%Y Cf. A001222, A007408, A124877, A124787.
%K nonn
%O 1,2
%A _M. F. Hasler_, Nov 11 2006
%E More terms from _R. J. Mathar_, May 18 2007
%E More terms from _Amiram Eldar_, Feb 09 2020