# Computed 2018-09-01 in GP/PARI CALCULATOR Version 2.9.4 (released) amd64 running linux (x86-64/GMP-6.1.2 kernel) 64-bit version, compiled: Dec 19 2017, gcc version 7.3.0 (Ubuntu 7.3.0-1ubuntu1), threading engine: pthread # with the following program: # # v064413 = readvec("b064413_upto65539_terms_only.txt"); \\ From b-file of A064413 prepared previously. # 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