OFFSET
0,9
FORMULA
G.f.: Sum_{k>=0} k! * x^(k*(3*k - 2)) / Product_{j=1..k} (1 - x^(6*j)).
EXAMPLE
a(21) = 6 because we have [13, 7, 1], [13, 1, 7], [7, 13, 1], [7, 1, 13], [1, 13, 7] and [1, 7, 13].
MATHEMATICA
nmax = 83; CoefficientList[Series[Sum[k! x^(k (3 k - 2))/Product[1 - x^(6 j), {j, 1, k}], {k, 0, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 23 2020
STATUS
approved