login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A109154
Column 1 of triangle A109152.
5
1, 2, 4, 14, 56, 262, 1348, 7574, 45616, 292442, 1978788, 14060654, 104403304, 807162070, 6476169012, 53781863598, 461194543136, 4075567775538, 37048022605956, 345883045698222, 3311798037556952, 32480125869403430
OFFSET
0,2
COMMENTS
Triangular matrix T=A109152 satisfies: T(n,k) = [T^2](n-1,k) for n>k+1>=1, with T(n,n) = 1 and T(n+1,n) = n+1 for n>=0.
PROG
(PARI) {a(n)=local(M=matrix(n+2, n+2)); M=M^0; for(i=1, n, M=matrix(n+2, n+2, r, c, if(r>=c, if(r==c, 1, if(r==c+1, c, (M^2)[r-1, c]))))); return(M[n+2, 2])}
CROSSREFS
Cf. A109152 (triangle), A109153 (column 0), A109155 (column 2), A109156 (row sums).
Sequence in context: A047990 A000939 A367558 * A030853 A306150 A030962
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 20 2005
STATUS
approved