A093786 Hook products of all partitions of 8.

448, 576, 576, 630, 630, 720, 720, 960, 1152, 1152, 1440, 1440, 1920, 1920, 2016, 2016, 2880, 2880, 5760, 5760, 40320, 40320

Offset

1,1

Links

Table of n, a(n) for n=1..22.

Formula

a(n) = 8!/A003872(23-n).

Maple

H:=proc(pa) local F, j, p, Q, i, col, a, A: F:=proc(x) local i, ct: ct:=0: for i from 1 to nops(x) do if x[i]>1 then ct:=ct+1 else fi od: ct; end: for j from 1 to nops(pa) do p[1][j]:=pa[j] od: Q[1]:=[seq(p[1][j], j=1..nops(pa))]: for i from 2 to pa[1] do for j from 1 to F(Q[i-1]) do p[i][j]:=Q[i-1][j]-1 od: Q[i]:=[seq(p[i][j], j=1..F(Q[i-1]))] od: for i from 1 to pa[1] do col[i]:=[seq(Q[i][j]+nops(Q[i])-j, j=1..nops(Q[i]))] od: a:=proc(i, j) if i<=nops(Q[j]) and j<=pa[1] then Q[j][i]+nops(Q[j])-i else 1 fi end: A:=matrix(nops(pa), pa[1], a): product(product(A[m, n], n=1..pa[1]), m=1..nops(pa)); end: with(combinat): rev:=proc(a) [seq(a[nops(a)+1-i], i=1..nops(a))] end: sort([seq(H(rev(partition(8)[q])), q=1..numbpart(8))]);

Mathematica

h[l_] := With[{n = Length[l]}, Total[l]!/Product[Product[1 + l[[i]] - j + Sum[If[l[[k]] >= j, 1, 0], {k, i + 1, n}], {j, 1, l[[i]]}], {i, 1, n}]];
g[n_, i_, l_] := If[n == 0 || i == 1, h[Join[l, Array[1&, n]]], If[i < 1, 0, Flatten@Table[g[n - i*j, i - 1, Join[l, Array[i&, j]]], {j, 0, n/i}]]];
T[n_] := g[n, n, {}];
Sort[8!/T[8]] (* Jean-François Alcover, Aug 12 2024, after Alois P. Heinz in A060240 *)

Crossrefs

Cf. A003872, A060240.

Row n=8 of A093784.

Sequence in context: A147628 A153486 A062040 * A251646 A235177 A346514

Adjacent sequences: A093783 A093784 A093785 * A093787 A093788 A093789

Keyword

fini,full,nonn

Author

Emeric Deutsch, May 17 2004

Status

approved