OFFSET
1,2
LINKS
PROG
(PARI)
up_to = 32768;
v064413 = readvec("b064413_upto65539_terms_only.txt"); \\ From precomputed file.
A064413(n) = v064413[n];
\\ Then its inverse A064664 is prepared:
m064664 = Map();
for(n=1, 65539, mapput(m064664, A064413(n), n));
A064664(n) = mapget(m064664, n);
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.
v323411 = DirInverse(vector(up_to, n, A064664(n)));
A323411(n) = v323411[n];
CROSSREFS
KEYWORD
sign
AUTHOR
Antti Karttunen, Jan 13 2019
STATUS
approved