OFFSET
0,2
COMMENTS
PROG
(PARI) {a(n)=local(A=matrix(2, 2), B); A[1, 1]=1; for (m=2, n+1, B=matrix(m, m); for (i=1, m, for (j=1, i, if(j==i, B[i, j]=j, if(j==1, B[i, j]=(A^3)[i-1, 1], B[i, j]=(A^3)[i-1, j])); )); A=B); return((A^2)[n+1, 2])}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 21 2005
STATUS
approved