login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

A346238
Sum of A005187 and its Dirichlet inverse.
2
2, 0, 0, 9, 0, 24, 0, 15, 16, 48, 0, 8, 0, 66, 64, 31, 0, 32, 0, 4, 88, 114, 0, 50, 64, 138, 64, 7, 0, -116, 0, 63, 152, 192, 176, 80, 0, 210, 184, 90, 0, -138, 0, -1, 80, 252, 0, 94, 121, -6, 256, -5, 0, 28, 304, 125, 280, 324, 0, 390, 0, 342, 136, 127, 368, -276, 0, -20, 336, -386, 0, 126, 0, 426, 24, -17, 418, -360, 0
OFFSET
1,1
LINKS
FORMULA
a(n) = A005187(n) + A346237(n).
PROG
(PARI)
up_to = 65537;
DirInverseCorrect(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = (-u[1]*sumdiv(n, d, if(d<n, v[n/d]*u[d], 0)))); (u) }; \\ Compute the Dirichlet inverse of the sequence given in input vector v.
A005187(n) = { my(s=n); while(n>>=1, s+=n); s; };
v346237 = DirInverseCorrect(vector(up_to, n, A005187(n)));
A346237(n) = v346237[n];
A346238(n) = (A005187(n)+A346237(n));
CROSSREFS
Cf. also A300244.
Sequence in context: A345065 A323364 A323403 * A347085 A347091 A347229
KEYWORD
sign
AUTHOR
Antti Karttunen, Jul 13 2021
STATUS
approved