Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A346255
Sum of A336849 and its Dirichlet inverse.
4
2, 0, 0, 9, 0, 30, 0, 27, 25, 42, 0, -15, 0, 66, 70, 81, 0, -25, 0, 63, 110, 78, 0, 45, 49, 102, 125, -33, 0, -140, 0, 243, 130, 114, 154, 625, 0, 138, 170, 189, 0, -440, 0, 117, 175, 174, 0, -81, 121, 147, 190, -51, 0, 625, 182, 99, 230, 186, 0, 315, 0, 222, 275, 729, 238, -260, 0, 171, 290, -308, 0, 15, 0, 246, 245, -69, 286
OFFSET
1,1
FORMULA
a(n) = A336849(n) + A346254(n).
PROG
(PARI)
up_to = 16384;
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.
A003961(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ From A003961
A336849(n) = { my(u=A003961(n)); (u/gcd(u, sigma(u))); };
v346254 = DirInverseCorrect(vector(up_to, n, A336849(n)));
A346254(n) = v346254[n];
A346255(n) = (A336849(n)+A346254(n));
CROSSREFS
KEYWORD
sign
AUTHOR
Antti Karttunen, Jul 19 2021
STATUS
approved