A332508 a(n) = Sum_{d|n} binomial(n+d-2, n-1).

1, 3, 7, 25, 71, 280, 925, 3561, 12916, 49346, 184757, 710255, 2704157, 10427747, 40119781, 155288897, 601080391, 2334714319, 9075135301, 35352181116, 137846759282, 538302226628, 2104098963721, 8233718962365, 32247603703576, 126412458920775, 495918551104687

Offset

1,2

Links

Seiichi Manyama, Table of n, a(n) for n = 1..1000

Formula

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

a(n) ~ 4^(n-1) / sqrt(Pi*n). - Vaclav Kotesovec, Aug 04 2022

Mathematica

Table[DivisorSum[n, Binomial[n + # - 2, n - 1] &], {n, 1, 27}]
Table[SeriesCoefficient[Sum[x^k/(1 - x^k)^n, {k, 1, n}], {x, 0, n}], {n, 1, 27}]

Prog

(PARI) a(n) = sumdiv(n, d, binomial(n+d-2, n-1)); \\ Michel Marcus, Feb 14 2020

Crossrefs

Cf. A000005, A000203, A007437, A059358, A073570, A101289, A308814, A332470.

Sequence in context: A148722 A148723 A148724 * A309911 A148725 A148726

Adjacent sequences: A332505 A332506 A332507 * A332509 A332510 A332511

Keyword

nonn

Author

Ilya Gutkovskiy, Feb 14 2020

Status

approved