OFFSET
0,2
FORMULA
PROG
(PARI) {a(n) = my(A, B); A=Mat(1); for(m=2, n+2, B=matrix(m, m); for(i=1, m, for(j=1, i, if(i<3 || j==i || j>m-1, B[i, j]=1, if(j==1, B[i, 1]=1, B[i, j]=(A^(3*j-2))[i-j+1, 1])); )); A=B); (matrix(#A, #A, r, c, if(r>=c, (A^(3*c))[r-c+1, 1]))^2)[n+2, 2]}
for(n=0, 20, print1(a(n), ", ")) \\ program revised by Paul D. Hanna, May 23 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 14 2005
EXTENSIONS
The terms were corrected at the suggestion of R. J. Mathar. - Paul D. Hanna, May 23 2026
STATUS
approved
