OFFSET
1,8
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
f[p_, e_] := Floor[e/2] + If[OddQ[e] || OddQ[DigitCount[e + 1, 2, 1]], 1, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) s(n) = n\2 + (n%2 || hammingweight(n+1)%2); \\ after Charles R Greathouse IV at A159481
a(n) = vecprod(apply(x -> s(x), factor(n)[, 2]));
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Amiram Eldar, Oct 27 2023
STATUS
approved