OFFSET
1,6
COMMENTS
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
FORMULA
Mobius transform of A010051, the characteristic function of the primes.
Row sums of triangle A143518.
a(n) = Sum_{a*b*c=n} omega(a)*mu(b)*mu(c). - Benedict W. J. Irwin, Mar 02 2022
EXAMPLE
MATHEMATICA
Table[Sum[MoebiusMu[n/d] Boole[PrimeQ@ d], {d, Divisors@ n}], {n, 89}] (* Michael De Vlieger, Jul 19 2017 *)
PROG
(Sage)
def A143519(n) :
D = filter(is_prime, divisors(n))
return add(moebius(n/d) for d in D)
[A143519(n) for n in (1..89)] # Peter Luschny, Feb 01 2012
(PARI) A143519(n) = sumdiv(n, d, isprime(d)*moebius(n/d)); \\ (After Luschny's Sage-code) - Antti Karttunen, Jul 19 2017
CROSSREFS
KEYWORD
sign
AUTHOR
Gary W. Adamson, Aug 22 2008
EXTENSIONS
More terms from R. J. Mathar, Jan 19 2009
STATUS
approved