%I#39 2022年9月8日08:45:03

%S 1,6,7,21,9,42,11,58,30,54,15147,17,66,63141,21180,23189,77,90,27，

%电话：406,54102106231,33378,35318105126,99630,41138119522,45，

%电话：462,47315270162,51987,86324147357,5763135638161198,631323

%N和tau^2（N）的N Dirichlet卷积。

%A000027和A035116的C狄利克雷卷积。

%A060724.-的C逆Mobius变换_R.J.Mathar，2011年10月15日

%H Vincenzo Librandi，n的表，n=1..5000的a（n）</a>

%F a（n）=总和{i|n，j|n}σ（i）*σ（j）/σ（lcm（i，j）），其中σ（n）=n的除数之和。

%F a（n）=和{i|d，j|d}σ（gcd（i，j））；

%F a（n）=和{d|n}d*tau（n/d）^2，其中tau（n）=n的除数。

%F与a（p^e）的乘法=（1-p^（3+e）-p^（2+e）+e^2+4*p^2+p^2*e^2+2*e-3*p+4*p^2*e-6*e*p-2*e*e^2*p）/（1-p）^3。

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

%F G.F.：Sum_｛k>=1｝tau（k）^2*x^k/（1-x^k）^2。-_伊利亚·古特科夫斯基，2018年11月2日

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

%F a（n）=和{d|n}τ（d^2）*σ（n/d），其中τ（n）=n的除数，σ（n

%t a[n_]：=和[DivisorSigma[1，i]*除数Sigma[1]，j]/DivisorSigma[1，LCM[i，j]]，{i，除数[n]}，{j，除数[n]}]；表[a[n]，{n，1，60}]（*Jean-François Alcover_，2013年3月26日*）

%o（PARI）a（n）=sumdiv（n，d，d*numdiv（n/d）^2）；\\_米歇尔·马库斯，2018年11月3日

%o（岩浆）[&+[d*#除数（底面（n/d））^2:d，单位为除数（n）]：n，单位为[1..60]]；//_Marius A.Burtea，2019年8月25日

%Y参考A000203、A060648。

%Y参考A000005、A060724、A064950、A062369、A0623658、A062380。

%K nonn，多个

%O 1,2号机组

%A_Vladeta Jovovic_，2001年7月7日