OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(p^e) = (1 + Sum_{k>=0} p^(2^k)), where e = Sum_{k>=0} e_k * 2^k, i.e., e_k is k-th bit of e.
MATHEMATICA
f[p_, e_] := 1 + Total[p^(2^(-1 + Position[Reverse@IntegerDigits[e, 2], 1] // Flatten))]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) a(n) = {my(f = factor(n), e); prod(k = 1, #f~, e = binary(f[k, 2]); 1 + sum(i = 1, #e, if(e[i], f[k, 1]^(2^(#e-i)), 0))); }
CROSSREFS
KEYWORD
nonn,mult,easy
AUTHOR
Amiram Eldar, Mar 14 2026
STATUS
approved
