OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..200
FORMULA
a(n) = A292746(2n,n).
a(n) ~ 2^(2*n) * n^(n-1/2) / (sqrt(2*Pi*(1-c)) * exp(n) * c^n * (2-c)^n), where c = -LambertW(-2*exp(-2)) = -A226775 = 0.40637573995995990767695812412483975821... - Vaclav Kotesovec, Sep 28 2017
EXAMPLE
a(2) = 8: 21a1b, 1a1a1a1b, 1a1a1b1a, 1a1a1b1b, 1a1b1a1a, 1a1b1a1b, 1a1b1b1a, 1a1b1b1b (the two kinds of 1's are denoted by 1a and 1b).
MAPLE
f:= (n, k)-> add(Stirling2(n, j), j=0..k):
b:= proc(n, i, k) option remember; `if`(n=0 or i<2,
f(n, k), add(b(n-i*j, i-1, k), j=0..n/i))
end:
a:= n-> b(2*n$2, n)-b(2*n$2, n-1):
seq(a(n), n=0..20);
MATHEMATICA
f[n_, k_] := Sum[StirlingS2[n, j], {j, 0, k}];
b[n_, i_, k_] := b[n, i, k] = If[n==0 || i<2, f[n, k], Sum[b[n - i*j, i-1, k], {j, 0, n/i}]];
a[n_] := b[2n, 2n, n] - b[2n, 2n, n-1];
a /@ Range[0, 20] (* Jean-François Alcover, Dec 12 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Sep 22 2017
STATUS
approved