OFFSET
1,35
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := DivisorSigma[0, e]; d[1] = 1;
d[n_] := Times @@ f @@@ FactorInteger[n];
pow2Q[n_] := n == 2^IntegerExponent[n, 2];
q[n_] := pow2Q[d[n]]; Log2[d /@ Select[Range[100], q]]
PROG
(PARI) d(n) = vecprod(apply(numdiv, factor(n)[, 2])); \\ A049419
ispow2(n) = n >> valuation(n, 2) == 1;
list(kmax) = {my(d1); for(k = 1, kmax, d1 = d(k); if(ispow2(d1), print1(valuation(d1, 2), ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Sep 22 2025
STATUS
approved
