login
A323893
Dirichlet inverse of A048673.
10
1, -2, -3, -1, -4, 4, -6, -2, -4, 5, -7, 3, -9, 7, 6, -4, -10, 8, -12, 4, 8, 8, -15, 8, -9, 10, -12, 6, -16, 5, -19, -8, 9, 11, 9, 8, -21, 13, 11, 11, -22, 11, -24, 7, 16, 16, -27, 20, -25, 18, 12, 9, -30, 32, 10, 17, 14, 17, -31, 6, -34, 20, 24, -16, 12, 14, -36, 10, 17, 20, -37, 16, -40, 22, 27, 12, 12, 20, -42, 28, -36, 23, -45, 12, 13
OFFSET
1,2
LINKS
PROG
(PARI)
up_to = 20000;
DirInverse(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = -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.
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
A048673(n) = (A003961(n)+1)/2;
v323893 = DirInverse(vector(up_to, n, A048673(n)));
A323893(n) = v323893[n];
CROSSREFS
Sequence in context: A324752 A298637 A034867 * A329721 A193790 A055446
KEYWORD
sign
AUTHOR
Antti Karttunen, Feb 08 2019
STATUS
approved