%I #8 Jan 04 2024 07:25:39
%S 0,1,1,1,1,1,3,1,1,1,1,1,1,1,1,1,1,3,1,3,1,1,1,5,1,1,1,1,1,1,3,1,1,1,
%T 1,1,1,1,3,1,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,3,1,1,1,1,1,
%U 5,1,1,1,1,3,1,1,1,1,1,1,1,1,1,1,1,3,1
%N The maximal exponent in the prime factorization of the exponentially odd numbers (A268335).
%C Differs from A368472 at n = 1, 154, 610, 707, 762, ... .
%H Amiram Eldar, <a href="/A368711/b368711.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A051903(A268335(n)).
%F a(n) is odd for n >= 2.
%F Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + 2 * Sum_{k>=1} (1 - Product_{p prime} (1 - 1/(p^(2*k-1)*(p^2+p-1)))) = 1.34877064483679975726... .
%t f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[AllTrue[e, OddQ], Max @@ e, Nothing]]; f[1] = 0; Array[f, 150]
%o (PARI) lista(kmax) = {my(e); print1(0, ", "); for(k = 2, kmax, e = factor(k)[,2]; if(vecprod(e)%2, print1(vecmax(e), ", ")));}
%Y Cf. A033150, A051903, A268335, A368472.
%Y Similar sequences: A368710, A368712, A368713.
%K nonn,easy
%O 1,7
%A _Amiram Eldar_, Jan 04 2024
|