|
nmax=32；系数列表[Series[Product[Product[Product[1/（1-x^（ijk）））^（i jk），{i，1，nmax}]，{j，1，nmax}]
nmax=32；系数列表[级数[乘积[1/（1-x^k）^（k和[DivisorSigma[0，d]，{d，Divisors[k]}]），{k，1，nmax}]，{x，0，nmax{]，x]
nmax=32；系数列表[Series[Exp[Sum[d^2 Sum[DivisorSigma[0，j]，{j，Divisors[d]}]，{d，Divisor[k]}]x^k/k，{k，1，nmax}]]，{x，0，nmax}]，x]
a[n]：=a[n]=如果[n==0，1，Sum[d^2 Sum[DivisorSigma[0，j]，{j，Divisors[d]}]，{d，Divisor[k]}]a[n-k]，{k，1，n}]/n]；表[a[n]，{n，0，32}]
nmax=50；A034718号=表[n*和[DivisorSigma[0，d]，{d，Divisors[n]}]，{n，1，nmax}]；s=1-x；Do[s*=和[二项式[A034718号[[k]]，j]*（-1）^j*x^（j*k），{j，0，nmax/k}]；s=展开[s]；s=取[s，Min[nmax+1，指数[s，x]+1，长度[s]]，{k，2，nmax}]；系数列表[系列[1/s，{x，0，nmax}]，x](*瓦茨拉夫·科泰索维奇2018年8月31日*）
| 