%I#32 2021年6月7日18:50:07

%S 0,1,2,7,331911297100100638766984701589893011039967031303132269，

%电话：17589153719254509227541393115823873564573130459613，

%电话：1124144767682215206776648944129654007606953861946878163539437728923181

%N三角形A104980的第1列；也等于三角形A104986的第0列，它等于A104980的矩阵对数。

%H G.C.Greubel，n表，n=0..440的a（n）</a>

%F From _Gary W.Adamson_，2011年7月14日：（开始）

%F设M=三角形A128175作为无限平方生产矩阵（删除第一个“1”）：

%F 1，1，0，0，0。。。

%F 2、2、1、0、0。。。

%F 4、4、3、1、0。。。

%F 8、8、7、4、1。。。

%F。。。

%F a（n）=M^（n-1）顶行项之和。示例：M^4的顶行=（71，71，38，10，1），总和=191=a（5）。（结束）

%F a（0）=1，a（n）=n*a（n-1）+和{j=1..n}A003319（j）*a（n-j），项1的偏移量为0_F.Chapoton_，2018年2月26日

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

%t a[n_]：=t[n，1]；

%t表[a[n]，{n，0，20}]（*Jean-François Alcover_，2018年8月9日*）

%o（PARI）{a（n）=如果（n<0，0，（矩阵（n+2，n+2、m、j），如果（m==j，1，如果（m==j+1，-m+1，-polcoeff（（1-1/总和（i=0，m，i！*x^i））/x+o（x^m），m-j-1）））^-1）[n+1,2]）}

%o（鼠尾草）

%o缓存函数

%o定义T（n，k）：

%o如果（k<0或k>n）：返回0

%o elif（k==n）：返回1

%o elif（k==n-1）：返回n

%o其他：对于（0..n-k-1）中的j，返回k*T（n，k+1）+总和（T（j，0）*T（n，j+k+1））

%o[T（n，1）for n in（0..30）]#_G.C.Greubel_，2021年6月7日

%Y参见A104980、A104982、A104886、A128175。

%K nonn公司

%0、3

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