%I#26 2021年1月27日10:09:38

%S 1,3,1,1,7,3,3,1,15,7,5,3,2,2,1,1,1,31,15,11,7,5,1,4,3,3，2,2，

%T 1,1,63,31,23,15,13,11,9,7,7,6,5,5,4,4,3,2,1,1,1127,63,43,31，

%U 25,21,19,15,14,12,11,10,9,9,8,7,7,6,5,5,4,3,2,2,1,1,1

%元素和可被k整除的[N]的非空子集的N个T（N，k）；三角形T（n，k），n>=1，1<=k<=n*（n+1）/2，按行读取。

%C T（n，k）是为所有n>=0，k>=1定义的。三角形只包含正项。如果k>n*（n+1）/2，T（n，k）=0。

%H Alois P.Heinz，<a href=“/A309402/b309402.txt”>行n=1..50，扁平</a>

%F和{k=1..n*（n+1）/2}k*T（n，k）=A309281（n）。

%对于n>=0，F T（n+1，n*（n+1）/2+1）=A000009（n）。

%e三角形T（n，k）开始：

%e 1；

%e 3、1、1；

%e第7、3、3、1、1和1条；

%e第15、7、5、3、3、2、2、1、1和1条；

%e 31、15、11、7、7、5、4、3、3、3、2、2、1、1、1；

%e第63、31、23、15、13、11、9、7、7、6、5、5、4、4、3、2、2、1、1和1条；

%e。。。

%p b:=proc（n，s）选项记忆`如果`（n=0，则添加（x^d，

%p d=数值[除数]（s），b（n-1，s）+b（n-l，s+n））

%p端：

%pT:=n->（p->seq（系数（p，x，i），i=1..度（p）））（b（n，0））：

%p序列（T（n），n=1..10）；

%t b[n_，s_]：=b[n，s]=如果[n==0，总和[x^d，

%t{d，除数[s]}]，b[n-1，s]+b[n-1，s+n]]；

%t t[n_]：=与[{p=b[n，0]}，表[系数[p，x，i]，

%t{i，1，指数[p，x]}]]；

%t阵列[t，10]//扁平（*_Jean-François Alcover_，2021年1月27日，以_Alois P.Heinz_*命名）

%Y列k=1表示A000225。

%Y行总和表示A309403。

%Y行长度给出A000217。

%Y T（n，n）给出A082550。

%Y行反向收敛到A000009。

%Y参考A309280和A309281。

%K nonn看tabf

%O 1,2号机组

%A _Alois P.Heinz，2019年7月28日