%I#24 2022年3月2日11:54:21

%S 1，-1，-3，-1,0,10,8,12,1，-28，-29，-67，-51，-28,79163256343273136，

%电话-351、-649、-1446、-1751、-1889、-1453、-1241924513876081063903、，

%U 100543143、-5799、-20521、-37217、-53057、-65661、-66086、-54430、-156483671719122732

%N乘积{k>=1}（1-x^k）^sigma（k）的展开。

%H Seiichi Manyama，n的表格，n=0..1000时的a（n）</a>

%F A061256的卷积逆。

%对于n>0，F a（0）=1，a（n）=-（1/n）*Sum_{k=1..n}A001001（k）*a（n-k）。

%F G.F.：exp（-总和{k>=1}σ_2（k）*x^k/（k*（1-x^k）））_伊利亚·古特科夫斯基，2018年10月29日

%p（数字理论）：

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

%p d*σ（d），d=除数（j）*b（n-j），j=1..n）/n）

%p端：

%p a:=proc（n）选项记忆`如果`（n=0，1，

%p-加（b（n-i）*a（i），i=0..n-1））

%p端：

%p序列（a（n），n=0..45）；#_Alois P.Heinz，2017年6月8日

%tb[n_]：=b[n]=如果[n==0，1，Sum[Sum[d*DivisorSigma[1，d]，{d，

%t除数[j]}]*b[n-j]，{j，1，n}]/n]；

%ta[n_]：=a[n]=如果[n==0，1，-和[b[n-i]*a[i]，{i，0，n-1}]]；

%t表[a[n]，{n，0，45}]（*_Jean-François Alcover_，2022年3月2日，在_Alois P.Heinz_*之后）

%Y参考A001001、A061256。

%Y积{k>=1}（1-x^k）^sigma_m（k）：A288098（m=0），此序列（m=1），A288389（m=2），A2 88392（m=3）。

%K符号

%0、3

%A _Seiichi Manyama，2017年6月8日