OFFSET
0,9
COMMENTS
A(n,k) is also the number of k-partite partitions of (n)^k into at most n k-tuples. A(2,2) = 5: [(2,2)], [(2,1),(0,1)], [(2,0),(0,2)], [(1,2),(1,0)], [(1,1),(1,1)].
EXAMPLE
A(2,2) = 5: (2*3)^2 = 36 has 5 factorizations into at most 2 factors: 36, 2*18, 3*12, 4*9, 6*6.
Square array A(n,k) begins:
1, 1, 1, 1, 1, 1, ...
1, 1, 1, 1, 1, 1, ...
1, 2, 5, 14, 41, 122, ...
1, 3, 19, 171, 1675, 16683, ...
1, 5, 74, 1975, 64182, 2203215, ...
1, 7, 248, 20096, 2213016, 268446852, ...
MATHEMATICA
b[n_, k_, i_] := b[n, k, i] = If[n>k, 0, 1] + If[PrimeQ[n] || i<2, 0, Sum[ If[d > k, 0, b[n/d, d, i-1]], {d, Divisors[n][[2 ;; -2]]}]]; A[0, _] = 1; A[1, _] = 1; A[_, 0] = 1; A[n_, k_] := With[{t = Times @@ Prime[ Range[k] ]}, b[t^n, t^n, n]]; Table[diag = Table[A[n-k, k], {k, n, 0, -1}]; Print[ diag]; diag, {n, 0, 10}] // Flatten (* Jean-François Alcover, Jan 08 2016, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Alois P. Heinz, Mar 27 2015
STATUS
approved