%2017年11月10日09:56:38
%S 1,1,1,1,1,1,1,9,1,9,1,9,1,9,1,1,9,1,1,9,
%T 1,17,1,9,9,9,9,1,17,1,9,9,9,1,17,1,9,
%U 17,1,36,9,9,1,44,1,9,9
%N对于小于或等于N的立方根的N的除数,除数函数的立方体之和。
%H Robert Israel,n的表,a(n)表示n=1..10000</a>
%H Sary Drappeau,<a href=“http://arxiv.org/abs/1302.4318“>Propriétés multiplicatories des entiers friables translateés</a>,arXiv:1307.4250[math.NT](见第9页)。
%H B.Landreau,<a href=“http://blms.oxfordjournals.org/content/21/4/366.摘录“>范德科尔普特定理的新证明,伦敦数学学会(1989)21(4):366-368。doi:10.1112/blms/21.4366,参见引理(3)第1页。
%F a(n)=(和{d|n}d<=n^(1/3))τ(d)^3。
%F如果p是素数,a(p^k)=A000537(1+楼层(k/3))_罗伯特·伊斯雷尔,2016年11月30日
%e a(7)=1,因为7的除数是1和7;只有1小于7的立方根,tau(1^3)=1,所以总和是1。
%e a(8)=9,因为8的除数是1,2,4,8;8的立方根是2,所以只有1和2是小于或等于立方根的除数,这些除数的立方是1和8,它们加起来是9。
%p f:=过程(n)加(数字理论:-tau(d)^3,d=选择(t->(t^3<=n),数字理论:-除数(n)))结束过程:
%p映射(f,[1..100]美元);#_罗伯特·伊斯雷尔,2016年11月30日
%t表格[selDivs=Select[Range[Floor[n^(1/3)]],IntegerQ[n/#]&];总和[DivisorSigma[0,selDivs[[m]]^3,{m,Length[selDivs]}],{n,100}](*_Alonso del Arte_,2013年7月21日*)
%o(PARI)a(n)=汇总(n,d,(d^3<=n)*numdiv(d)^3)\\马库斯,2013年7月21日
%Y参考A000537、A007425、A224834。
%K nonn公司
%O 1,8型
%马库斯,2013年7月21日
|