# This b-file was computed by Antti Karttunen, 2021-11-26 # 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: # # up_to = 20000; # v064413 = readvec("b064413_upto65539_terms_only.txt"); \\ Data prepared with _Chai Wah Wu_'s Dec 08 2014 Python-program given in A064413. # A064413(n) = v064413[n]; # \\ Then its inverse A064664 was prepared: # m064664 = Map(); # for(n=1,65539,mapput(m064664,A064413(n),n)); # A064664(n) = mapget(m064664,n); # DirInverseCorrect(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = (-u[1]*sumdiv(n, d, if(d