A321260 a(n) = [x^n] Product_{k>=1} 1/(1 - x^k)^(sigma_n(k)-k^n).

1, 0, 1, 1, 18, 2, 861, 132, 106024, 40910, 72980055, 6838271, 228282942581, 27620223647, 2050169324675668, 352809815149813, 87174966874755673105, 6798293425492905407, 18318448554980083512011863, 1187839217207171380193247, 11258918803635775614062752424535

Offset

0,5

Links

Table of n, a(n) for n=0..20.

Formula

a(n) = [x^n] exp(Sum_{k>=1} sigma_(n+1)(k)*x^(2*k)/(k*(1 - x^k))).

Mathematica

Table[SeriesCoefficient[Product[1/(1 - x^k)^(DivisorSigma[n, k] - k^n), {k, 1, n}], {x, 0, n}], {n, 0, 20}]
Table[SeriesCoefficient[Exp[Sum[DivisorSigma[n + 1, k] x^(2 k)/(k (1 - x^k)), {k, 1, n}]], {x, 0, n}], {n, 0, 20}]

Crossrefs

Cf. A283333, A318783, A318784, A319647, A321258, A321261.

Sequence in context: A040321 A040322 A040323 * A052453 A040315 A040316

Adjacent sequences: A321257 A321258 A321259 * A321261 A321262 A321263

Keyword

nonn

Author

Ilya Gutkovskiy, Nov 01 2018

Status

approved