%I#14 2022年2月12日09:29:35

%S 151319684262619531261009789240353608134479872387440172，

%电话：100195363823579476925170120704106044993742070140090438445332184，

%电话：6885369446411858787649819875680823632268769778051300026265679432041987212096271659961801526614642647101800448

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

%C n的平方除数的除数补码的9次方的和。

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

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

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

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

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

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

%F和{k=1..n}a（k）~n^10*zeta（10）/（10*zeta（20））=3273645375*n^10/（349222*Pi^10）。

%F和{k>=1}1/a（k）=Product_{primes p}（1+p^9/（p^18-1））=1.00204575331916689985388864168116922608947780516939765639888137700557…（结束）

%tf[p_，e_]：=p^（9*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^9）；

%o（PARI）用于（n=1100，打印1（目录（p=2，n，（1+X）/（1-p^9*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=1），A065958（k=2），AO65959（k=3），C065960（k=4），A351300。

%K nonn，多个

%O 1,2号机组

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