OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..500
MATHEMATICA
T[n_, k_]:= T[n, k]= If[k==0, 1, If[k==n, 0, T[n-1, k] - 3*T[n-1, k-1] + 3*T[n, k-1] + T[n+1, k-1] ]];
Table[T[n, n-3], {n, 3, 30}] (* G. C. Greubel, Feb 18 2021 *)
PROG
(PARI) my(x='x+O('x^33)); Vec((serreverse(x*(1-3*x+sqrt((1-3*x)*(1-7*x)))/2/(1-3*x))/x)^3)
(Sage)
@CachedFunction
def T(n, k):
if (k<0 or n<k): return 0
elif (k==0): return 1
elif (k==n): return 0
else: return T(n-1, k) - 3*T(n-1, k-1) + 3*T(n, k-1) + T(n+1, k-1)
[T(n, n-3) for n in (3..30)] # G. C. Greubel, Feb 18 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Apr 26 2006
STATUS
approved