OFFSET
1,6
LINKS
FORMULA
Additive with a(p^e) = A010054(e). - Antti Karttunen, Jul 08 2017
Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B + C), where B is Mertens's constant (A077761) and C = -P(2) + Sum_{k>=2} (P(k*(k+1)/2) - P(k*(k+1)/2 + 1)) = -0.34517646457715166126..., where P(s) is the prime zeta function. - Amiram Eldar, Sep 28 2023
EXAMPLE
The prime-factorization of 360 is 2^3 *3^2 *5^1. There are two exponents in this factorization which are triangular numbers, 1 and 3. So a(360) = 2.
MATHEMATICA
f[n_] := Length @ Select[Last /@ FactorInteger[n], IntegerQ[Sqrt[8# + 1]] &]; Table[f[n], {n, 110}] (* Ray Chandler, Nov 19 2006 *)
PROG
(PARI)
A010054(n) = issquare(8*n + 1); \\ This function from Michael Somos, Apr 27 2000.
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Leroy Quet, Nov 18 2006
EXTENSIONS
Extended by Ray Chandler, Nov 19 2006
STATUS
approved