OFFSET
1,5
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[n_] := Module[{p = Times @@ (1 + FactorInteger[n][[;; , 2]]), e}, e = IntegerExponent[p, 2]; If[p == 2^e, e, Nothing]]; f[1] = 0; Array[f, 150]
PROG
(PARI) lista(kmax) = {my(p, e); for(k = 1, kmax, p = vecprod(apply(x -> x + 1, factor(k)[, 2])); e = valuation(p, 2); if(p >> e == 1, print1(e, ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, May 02 2024
STATUS
approved