%I#46 2020年9月14日02:58:31

%S 1,3,4,7,6,12,8,7,13,18,12,28,14,24,7,18,39,20,42,32,36,24,28,31，

%第42,13,56,30,7,7,48,54,48,91,38,60,56,42,42,96,44,84,78,72,48,28页，

%U 57,93,72,98,54,39,72,56,80,90,60168,62,96104,7,84144,68126,96144,72

%N的立方因子之和。

%C n的立方核的除数之和（见第一个公式）。

%H Reinhard Zumkeller，n的表，n=1..10000的a（n）</a>

%F a（n）=A000203（A007948（n））。

%F a（n）<=A073183（n）。

%F与a（p）=1+p相乘，a（p^e）=1+p+p^2，e>1_Christian G.Bower_，2005年5月18日

%F a（n）=总和（A212793（A027750（n，k））*A027750_Reinhard Zumkeller，2012年5月27日

%F Dirichlet g.F.：zeta（s）*zeta（s-1）/zeta（3s-3）_R.J.Mathar，2011年4月12日

%F和{k=1..n}a（k）~Pi^2*n^2/（12*Zeta（3））_Vaclav Kotesovec_，2019年2月1日

%e 56的除数是{1，2，4，7，8，14，28，56}，8=2^3和56=7*2^3，因此a（56）=1+2+4+7+14+28=56。

%p charFfree：=进程（n，t）局部f；对于ifactors（n）[2]中的f，如果op（2，f）>=t，则返回0；结束条件：；end-do：返回1；结束过程：

%p A073185:=过程（n）加（d*charFfree（d，3），d=数值[除数]（n））；结束程序：#R.J.Mathar，2011年4月12日

%t nn=71；f[list_，i_]：=列表[[i]]；a=表格[If[Max[FactorInteger[n][[All，2]]<=2，n，0]，{n，1，nn}]；b=表[1，{nn}]；选择[Table[DirichletConvolve[f[a，n]，f[b，n]、n，m]，{m，1，nn}]，#>0&]（*_Geoffrey Critzer_，2015年3月22日*）

%t f[p_，e_]：=1+p+如果[e>1，p^2，0]；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2020年9月14日*）

%o（哈斯克尔）

%o a073185=总和。过滤器（（==1）。a212793）。a027750_低

%o——Reinhard Zumkeller，2012年5月27日

%o（PARI）a（n）={my（f=factor

%Y参见A000203、A073184、A004709、A073182、A073181、A048250。

%K nonn，多个

%O 1,2号机组

%A _Reinhard Zumkeller，2002年7月19日

%E 2015年2月6日，阿尔瓦尔·伊比亚斯删除了错误评论