OFFSET
1,2
LINKS
PROG
(PARI)
up_to = 11550;
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.
A032742(n) = if(1==n, n, n/vecmin(factor(n)[, 1]));
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
v359427 = DirInverseCorrect(vector(up_to, n, A358764(n)));
A359427(n) = v359427[n];
CROSSREFS
KEYWORD
sign,base
AUTHOR
Antti Karttunen, Jan 02 2023
STATUS
approved