cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A113365 Matrix cube of triangle A113350.

This page as a plain text file.
%I A113365 #5 Jun 13 2017 23:21:20
%S A113365 1,6,1,39,12,1,327,138,18,1,3556,1830,297,24,1,48659,28805,5349,516,
%T A113365 30,1,812462,535004,109095,11724,795,36,1,16136404,11568197,2529909,
%U A113365 292894,21795,1134,42,1,373415239,287143993,66345668,8117624,643790,36402,1533
%N A113365 Matrix cube of triangle A113350.
%F A113365 Column k of A113350^3 = column 1 of A113340^(2*k+2) for k>=0.
%e A113365 Triangle begins:
%e A113365 1;
%e A113365 6,1;
%e A113365 39,12,1;
%e A113365 327,138,18,1;
%e A113365 3556,1830,297,24,1;
%e A113365 48659,28805,5349,516,30,1;
%e A113365 812462,535004,109095,11724,795,36,1;
%e A113365 16136404,11568197,2529909,292894,21795,1134,42,1;
%e A113365 373415239,287143993,66345668,8117624,643790,36402,1533,48,1; ...
%o A113365 (PARI) T(n,k)=local(A,B);A=matrix(1,1);A[1,1]=1;for(m=2,n+1,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^(2*j-1))[i-j+1,1]));));A=B); (matrix(#A,#A,r,c,if(r>=c,(A^(2*c))[r-c+1,1]))^3)[n+1,k+1]
%Y A113365 Cf. A113340, A113350, A113346 (column 0), A113366 (column 1), A113367 (column 2); A113355 (=A113350^2), A113360 (=A113340^3).
%K A113365 nonn,tabl
%O A113365 0,2
%A A113365 _Paul D. Hanna_, Nov 09 2005