|
| |
| |
|
|
|
0, 2, 3, -2, 5, -5, 7, 0, -3, -7, 11, 2, 13, -9, -8, 0, 17, 3, 19, 2, -10, -13, 23, 0, -5, -15, 0, 2, 29, 10, 31, 0, -14, -19, -12, 0, 37, -21, -16, 0, 41, 12, 43, 2, 3, -25, 47, 0, -7, 5, -20, 2, 53, 0, -16, 0, -22, -31, 59, -2, 61, -33, 3, 0, -18, 16, 67, 2, -26, 14, 71, 0, 73, -39, 5, 2, -18, 18, 79, 0, 0, -43, 83, -2, -22, -45, -32, 0
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,2
|
|
|
COMMENTS
| Equals row sums of triangle A143517 [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Aug 22 2008]
|
|
|
FORMULA
| A054525 * A061397 = Mobius transform of [0, 2, 3, 0, 5, 0, 7, 0, 0, 0, 11, 0, 13,...].
|
|
|
EXAMPLE
| a(4) = -2 = (0, -1, 0, 1) dot (0, 2, 3, 0), where (0, -1, 0, 1) = row 4 of the Mobius triangle A054525 and (0, 2, 3, 0) = the first 4 terms of A061397.
|
|
|
MAPLE
| A061397 := proc(n) if isprime(n) then n; else 0 ; fi ; end: A054525 := proc(n, k) if n mod k = 0 then numtheory[mobius](n/k); else 0; fi ; end: A137851 := proc(n) local k ; add(A061397(k)* A054525(n, k), k=1..n) ; end: seq(A137851(n), n=1..120) ; - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), May 23 2008
|
|
|
PROG
| (Sage)
def A137851(n) :
D = filter(is_prime, divisors(n))
return add(d*moebius(n/d) for d in D)
[A137851(n) for n in (1..88)] # Peter Luschny, Feb 01 2012
|
|
|
CROSSREFS
| Cf. A061397, A054525, A143517, A143519.
Sequence in context: A144307 A144310 A035361 * A141346 A095402 A086294
Adjacent sequences: A137848 A137849 A137850 * A137852 A137853 A137854
|
|
|
KEYWORD
| easy,sign
|
|
|
AUTHOR
| Gary W. Adamson (qntmpkt(AT)yahoo.com), Feb 14 2008
|
|
|
EXTENSIONS
| More terms from R. J. Mathar (mathar(AT)strw.leidenuniv.nl), May 23 2008
|
| |
|
|
