OFFSET
1,6
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
For n = 8 = 2^3, the representation of 3 in factorial base is 11, i.e., 3 = 1! + 2!, so 8 = (2^(1!))^1 * (2^(2!))^1 and a(8) = 1 + 1 = 2.
For n = 16 = 2^4, the representation of 4 in factorial base is 20, i.e., 4 = 2 * 2!, so 16 = (2^(2!))^2 and a(16) = 1.
MATHEMATICA
fdignum[n_] := Module[{k = n, m = 2, r, s = 0}, While[{k, r} = QuotientRemainder[k, m]; k != 0 || r != 0, If[r > 0, s++]; m++]; s]; f[p_, e_] := fdignum[e]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) fdignum(n) = {my(k = n, m = 2, r, s = 0); while([k, r] = divrem(k, m); k != 0 || r != 0, if(r > 0, s ++); m++); s; }
a(n) = {my(e = factor(n)[, 2]); sum(i = 1, #e, fdignum(e[i])); }
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Amiram Eldar, Oct 08 2024
STATUS
approved