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.

A103241 Unreduced numerators of the elements T(n,k)/(n-k)!, read by rows, of the triangular matrix P^-1, which is the inverse of the matrix defined by P(n,k) = (-k^3)^(n-k)/(n-k)! for n >= k >= 1.

This page as a plain text file.
%I A103241 #21 Jul 21 2024 00:27:56
%S A103241 1,1,1,15,8,1,1024,368,27,1,198581,53672,2727,64,1,85102056,18417792,
%T A103241 710532,11904,125,1,68999174203,12448430408,386023509,4975936,38375,
%U A103241 216,1,95264160938080,14734002979456,381535651512,3977848832,23945000
%N A103241 Unreduced numerators of the elements T(n,k)/(n-k)!, read by rows, of the triangular matrix P^-1, which is the inverse of the matrix defined by P(n,k) = (-k^3)^(n-k)/(n-k)! for n >= k >= 1.
%C A103241 Define a triangular matrix P where P(n,k) = (-k^3)^(n-k)/(n-k)!, then M = P*D*P^-1 = A102098 satisfies M^3 = SHIFTUP(M) where D is the diagonal matrix consisting of {1,2,3,...}. The operation SHIFTUP(M) shifts each column of M up 1 row.
%C A103241 Essentially equal to square array A082170 as a triangular matrix. The first column is A082162 (enumerates acyclic automata with 3 inputs).
%F A103241 For n > k >= 1: 0 = Sum_{m=k..n} C(n-k, m-k)*(-m^3)^(n-m)*T(m, k).
%F A103241 For n > k >= 1: 0 = Sum_{j=k..n} C(n-k, j-k)*(-k^3)^(j-k)*T(n, j).
%e A103241 Rows of unreduced fractions T(n,k)/(n-k)! begin:
%e A103241   [1/0!],
%e A103241   [1/1!, 1/0!],
%e A103241   [15/2!, 8/1!, 1/0!],
%e A103241   [1024/3!, 368/2!, 27/1!, 1/0!],
%e A103241   [198581/4!, 53672/3!, 2727/2!, 64/1!, 1/0!],
%e A103241   [85102056/5!, 18417792/4!, 710532/3!, 11904/2!, 125/1!, 1/0!], ...
%e A103241 forming the inverse of matrix P where P(n,k) = A103246(n,k)/(n-k)!:
%e A103241   [1/0!],
%e A103241   [-1/1!, 1/0!],
%e A103241   [1/2!, -8/1!, 1/0!],
%e A103241   [-1/3!, 64/2!, -27/1!, 1/0!],
%e A103241   [1/4!, -512/3!, 729/2!, -64/1!, 1/0!], ...
%o A103241 (PARI) {T(n,k)=my(P);if(n>=k&k>=1, P=matrix(n,n,r,c,if(r>=c,(-c^3)^(r-c)/(r-c)!))); return(if(n<k||k<1,0,(P^-1)[n,k]*(n-k)!))}
%Y A103241 Cf. A103246, A102098, A082170, A082162.
%K A103241 nonn,tabl,frac
%O A103241 1,4
%A A103241 _Paul D. Hanna_, Feb 02 2005