login
A167017
Column 1 of triangle T=A167015: a(n) = T(n+1,1) = [T^(n+1)](n-1,0) for n>0 with a(0)=1.
4
1, 1, 3, 14, 105, 1080, 13916, 212634, 3723198, 73047825, 1581256303, 37350588290, 954796686233, 26246234745486, 771853099673490, 24180166328032252, 803996438505301917, 28283254774663359546, 1049653570430645994556
OFFSET
0,3
COMMENTS
In triangle T=A167015, row n of T equals the initial n terms in column 0 of matrix power T^n, reversed and with an appended '1'.
PROG
(PARI) {A167015(n, k)=local(M, N); if(n==k||n==k+1, 1, if(n==k+2, k+2, N=matrix(n, n, r, c, if(r>=c, A167015(r-1, c-1))); M=matrix(n+1, n+1, r, c, if(r>=c, if(r<=n, N[r, c], (N^n)[n-k, 1]))); M[n+1, k+1]))}
{a(n)=A167015(n+1, 1)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 27 2009
STATUS
approved