OFFSET
1,3
LINKS
FORMULA
PROG
(PARI)
v064413 = readvec("b064413_upto65539_terms_only.txt"); \\ From b-file of A064413 prepared beforehand.
A064413(n) = v064413[n];
m064664 = Map();
for(n=1, 65539, mapput(m064664, A064413(n), n));
A064664(n) = mapget(m064664, n);
up_to = (2^14);
DirSqrt(v) = {my(n=#v, u=vector(n)); u[1]=1; for(n=2, n, u[n]=(v[n]/v[1] - sumdiv(n, d, if(d>1&&d<n, u[d]*u[n/d], 0)))/2); u}; \\ From A317937.
v318664_65 = DirSqrt(vector(up_to, n, A064664(n)));
A318664(n) = numerator(v318664_65[n]);
A318665(n) = denominator(v318664_65[n]);
CROSSREFS
KEYWORD
sign,frac
AUTHOR
Antti Karttunen, Sep 01 2018
STATUS
approved