|
Search: id:A093769
|
|
|
| A093769 |
|
Hook products of all partitions of 7. |
|
+0
2
|
|
| 144, 144, 240, 240, 252, 336, 336, 360, 360, 360, 360, 840, 840, 5040, 5040
(list; graph; listen)
|
|
|
OFFSET
|
1,1
|
|
|
COMMENT
|
a(n)=7!/A003871(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(7)[q])), q=1..numbpart(7))]);
|
|
CROSSREFS
|
Cf. A003871.
Sequence in context: A057404 A093159 A056628 this_sequence A124512 A101936 A044868
Adjacent sequences: A093766 A093767 A093768 this_sequence A093770 A093771 A093772
|
|
KEYWORD
|
fini,full,nonn
|
|
AUTHOR
|
Emeric Deutsch (deutsch(AT)duke.poly.edu), May 17 2004
|
|
|
Search completed in 0.002 seconds
|
