#2018-08-11在GP/PARI CALCULATOR版本2.9.4(发布)和运行linux(x86-64/GMP-6.1.2内核)64位版本的d64中计算,编译时间:2017年12月19日,gcc版本7.3.0(Ubuntu 7.3.0-1ubuntu1),线程引擎:pthread#使用以下程序:##A005187(n)={my(s=n);而(n>>=1,s+=n),s;};#memo=地图();#A317927perA317928(n)=如果(1==n,n,if(mapisdefined(memo,n),mapget(memo(n)),my(v=(A005187(n)-sumdiv(n,d,if)((d>1)&&(d)-->