A359355 a(n) = A359107(2*n, n) = Sum_{j=0..n} Stirling2(2*n, j) = Sum_{j=0..n} A048993(2*n, j).
1, 1, 8, 122, 2795, 86472, 3403127, 164029595, 9433737120, 635182667816, 49344452550230, 4371727233798927, 437489737355466560, 49048715505983309703, 6116937802946210183545, 843220239072837883168510, 127757559136845878072576947, 21166434937698025552654090472
Offset: 0
Keywords
Links
- Alois P. Heinz, Table of n, a(n) for n = 0..288
Programs
-
Maple
b:= proc(n) option remember; `if`(n=0, 1, add(expand(b(n-j)*binomial(n-1, j-1)*x), j=1..n)) end: a:= n-> (p-> add(coeff(p, x, i), i=0..n))(b(2*n, 0)): seq(a(n), n=0..17); # Alois P. Heinz, Jun 13 2023 -
Mathematica
b[n_] := b[n] = If[n == 0, 1, Sum[Expand[b[n-j]*Binomial[n-1, j-1]*x], {j, 1, n}]]; a[n_] := With[{p = b[2*n]}, Sum[Coefficient[p, x, i], {i, 0, n}]]; Table[a[n], {n, 0, 17}] (* Jean-François Alcover, Jul 01 2025, after Alois P. Heinz *) -
PARI
a(n) = sum(j=0, n, stirling(2*n, j, 2)); \\ Michel Marcus, Dec 27 2022
Formula
a(n) = A102661(2n,n) for n >= 1. - Alois P. Heinz, Jun 13 2023
Comments