%I#20 2022年11月15日17:27:11

%S 1257656265792390626168643457648021684275243053282100390882，

%电话：214358882431727104815730722148155411425632878124311744512，

%电话：697575744211064693474169835630422570006579237828630724550902326747831098528211052213862415258281250

%N a（N）=N^8*乘积{p|N，p素数}（1+1/p^8）。

%C n的无平方除数的除数补的8次幂之和。

%H Sebastian Karlsson，n表，n（n）表示n=1..10000</a>

%F a（n）=和{d|n}d^8*mu（n/d）^2。

%Fa（n）=n^8*Sum_{d|n}mu（d）^2/d^8。

%F与a（p^e）的乘积=p^（8*e）+p^_Sebastian Karlsson_，2022年2月8日

%F From _Vaclav Kotesovec_，2022年2月12日：（开始）

%F Dirichlet g.F.：zeta（s）*zeta（s-8）/zeta（2*s）。

%F和{k=1..n}a（k）~n^9*zeta（9）/（9*zeta（18））=4331032831125*n^9*zeta（九）/（43867*Pi^18）。

%F和{k>=1}1/a（k）=积{素数p}（1+p^8/（p^16-1））=1.004062071480173688638170669970682370243496458304179434830922739661777…（结束）

%F a（n）=J_16（n）/J_8（n）=J_16_恩里克·佩雷斯·埃雷罗，2022年11月14日

%tf[p_，e_]：=p^（8*e）+p^；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；数组[a，25]（*_Amiram Eldar_，2022年2月8日*）

%o（PARI）a（n）=sumdiv（n，d，moebius（n/d）^2*d^8）；

%o（PARI）用于（n=1100，打印1（目录（p=2，n，（1+X）/（1-p^8*X））[n]，“，”））\\_Vaclav Kotesovec_，2022年2月12日

%Y参考A008683（亩）。

%形式为n^k*Product_{p|n，pprime}（1+1/p^k）的Y序列，k=0..10:A034444（k=0），A001615。

%K nonn，多个

%O 1,2号机组

%韦斯利·伊万·赫特，2022年2月6日