A293574 a(n) = Sum_{k=0..n} n^(n-k)*binomial(n+k-1,k).

1, 2, 11, 82, 787, 9476, 139134, 2422218, 48824675, 1118286172, 28679699578, 814027423892, 25330145185646, 857375286365768, 31360145331198428, 1232586016712594010, 51805909208539809315, 2318588202311267591852, 110085368092924083334626, 5526615354023679440754396, 292501304641192746350100410

Offset

0,2

Comments

a(n) is the n-th term of the main diagonal of iterated partial sums array of powers of n (see example).

Links

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

Formula

a(n) = [x^n] 1/((1 - x)^n*(1 - n*x)).

a(n) ~ exp(1) * n^n. - Vaclav Kotesovec, Oct 16 2017

Example

For n = 2 we have:
----------------------------
0   1   [2]   3    4     5
----------------------------
1,  2,   4,   8,  16,   32, ... A000079 (powers of 2)
1,  3,   7,  15,  31,   63, ... A126646 (partial sums of A000079)
1,  4, [11], 26,  57,  120, ... A000295 (partial sums of A126646)
----------------------------
therefore a(2) = 11.

Mathematica

Join[{1}, Table[Sum[n^(n - k) Binomial[n + k - 1, k], {k, 0, n}], {n, 1, 20}]]
Table[SeriesCoefficient[1/((1 - x)^n (1 - n x)), {x, 0, n}], {n, 0, 20}]
Join[{1, 2}, Table[n^(2 n)/(n - 1)^n - Binomial[2 n, n + 1] Hypergeometric2F1[1, 2 n + 1, n + 2, 1/n]/n, {n, 2, 20}]]

Prog

(PARI) a(n) = sum(k=0, n, n^(n-k)*binomial(n+k-1, k)); \\ Michel Marcus, Oct 12 2017

Crossrefs

Cf. A000312, A001700, A031973, A032443, A100192, A100193.

Sequence in context: A215654 A330677 A209094 * A322644 A243408 A352655

Adjacent sequences: A293571 A293572 A293573 * A293575 A293576 A293577

Keyword

nonn

Author

Ilya Gutkovskiy, Oct 12 2017

Status

approved