OFFSET
0,2
COMMENTS
Equals column 1 of table A125781. Equals row sums and column 0 (shifted) of triangle A127420. - Paul D. Hanna, Feb 11 2007
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..300
PROG
(PARI) \\ Using Matrix Products
{A091351(n, k) = my(M=matrix(2, 2, r, c, if(r>=c, 1)), D, L); for(i=2, n,
\\ D = M shift down 1 row
D = matrix(#M+1, #M+1, r, c, if(r==1, 0, if(c>#M, 0, M[r-1, c])) );
\\ L = New matrix, shifted left 1 column
L = D + D^2;
\\ New matrix = L shifted right 1 column
M = matrix(#M+1, #M+1, r, c, if(c==1, 1, L[r, c-1]) ));
TM = M; M[n+1, k+1]}
\\ Print up to n terms of column 1 of triangle A091351
{upto(n) = A091351(n+1, 1); for(r=0, n, print1(TM[r+2, 2], ", "))}
upto(25) \\ Paul D. Hanna, Jul 15 2026
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Paul D. Hanna, Jan 02 2004
STATUS
approved
