A372505 a(n) = log_2(A368473(n)).

0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 2, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 2, 2, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1

Offset: 1

Author

Amiram Eldar, May 04 2024

Keywords

• nonn
• easy

Comments

The first position of k, for k = 0, 1, ..., is 1, 4, 15, 126, 1134, ..., which is the position of A085629(2^k) in A138302.

Links

• Amiram Eldar, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A005361, A085629, A138302, A368473.

Cf. A369934, A370078, A372467, A372469.

Programs

Mathematica

f[n_] := Module[{p = Times @@ FactorInteger[n][[;; , 2]], e}, e = IntegerExponent[p, 2]; If[p == 2^e, e, Nothing]]; Array[f, 150]

PARI

lista(kmax) = {my(p, e); for(k = 1, kmax, p = vecprod(factor(k)[, 2]); e = valuation(p, 2); if(p >> e == 1, print1(e, ", ")));}

Formula

a(n) = log_2(A005361(A138302(n))).