%I#33 2022年6月19日15:23:21
%S 1,5,282613126466823544167774773874205171000003130,
%电话285311706128916100495200302875106592254111112006826381564,
%电话:43789389038086252818446744073726329093827240261866367641783964080752969250426011978419655660313589123980
%f(N)=N^N(A000312)的N逆Moebius变换。
%H Seiichi Manyama,n表,n=1..386的a(n)
%F a(n)=总和。
%F G.F.:和{n>=1}n^n*x^n/(1-x^n)_Paul D.Hanna,2009年10月27日
%F A023879的对数导数。-_Paul D.Hanna,2012年9月5日
%e n=6:除数={1,2,3,6};1^1+2^2+3^3+6^6=1+4+27+46656=46688=a(6)。
%t a[n_]:=除数和[n,#^#&];阵列[a,19](*_Jean-François Alcover_,2015年12月23日*)
%o(PARI)向量(17,n,sumdiv(n,d,d^d))
%o(PARI){a(n)=polcoeff(总和(m=1,n,m^m*x^m/(1-x^m+x*o(x^n))),n)}\\_Paul D.Hanna,2009年10月27日
%o(PARI)a(n)=汇总(n,d,d^d);\\_Joerg Arndt_,2013年4月14日
%o(Python)
%o来自sympy导入除数
%o def A062796(n):返回和(d**d表示d的除数(n,生成器=True))#_Chai Wah Wu_,2022年6月19日
%Y参考A000312、A023879。
%K nonn公司
%O 1,2号机组
%A _Labos Elemer,2001年7月19日
|