%I#11 2018年6月29日22:16:12

%S 1,4,7,17,25,58,8717728952886015502486425769101147418335，

%电话：29941473317581911888718733829078445290469605810712341632947，

%电话：2487504375961356764248512310127449031897583928194293416911576151639490379785

%N和{k>=1}（-1+乘积{j>=1}1/（1-x^（k*j））^j）的展开。

%A000219的C Moebius逆变换。

%H Alois P.Heinz，n的表格，n=1..10000的a（n）</a>

%H N.J.A.Sloane，转换</a>

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/PlanePartition.html“>平面分区</a>

%F G.F.：总和{k>=1}A000219（k）*x^k/（1-x^k）。

%F a（n）=和{d|n}A000219（d）。

%p b:=proc（n）选项记忆`如果`（n=0，1，添加(

%p b（n-j）*numtheory[σ][2]（j），j=1..n）/n）

%p端：

%p a:=n->加（b（d），d=numtheory[除数]（n））：

%p序列（a（n），n=1..40）；#_Alois P.Heinz，2018年6月21日

%t最大值=37；Rest[系数列表[级数[和[-1+积[1/（1-x^（kj））^j，{j，1，nmax}]，{k，1，nmax}]

%tb[n_]：=b[n]=系列系数[乘积[1/（1-x^k）^k，{k，1，n}]，{x，0，n}]；a[n]：=a[n]=级数系数[和[b[k]x^k/（1-x^k），{k，1，n}]，{x，0，n}]；表[a[n]，{n，37}]

%tb[0]=1；b[n]：=b[n]=和[b[n-j]除数Sigma[2，j]，{j，n}]/n；a[n]：=a[n]=和[b[d]，{d，除数[n]}]；表[a[n]，{n，37}]

%Y参见A000219、A047966、A0479968、A300275、A302550。

%K nonn公司

%O 1,2号机组

%A _Ilya Gutkovskiy_，2018年6月20日