%我#2020年9月14日02:55:58

%S 1,9,28,651262523445137301134133218202198309635284097，

%电话：4914657068608190963211988121681436415626197821968422360，

%电话：2439031752297923276937296442264334447450506546174061544638689228668879508

%N的酉因子的立方体和。

%n的幺正除数是一个除数d，使得gcd（d，n/d）=1。

%H Amiram Eldar，n表，n的a（n）表示n=1..10000（Harvey P.Dale的术语1..1000）

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

%F如果n=乘积（p_j^k_j），则a（n）=乘积（1+p_j^（3*k_j））。-_伊利亚·古特科夫斯基，2018年11月4日

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

%e 6的幺正因子是1、2、3和6，因此a（6）=252。

%t scud[n_]：=总数[Select[Divisors[n]，CoprimQ[#，n/#]&]^3]；阵列[scud，40]（*哈维·P·戴尔，2016年10月16日*）

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

%o（PARI）A034677_vec（长度）={

%o a000012=方向（p=2，len，1/（1-X））；

%o a000578=方向（p=2，长度，1/（1-p^3*X））；

%o a000578x=方向（p=2，长度，1-p^3*X^2）；

%o dirmul（目录（a000012、a000578），a000578x）

%o}（o）

%o A034677_vec（70）/*通过D.g.f.，R.J.Mathar_，2011年3月5日*/

%Y参考A034444、A034448。

%A286880的Y行n=3。

%K非n，复数

%O 1,2号机组

%A _弗里德曼_