A342403 a(1) = 1; a(n) = -Sum_{d|n, d < n} d * a(d).

1, -1, -1, 1, -1, 4, -1, -3, 2, 6, -1, -24, -1, 8, 7, 21, -1, -38, -1, -58, 9, 12, -1, 288, 4, 14, -16, -108, -1, -180, -1, -315, 13, 18, 11, 930, -1, 20, 15, 1126, -1, -314, -1, -256, -116, 24, -1, -6960, 6, -154, 19, -354, -1, 1078, 15, 2940, 21, 30, -1, 6664, -1, 32, -198, 9765, 17

Offset

1,6

Links

Table of n, a(n) for n=1..65.

Mathematica

a[1] = 1; a[n_] := a[n] = -Sum[d a[d], {d, Most @ Divisors[n]}]; Table[a[n], {n, 1, 65}]

Prog

(PARI) a(n) = if (n==1, 1, -sumdiv(n, d, if (d<n, d*a(d)))); \\ Michel Marcus, Mar 11 2021
(Python)
from sympy import divisors
def A342403(n): return 1 if n == 1 else -sum(d*A342403(d) for d in divisors(n) if d < n) # Chai Wah Wu, Mar 11 2021

Crossrefs

Cf. A008683, A055615, A165552, A321384.

Sequence in context: A197840 A360995 A344439 * A019633 A067277 A177951

Adjacent sequences: A342400 A342401 A342402 * A342404 A342405 A342406

Keyword

sign

Author

Ilya Gutkovskiy, Mar 10 2021

Status

approved