OFFSET
0,3
COMMENTS
a(n) is the number of partitions of n*(9*n+1)/2 into n distinct parts <=9*n.
EXAMPLE
a(2) = 9 because there are 9 2-element subsets that can be chosen from {1,2,...,18} having element sum 19: {1,18}, {2,17}, {3,16}, {4,15}, {5,14}, {6,13}, {7,12}, {8,11}, {9,10}.
MAPLE
b:= proc(n, i, t) option remember;
`if`(i<t or n<t*(t+1)/2 or n>t*(2*i-t+1)/2, 0,
`if`(n=0, 1, b(n, i-1, t) +`if`(n<i, 0, b(n-i, i-1, t-1))))
end:
a:= n-> b(n*(9*n+1)/2, 9*n, n):
seq(a(n), n=0..20);
MATHEMATICA
b[n_, i_, t_] /; i<t || n<t(t+1)/2 || n>t(2i-t+1)/2 = 0; b[0, _, _] = 1;
b[n_, i_, t_] := b[n, i, t] = b[n, i-1, t] + If[n<i, 0, b[n-i, i-1, t-1]];
a[n_] := b[n(9n+1)/2, 9n, n];
a /@ Range[0, 10] (* Jean-François Alcover, Dec 07 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jan 18 2012
STATUS
approved