OFFSET
0,6
LINKS
Alois P. Heinz, Antidiagonals n = 0..52
FORMULA
G.f. of column k: Product_{j>=1} 1/(1-x^j)^(j^(k*j)).
EXAMPLE
Square array begins:
1, 1, 1, 1, ...
1, 1, 1, 1, ...
2, 5, 17, 65, ...
3, 32, 746, 19748, ...
5, 298, 66418, 16799044, ...
MAPLE
with(numtheory):
A:= proc(n, k) option remember; `if`(n=0, 1, add(add(
d*d^(k*d), d=divisors(j))*A(n-j, k), j=1..n)/n)
end:
seq(seq(A(n, d-n), n=0..d), d=0..10); # Alois P. Heinz, Mar 15 2017
MATHEMATICA
A[n_, k_] := If[n==0, 1, Sum[Sum[d*d^(k*d), {d, Divisors[j]}] *A[n - j, k], {j, n}] / n]; Flatten[Table[A[d - n, n], {d, 0, 10}, {n, d, 0, -1}]] (* Indranil Ghosh, Mar 17 2017 *)
PROG
(PARI) A(n, k) = if(n==0, 1, sum(j=1, n, sumdiv(j, d, d*d^(k*d)) * A(n - j, k))/n);
{for(d=0, 10, for(n=0, d, print1(A(n, d - n), ", "); ); print(); ); } \\ Indranil Ghosh, Mar 17 2017
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Seiichi Manyama, Mar 14 2017
STATUS
approved