cp's OEIS Frontend

A104446 Square of triangular matrix A104445, read by rows, where X=A104445 satisfies: SHIFT_LEFT_UP(X) = X^2 - X + I.

%I A104446 #7 Jun 13 2017 01:29:48
%S A104446 1,2,1,3,2,1,5,5,2,1,10,13,7,2,1,25,39,25,9,2,1,78,139,100,41,11,2,1,
%T A104446 296,587,459,205,61,13,2,1,1330,2897,2418,1149,366,85,15,2,1,6935,
%U A104446 16462,14506,7233,2421,595,113,17,2,1,41352,106301,98161,50905,17706,4535,904
%N A104446 Square of triangular matrix A104445, read by rows, where X=A104445 satisfies: SHIFT_LEFT_UP(X) = X^2 - X + I.
%C A104446 Column 0: T(n,0) = 1 + A091352(n-1) for n>0. Column 1 is A104447. Row sums form A104448.
%F A104446 T(n, k) = A104445(n, k) + A104445(n+1, k+1) - I(n, k), where I=identity matrix. T(n, k) = A091351(n-1, k) + A091351(n, k+1) - I(n, k), for n>k>=0.
%e A104446 Rows begin:
%e A104446 1;
%e A104446 2,1;
%e A104446 3,2,1;
%e A104446 5,5,2,1;
%e A104446 10,13,7,2,1;
%e A104446 25,39,25,9,2,1;
%e A104446 78,139,100,41,11,2,1;
%e A104446 296,587,459,205,61,13,2,1;
%e A104446 1330,2897,2418,1149,366,85,15,2,1
%e A104446 6935,16462,14506,7233,2421,595,113,17,2,1; ...
%o A104446 (PARI) T(n,k)=local(A=Mat(1),B);for(m=1,n,B=A^2-A+A^0; A=matrix(m+1,m+1);for(i=1,m+1, for(j=1,i, if(i<2 || j==i,A[i,j]=1,if(j==1,A[i,j]=1,A[i,j]=B[i-1,j-1]))))); return((A^2)[n+1,k+1])
%Y A104446 Cf. A104445, A091351, A091352, A104447, A104448.
%K A104446 nonn,tabl
%O A104446 0,2
%A A104446 _Paul D. Hanna_, Mar 08 2005