OFFSET
1,3
MAPLE
b:= proc(n, i) option remember; `if`(n=0, {1}, `if`(i<1, {},
{seq(map(x-> x*i!^j*j!, b(n-i*j, i-1))[], j=0..n/i)})) end:
a:= n-> nops(b(n, n)):
seq(a(n), n=1..40); # Alois P. Heinz, Apr 13 2012
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0, {1}, If[i < 1, {}, Table[# i!^j j!& /@ b[n - i j, i - 1], {j, 0, n/i}] // Flatten // Union]];
a[n_] := Length[b[n, n]];
Array[a, 40] (* Jean-François Alcover, Nov 09 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladeta Jovovic, Feb 23 2005
EXTENSIONS
More terms from David Wasserman, Apr 11 2008
STATUS
approved