%I#9 2019年4月3日03:00:15

%S 1,1,3,7,11,25,51,87129286462760131221023470598884013884，

%电话：22577335455596185341126705194317293621435040641472971503，

%电话：1462483210816131244894474579654580995619231351867819809034283876254028663157039233

%N乘积{k>=1}（1+sigma（k）*x^k）的展开式，其中sigma是k（A000203）的除数之和。

%H Vaclav Kotesovec，<a href=“/A316962/b316962.txt”>n表，n=0..10000时的a（n）</a>

%F G.F.：exp（Sum_｛k>=1｝Sum_｛j>=1｝（-1）^（k+1）*西格玛（j）^k*x^（j*k）/k）。

%p与（数字理论）：a：=系列（mul（1+sigma（k）*x^k，k=1..100），x=0.39）：seq（系数（a，x，n），n=0..38）；#_保罗·拉瓦（Paolo P.Lava），2019年4月2日

%t nmax=38；系数列表[系列[积[（1+DivisorSigma[1，k]x^k），{k，1，nmax}]，{x，0，nmax{]，x]

%t nmax=38；系数列表[Series[Exp[Sum[Sum[（-1）^（k+1）Divisor Sigma[1，j]^k x ^（j k）/k，{j，1，nmax}]，{k，1，nmax}]]，{x，0，nmax}]，x]

%t a[n_]：=a[n]=如果[n==0，1，Sum[Sum[（-1）^（k/d+1）d除数Sigma[1，d]^（k/d），{d，除数[k]}]a[n-k]，{k，1，n}]/n]；表[a[n]，{n，0，38}]

%Y参见A000203、A180305、A192065、A279786、A316961。

%K nonn公司

%0、3

%A _Ilya Gutkovskiy_，2018年7月17日