OFFSET
5,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 5..750
MAPLE
b:= proc(n, k, t) option remember; `if`(t=0, 1, `if`(t=1,
A(n-1, k), add(A(j, k)*b(n-j-1, k, t-1), j=0..n-2)))
end:
A:= proc(n, k) option remember; `if`(n=0, 1,
-add(binomial(k, j)*(-1)^j*b(n+1, k, 2^j), j=1..k))
end:
T:= proc(n, k) option remember;
add(A(n, k-i)*(-1)^i*binomial(k, i), i=0..k)
end:
a:= n-> T(n, 5):
seq(a(n), n=5..25);
MATHEMATICA
b[n_, k_, t_] := b[n, k, t] = If[t == 0, 1, If[t == 1, A[n - 1, k], Sum[A[j, k]*b[n - j - 1, k, t - 1], {j, 0, n - 2}]]];
A[n_, k_] := A[n, k] = If[n == 0, 1, -Sum[Binomial[k, j]*(-1)^j*b[n + 1, k, 2^j], {j, 1, k}]];
T[n_, k_] := Sum[A[n, k - i]*(-1)^i*Binomial[k, i], {i, 0, k}];
a[n_] := T[n, 5];
a /@ Range[5, 25] (* Jean-François Alcover, Dec 11 2020, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, May 29 2015
STATUS
approved