%I#9 2021年6月7日09:27:10

%S 1,2,1,8,4,1,42,20,61266120,38,8,11954836270,62,10,1162706616，

%电话：2150516,92,12,115121857619030468882128,14,155134573672，

%电话：1856744651690501392170,16,1174141146159976198231050232499994159562070218,18

%N三角形A104980的矩阵正方形。

%C三角矩阵A104980满足：当p>=0时，SHIFT_LEFT（A104980^p的第0列）=p*（A104990的第p+1列）。

%H G.C.Greubel，<a href=“/A104988/b104988.txt”>三角形的n=0..50行，展平</a>

%F T（n+1，0）=2*A104980（n+3，3）=2*A104982（n）对于n>=0。

%e三角形开始：

%e 1；

%e 2，1；

%e 8、4、1；

%e第42、20、6、1条；

%e第266、120、38、8、1页；

%e 1954、836、270、62、10、1；

%e 16270、6616、2150、516、92、12、1；

%e 151218、58576、19030、4688、882、128、14、1；

%e 1551334、573672、185674、46516、9050、1392、170、16、1；

%e 17414114、6159976、1982310、502324、99994、15956、2070、218、18、1；

%t t[n_，k_]：=t[n，k]=如果[k<0|k>n，0，如果[k==n，1，如果[k==n-1，n，k*t[n、k+1]+和[t[j，0]*t[n，j+k+1]，{j，0，n-k-1}]]；（*t=A104980*）

%t M：=使用[{q=20}，表[If[j>i，0，t[i，j]]，{i，0、q}，{j，0，q}]]；

%t表[MatrixPower[M，2][[n+1，k+1]]，{n，0,12}，{k，0，n}]//扁平（*_G.C.Greubel_，2021年6月7日*）

%o（PARI）T（n，k）=如果（n<k||k<0，0，（矩阵（n+1，n+1，m，j，if（m==j，1，if）（m==j+1，-m+1，-polcoeff（（1-1/总和（i=0，m，i！*x^i））/x+o（x^m），m-j-1）））^-2）[n+1，k+1]

%Y参见A104980、A104982（第0列）、A104899（行总和）。

%K nonn，表

%0、2

%A·保罗·D·汉纳，2005年4月10日