OFFSET
1,6
COMMENTS
Let n = Product p^e be the canonical prime factorization of n. The factorization of n that is based on the factorial-base representation of the exponents is done by factoring each prime power p^e into prime powers, p^e = Product_{k} (p^(k!))^d_k where e = Sum_{k>=1} d_k * k! is the factorial-base representation of e. So, the factors in this factorization are prime powers with exponents that are factorial numbers. Each factor in the factorization, p^(k!), can have a multiplicity d in the range [1, k], so this factorization of n is a product of numbers of the form (p^(k!))^d.
The number of factors counted with multiplicity (the sum of the multiplicities d in (p^(k!))^d) is given by this sequence (analogous to A001222 for the canonical prime factorization).
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) = 2.
MATHEMATICA
fdigsum[n_] := Module[{k = n, m = 2, r, s = 0}, While[{k, r} = QuotientRemainder[k, m]; k != 0 || r != 0, s += r; m++]; s]; f[p_, e_] := fdigsum[e]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) fdigsum(n) = {my(k = n, m = 2, r, s = 0); while([k, r] = divrem(k, m); k != 0 || r != 0, s += r; m++); s; }
a(n) = {my(e = factor(n)[, 2]); sum(i = 1, #e, fdigsum(e[i])); }
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Amiram Eldar, Oct 08 2024
STATUS
approved