%I#20 2023年2月4日09:53:16

%S 1,1,6,1,10,27,1,14,69108,1,18123404405,1,2219589221551458,1，

%电话2627317165845108305103,1,30375273213525360425224117496,1，

%电话：344774324245759977421364724464859049、18603604454520814670521512325041120599830

%反对角线读取的表：T（N，k）是秩为1的幂等N X N 0..k矩阵的个数。

%C表启动

%C1,1,1,1,1，1,1,1,1,1。。。

%C 6、10、14、18、22、26、30、34、38。。。

%C 27、69、123、195、273、375、477、603。。。

%C 108、404、892、1716、2732、4324、6060。。。

%C 405、2155、5845、13525、24575、44545。。。

%C 1458、10830、36042、99774、208146。。。

%C 510352241213647705215。。。

%C 17496，244648。。。

%C 59049。。。

%C。。。

%H Robert Israel，n表，n=1..10011的a（n）</a>

%e n=3，k=4的一些解：

%e 1 0 0 0 4 4 0 0 0 0 2 1 2 1 0 0 00 1 0

%e 0 0 0 0 1 1 3 1 0 0 0 00 0 0 0零0 1 0

%e 10 0 0 0 0 00 0 0 0 2 1 0 0 0 1 4 1 0 0

%p f：=进程（n，k）

%p局部tot，a1，a0，a2，m，u；

%p总计：=0；

%a1从1到n的p do

%p表示a0从0到n-a1 do

%p a2：=n-a1-a0；

%p如果a0=0，则tot:=tot+n/（a1！*a2！）*a1*（k-1）^a2

%p elif a2=0，然后tot:=tot+n/（a0！*a1！）*a1*（k+1）^a0

%p其他

%p u：=n/（a0！*a1！*a2！）*a1；

%p代表从2到kdo的m

%p tot:=tot+u*（（m-1）^a2-（m-2）^a2）*（楼层（k/m）+1）^a0

%日期

%功率因数

%p od od；

%p总计

%p端程序：

%p序列（序列（f（i，j-i），i=1..j-1），j=2..20）；#_罗伯特·伊斯雷尔，2019年12月15日

%t取消[电源]保护；0^0=1；保护[电源]；

%t f[n_，k_]：=模[{tot，a1，a0，a2，m，u}，tot=0；对于[a1=1，a1<=n，a1++，对于[a0=0，a0<=n-a1，a0++，a2=n-a1-a0；其中[a0==0，tot=tot+n！/）*a1*（k+1）^a0，真，u=n！/（a0！*a1！*a2！）*a1；对于[m=2，m<=k，m++，tot=tot+u*（（m-1）^a2-（m-2）^a2）*（楼层[k/m]+1）^a0]]]]；tot]；

%t表[表[f[i，j-i]，{i，1，j-1}]，{j，2，20}]//扁平（*Jean-François Alcover_，2023年2月4日，在_Robert Israel_*之后）

%Y列1为A027471（n+1）。

%K nonn，表

%氧1,3

%A R.H.Hardin，2013年4月9日