|
|
|
|
0, 1, 1, -1, 1, -2, 1, 0, -1, -2, 1, 1, 1, -2, -2, 0, 1, 1, 1, 1, -2, -2, 1, 0, -1, -2, 0, 1, 1, 3, 1, 0, -2, -2, -2, 0, 1, -2, -2, 0, 1, 3, 1, 1, 1, -2, 1, 0, -1, 1, -2, 1, 1, 0, -2, 0, -2, -2, 1, -1, 1, -2, 1, 0, -2, 3, 1, 1, -2, 3, 1, 0, 1, -2, 1, 1, -2, 3, 1, 0, 0, -2, 1, -1, -2, -2, -2, 0, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,6
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Mobius transform of A010051, the characteristic function of the primes.
|
|
EXAMPLE
|
a(4) = -1 since row 4 of triangle A043518 = (0, -1, 0, 0).
a(4) = -1 = (0, -1, 0, 1) dot (0, 1, 1, 0), where (0, -1, 0, 1) = row 4 of A054525 and A010051 = (0, 1, 1, 0, 1, 0, 1, 0, ...).
|
|
MATHEMATICA
|
Table[Sum[MoebiusMu[n/d] Boole[PrimeQ@ d], {d, Divisors@ n}], {n, 89}] (* Michael De Vlieger, Jul 19 2017 *)
|
|
PROG
|
(Sage)
D = filter(is_prime, divisors(n))
return add(moebius(n/d) for d in D)
(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
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|