OFFSET
0,8
FORMULA
T(n, 0) = 1. T(n, 1) = a(n) if n>0, T(n, 2) = a(n+1)*a(n-1) - a(n)^2 if n>1, T(n, 3) = det([a(n-2), a(n-1), a(n); a(n-1), a(n), a(n+1); a(n), a(n+1), a(n+2)]) if n>2 where a(n) = A004148(n).
EXAMPLE
Table is:
n\k 0 1 2 3 4 5 6
-- -- -- -- -- -- -- --
0 | 1
1 | 1 1
2 | 1 1 1
3 | 1 2 0 0
4 | 1 4 0 0 -1
5 | 1 8 4 -2 -1 -1
6 | 1 17 7 3 -3 -1 -1
PROG
(PARI) {T(n, k) = my(m); if( k<0 || k>n, 0, matdet( matrix(k, k, i, j, polcoeff( (1 - x + x^2 - sqrt(1 - 2*x - x^2 + x^3*(-2 + x + O(x^(m=i+j+n-k-1))))) / (2*x^2), m))))};
CROSSREFS
KEYWORD
sign,tabl
AUTHOR
Michael Somos, Oct 04 2006
STATUS
approved