login
A109283
Column 0 of triangle A109282.
4
1, 1, 3, 15, 96, 735, 6447, 63120, 677739, 7878921, 98216466, 1302951903, 18283530882, 270031107630, 4180037999178, 67581115214649, 1137710941143720, 19890998444408631, 360327514948457178, 6749464344626557092
OFFSET
0,3
COMMENTS
Triangular matrix T=A109282 satisfies: T(n,k) = [T^3](n-1,k) for n>k+1>=1, with T(n,n) = 1 and T(n+1,n) = n+1 for n>=0; also, T^(m+3) = SHIFT_UP(T^(m+1) - T^m) - D*T^m for all m where diagonal matrix D = [0,1,2,3,...] and SHIFT_UP shifts each column up 1 row.
PROG
(PARI) {a(n)=local(M=matrix(n+1, n+1)); M=M^0; for(i=1, n, M=matrix(n+1, n+1, r, c, if(r>=c, if(r==c, 1, if(r==c+1, c, (M^3)[r-1, c]))))); return(M[n+1, 1])}
CROSSREFS
Cf. A109153, A109282 (triangle), A109284 (column 1), A109285 (column 2), A109286 (row sums).
Sequence in context: A304072 A076301 A112913 * A370210 A079689 A262751
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 24 2005
STATUS
approved