%I#24 2023年3月24日19:07:51

%S 1,1,1,9,1,1,1,9,28,1,1,9/1,1,1,1,73,1,28,1,9,1,1,1,1,1,9126,1,28,9,11,1,1，

%电话：73,1,1,1252,1,1,1,9,1,1,9,1,1,9,28,1,1,73344126,1,9,1,28,1,91,1,1,9，

%U 1,1,28585,1,1,1,9,1,1252,1,1126,9,11,1,1,73

%N满足d^2|N的除数d的立方体d^3的和。

%C莫比乌斯变换是1，0，0，8，0，0,0，27,0，0…=n ^（3/2）*A010052（n）。

%H Winston de Greef，n表，n=1..10000的a（n）</a>

%H A.Dixit、B.Maji和A.Vatwani，<A href=“https://arxiv.org/abs/2303.09937“>广义除数函数sigma_z^k（n）的Voronoi求和公式，arXiv:2303.09937[math.NT]，2023，sigma（z=3，k=2，n）。

%F a（n）=和{d^2|n}d^3。

%F与a（p^e）的乘积=（p^（3*（floor（e/2）+1））-1）/（p^3-1）_Amiram Eldar，2023年3月24日

%p gsigma：=进程（n，z，k）

%p局部a、d；

%p a：=0；

%数理论中d的p（除数）（n）do

%p如果modp（n，d^k）=0，则

%p a：=a+d ^z；

%p end if；

%p端do：

%p a；

%p端程序：

%p序列（gsigma（n，3,2），n=1..80）；

%t f[p_，e_]：=（p^（3*（楼层[e/2]+1））-1）/（p^3-1）；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2023年3月24日*）

%o（PARI）a（n）=sumdiv（n，d，if（issquare（d），sqrtint（d）^3））；\\_米歇尔·马库斯，2023年3月24日

%o（Python）

%o来自math导入prod

%o来自症状输入因子int

%o定义A361794（n）：返回prod（（p**（3*（e>>1）+3）-1）//（p**3-1）for p，e in factorint（n）.items（））#_Chai Wah Wu_，2023年3月24日

%Y参考A010052、A035316、A069290、A351308、A333843和A333844。

%K nonn，mult，容易

%O 1,4型

%A.R.J.Mathar_，2023年3月24日