%I#14 2022年2月12日09:26:45

%S 1,6573041601562647450117650266240532170101156901771562，

%电话：3036800482681076472501140698017039360241375703459105047045882，

%电话：65004160858845001151530148035890194355200441562503137426503879519304894240005948232741453700

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

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

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

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

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

%F与a（p^e）的乘积=p^（6*e）+p^_塞巴斯蒂安·卡尔松，2022年2月8日

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

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

%F求和{k=1..n}a（k）~n^7*zeta（7）/（7*zeta（14））=2606175*n^7*zeta（七）/（2*Pi^14）。

%F和{k>=1}1/a（k）=积{素数p}（1+p^6/（p^12-1））=1.017096592895596077024247497998914920118274875188346774241790304…（结束）

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

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

%o（PARI）用于（n=1100，打印1（目录（p=2，n，（1+X）/（1-p^6*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日