# Computed 2018-08-12 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: # # A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415 # memo317835 = Map(); # A317835aux(n) = if(1==n,n,if(mapisdefined(memo317835,n),mapget(memo317835,n),my(v = (A003415(n)-sumdiv(n,d,if((d>1)&&(d