OFFSET
0,3
COMMENTS
A k-packed matrix of size n X n is a matrix with entries in the alphabet A_k = {0,1, ..., k} such that each row and each column contains at least one nonzero entry.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..662 (rows 0..12 flattened)
H. Cheballah, S. Giraudo, R. Maurice, Combinatorial Hopf algebra structure on packed square matrices, arXiv preprint arXiv:1306.6605 [math.CO], 2013.
FORMULA
From Andrew Howroyd, Sep 20 2017: (Start)
T(n, k) = Sum_{i=0..n} Sum_{j=0..n} (-1)^(i+j) * binomial(n,i) * binomial(n,j) * binomial(i*j,k) * 2^k.
T(n, k) = 0 for n > k.
T(n, n) = A000165(n).
(End)
EXAMPLE
Triangle begins:
1
0 2
0 0 8 32 16
0 0 0 48 720 2880 4992 4608 2304 512
...
MATHEMATICA
p[k_, n_, l_] := Sum[(-1)^(i+j)*Binomial[n, i]*Binomial[n, j]*Binomial[i*j, l]*k^l, {i, 0, n}, {j, 0, n}];
T[n_, k_] := p[2, n, k];
Table[T[n, k], {n, 0, 5}, {k, 0, n^2}] // Flatten (* Jean-François Alcover, Oct 08 2017, translated from PARI *)
PROG
(PARI) \\ T(n, k) = p(2, n, k) (see Cheballah et al. ref).
p(k, n, l) = {sum(i=0, n, sum(j=0, n, (-1)^(i+j) * binomial(n, i) * binomial(n, j) * binomial(i*j, l) * k^l))}
for (n=0, 5, for(k=0, n^2, print1(p(2, n, k), ", ")); print); \\ Andrew Howroyd, Sep 20 2017
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
N. J. A. Sloane, Nov 09 2013
EXTENSIONS
Terms a(18) and beyond from Andrew Howroyd, Sep 20 2017
STATUS
approved