OFFSET
1,6
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
s[0] = 0; s[n_] := s[n] = s[Floor[n/4]] + If[OddQ[Mod[n, 4]], 1, 0]; f[p_, e_] := s[e]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) s(e) = if(e>3, s(e\4)) + e%2 \\ after Charles R Greathouse IV at A139351
a(n) = vecsum(apply(s, factor(n)[, 2]));
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Oct 05 2023
STATUS
approved