OFFSET
1,4
COMMENTS
The maximum exponent in the prime factorization of 1 is considered to be A051903(1) = 0.
LINKS
FORMULA
a(n) = Sum_{d|n} [m(d) = m(n)], where m(n) = A051903(n) and [] is the Iverson bracket.
If n = Product_{i} p_i^e_i (where p_i are distinct primes), then a(n) = Product_{i} (e_i + [e_i < Max_{i}(e_i)]).
a(n) <= 1 if and only if n is squarefree (A005117), and a(n) = 0 only for n = 1.
a(n) < A000005(n).
MATHEMATICA
a[n_] := Module[{e = FactorInteger[n][[;; , 2]], m}, m = Max[e]; Times@@ ((Min[#, m-1] & /@ e) + 1)]; a[1] = 0; Array[a, 100]
PROG
(PARI) a(n) = if(n == 1, 0, my(e = factor(n)[, 2], m = vecmax(e)); vecprod(apply(x -> 1 + min(x, m-1), e)));
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Sep 11 2024
STATUS
approved