OFFSET
1,13
LINKS
Alois P. Heinz, Antidiagonals n = 1..141, flattened
Eric Weisstein's World of Mathematics, Power Tower
Wikipedia, Knuth's up-arrow notation
Wikipedia, Tetration
FORMULA
EXAMPLE
Square array A(n,k) begins:
1, 1, 1, 1, 1, 1, 1, 1, ...
0, 1, 1, 1, 1, 1, 1, 1, ...
0, 1, 3, 3, 3, 3, 3, 3, ...
0, 2, 8, 14, 14, 14, 14, 14, ...
0, 2, 36, 72, 96, 96, 96, 96, ...
0, 9, 159, 489, 729, 849, 849, 849, ...
0, -6, 932, 3722, 6842, 8642, 9362, 9362, ...
0, 118, 5627, 33641, 71861, 102941, 118061, 123101, ...
MAPLE
f:= proc(n) f(n):= `if`(n=0, 1, (x+1)^f(n-1)) end:
A:= (n, k)-> (n-1)!*coeff(series(f(k), x, n+1), x, n):
seq(seq(A(n, 1+d-n), n=1..d), d=1..14);
# second Maple program:
b:= proc(n, k) option remember; `if`(n=0, 1, `if`(k=0, 0,
-add(binomial(n-1, j)*b(j, k)*add(binomial(n-j, i)*
(-1)^i*b(n-j-i, k-1)*(i-1)!, i=1..n-j), j=0..n-1)))
end:
A:= (n, k)-> b(n, min(k, n))/n:
seq(seq(A(n, 1+d-n), n=1..d), d=1..14);
MATHEMATICA
b[n_, k_] := b[n, k] = If[n == 0, 1, If[k == 0, 0, -Sum[Binomial[n - 1, j]*b[j, k]*Sum[Binomial[n - j, i]*(-1)^i*b[n - j - i, k - 1]*(i - 1)!, {i, 1, n - j}], {j, 0, n - 1}]]];
A[n_, k_] := b[n, Min[k, n]]/n;
Table[A[n, 1 + d - n], {d, 1, 14}, {n, 1, d}] // Flatten (* Jean-François Alcover, May 25 2018, translated from 2nd Maple program *)
CROSSREFS
KEYWORD
sign,tabl
AUTHOR
Alois P. Heinz, Nov 12 2017
STATUS
approved