cp's OEIS Frontend

A286582 a(n) = A001222(A048673(n)).

%I A286582 #16 May 15 2020 16:08:28
%S A286582 0,1,1,1,2,3,2,2,1,1,1,1,2,1,3,1,2,2,3,5,3,3,2,3,2,2,3,3,4,1,1,2,2,1,
%T A286582 2,1,2,2,1,2,2,1,4,1,4,3,3,2,1,2,5,2,3,3,2,1,2,1,1,2,2,4,3,2,4,3,4,2,
%U A286582 1,3,1,3,4,2,2,4,5,7,3,3,1,2,3,4,1,1,3,5,2,1,2,1,2,5,4,6,2,3,1,2,3,2,4,3,1,1
%N A286582 a(n) = A001222(A048673(n)).
%H A286582 Antti Karttunen, <a href="/A286582/b286582.txt">Table of n, a(n) for n = 1..10000</a>
%F A286582 a(n) = A001222(A048673(n)).
%o A286582 (PARI)
%o A286582 A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Using code of _Michel Marcus_
%o A286582 A048673(n) = (A003961(n)+1)/2;
%o A286582 A286582(n) = bigomega(A048673(n));
%o A286582 (Scheme) (define (A286582 n) (A001222 (A048673 n)))
%o A286582 (Python)
%o A286582 from sympy import factorint, nextprime, primefactors, prod
%o A286582 def a001222(n): return 0 if n==1 else a001222(n//primefactors(n)[-1]) + 1
%o A286582 def a048673(n):
%o A286582     f = factorint(n)
%o A286582     return 1 if n==1 else (1 + prod(nextprime(i)**f[i] for i in f))//2
%o A286582 def a(n): return a001222(a048673(n))
%o A286582 print([a(n) for n in range(1, 51)]) # _Indranil Ghosh_, Jun 12 2017
%Y A286582 Cf. A001222, A048673, A091304, A278224, A286583, A286584, A286585.
%K A286582 nonn
%O A286582 1,5
%A A286582 _Antti Karttunen_, May 31 2017