The On-Line Encyclopedia of Integer Sequences (OEIS)

Revision History for A113373

(Bold, blue-underlined text is an addition; faded, red-underlined text is a ~~deletion~~.)
Showing all changes.

A113373

Column 3 of triangle A113370, also equals column 0 of A113370^10.
(history; published version)

#5 by Charles R Greathouse IV at Tue Jun 13 23:21:07 EDT 2017

STATUS

editing

approved

#4 by Charles R Greathouse IV at Tue Jun 13 23:20:58 EDT 2017

PROG

(PARI) {a(n)=local(A, B); A=Mat(1); for(m=2, n+4, 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); A[n+4, 4]}

STATUS

approved

editing

#3 by Russ Cox at Fri Mar 30 18:36:52 EDT 2012

AUTHOR

_Paul D. Hanna (pauldhanna(AT)juno.com), _, Nov 13 2005

Discussion

Fri Mar 30

18:36

OEIS Server: https://oeis.org/edit/global/213

#2 by N. J. A. Sloane at Sat Nov 10 03:00:00 EST 2007

KEYWORD

nonn,new

nonn

AUTHOR

Paul D . Hanna (pauldhanna(AT)juno.com), Nov 13 2005

#1 by N. J. A. Sloane at Tue Jan 24 03:00:00 EST 2006

NAME

Column 3 of triangle A113370, also equals column 0 of A113370^10.

DATA

1, 10, 190, 4945, 163705, 6617605, 317416204, 17677189426, 1123660618048, 80411135580587, 6405386721320372, 562632770359355432, 54061761393693393203, 5643864294810397451552, 636388356169482970183598

OFFSET

0,2

COMMENTS

Column k of A113370 = column 0 of A113370^(3*k+1) for k>=0.

PROG

(PARI) {a(n)=local(A, B); A=Mat(1); for(m=2, n+4, 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); A[n+4, 4]}

CROSSREFS

Cf. A113370, A113371 (column 1), A113372 (column 2).

KEYWORD

nonn

AUTHOR

Paul D Hanna (pauldhanna(AT)juno.com), Nov 13 2005

STATUS

approved