A133723 Column 4 of triangle in A133721.

0, 0, 0, 1, 1, 1, 1, 1, 13, 7, 3, 1, 81, 25, 6, 1, 325, 65, 10, 1, 995, 140, 15, 1, 2541, 266, 21, 1, 5698, 462, 28, 1, 11586, 750, 36, 1, 21825, 1155, 45, 1, 38665, 1705, 55, 1, 65131, 2431, 66, 1, 105183, 3367, 78, 1, 163891, 4550, 91, 1

Offset: 1

Author

N. J. A. Sloane, Dec 30 2007

Keywords

• nonn

Programs

Maple

A133723 := proc(n)
        A133721(n,4) ;
end proc:
seq(A133723(n),n=1..60) ; # R. J. Mathar, Nov 23 2011

Mathematica

A133713[l_, cl_] := Module[{g, k, s}, g = 1; For[k = 1, k <= cl + 1, k++, s = Sum[Binomial[Binomial[l, k + 1] + i - 1, i]*t^(i*k), {i, 0, Ceiling[ cl/k]}]; g = g*s]; SeriesCoefficient[g, {t, 0, cl}]];
a[m_, n_] := A133713[Ceiling[m/n], n*Ceiling[m/n] - m];
Table[a[m, 4], {m, 1, 60}] (* Jean-François Alcover, Apr 03 2020, after R. J. Mathar *)