%I#19 2018年12月10日08:52:56

%S 1,1,1,1,3,1,4,7,1,1,7,13,15,1,1,6,35,40,31,1,12,31155121,63,1，

%电话：1,8,91156651364127,1,15,5760078126671093255,1,13155，

%电话：40037513906107953280511,1,1813013952801229321953143598411023,1

%反对角线读取的N数组：T（N，k）是一般k维格中索引N的子格数（N>=1，k>=1）。

%D Günter Scheja，Uwe Storch，代数Lehrbuch，第2版。BG Teubner，斯图加特，1988年。【§63，Aufg.13】

%Hálvar Ibeas，<a href=“/A160870/b160870.txt”>前100名反糖尿病患者，扁平</a>

%H Michael Baake，<a href=“http://arxiv.org/abs/math/0605222“>维度d≤4的重合问题的求解，arXiv:math/0605222[math.MG]，2006。【附录A】

%H B.Gruber，<a href=“https://doi.org/10.107/S0108767397009781“>子晶格数的替代公式</a>，《晶体学报》A53（1997）807-808。

%H J.H.Kwak和J.Lee，<a href=“https://doi.org/10.1142/9789812799890_0005“>图覆盖、表面分支覆盖和相关群论的枚举，摘自组合数学和计算数学（Pohang，2000），S.Hong等人，《世界科学》，新加坡，2001年，第97-161页。

%H Yi Ming Zou，<a href=“https://arxiv.org/abs/math/0610684“>高斯二项式和子格数，arXiv:math/0610684[math.CO]，2006。

%H Yi Ming Zou，<a href=“https://doi.org/10.1107/S01076730602455X“>高斯二项式和子格数，《晶体学报》A62（2006）409-410。

%F T（n，1）=1；T（1，k）=1；T（n，k）=和{d|n}d*T（d，k-1）。

%F自2015年10月31日起：（开始）

%F T（n，k）=和{d|n}（n/d）^（k-1）*T（d，k-1）。

%F T（乘积（p^e），k）=乘积（Gaussian_poly[e+k-1，e]_p）。

%F（结束）

%e阵列开始：

%e 1,1,1,1,1,1,1,1,1,1,1,1,,1,1,11,1,1,1_1,1,1,1，1,1,1,1-1,1,1,1，1,1,1-1，。。。

%电子邮箱：1,3,7,15,312,63127255511023204740958191163833276765535，。。。

%电子邮箱：1,4,13,40121364109332809841295248857573265727207971612391484，。。。

%电子邮箱：1,7,35155651266710795434517425169802727941551180715，。。。

%电子邮箱：1,6,311567813906195319765648828124414061220703161035156，。。。

%e。。。

%tT[_，1]=1；T[1，_]=1；T[n_，k_]：=T[n，k]=除数和[n，（n/#）^（k-1）*T[#，k-1]&]；表[T[n-k+1，k]，{n，1，11}，{k，1，n}]//Flatten（*_Jean-François Alcover_，2015年12月4日*）

%o（PARI）：

%o阿杜（M）=

%o｛/*由AntiDiagonals读取，向上*/

%o局部（N=材料尺寸（M）[1]）；

%o表示（n=1，n，表示（j=0，n-1，打印1（M[n-j，j+1]，“，”））；

%o}（o）

%o T（n，k）=

%o（o）{

%o如果（（n==1）||（k==1，返回（1））；

%o返回（sumdiv（n，d，d*T（d，k-1）））；

%o}（o）

%o M=矩阵（15,15，n，k，T（n，k））/*方阵*/

%o adu（M）/*序列*/

%Y列：A000203、A001001、A038991、A0389962、A03899.3、A0389994、A03899、A03896、A03897。

%Y行：A000012、A000225、A003462、A006095、A00346、A160869、A023000、A00609。

%Y换位数组：A128119。

%K non，tabl，简单

%O 1,5型

%A _N.J.A.Sloane，2009年11月19日