# Computed 2018-08-11 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: # # A002487(n) = { my(a=1, b=0); while(n>0, if(bitand(n, 1), b+=a, a+=b); n>>=1); (b); }; \\ From Charles R Greathouse IV, Oct 05 2008 # memo = Map(); # A317931perA317932(n) = if(1==n,n,if(mapisdefined(memo,n),mapget(memo,n),my(v = (A002487(n)-sumdiv(n,d,if((d>1)&&(d