%我#36 2021年3月27日22:47:39

%S 0,0,1,0,0,0，0,0'，0,1,2,0,00,0',0,2,0,1,0，0，0,1，0,0，

%T 0,3,0,0,0-0,0,1,0,0，0,0'，0,0，

%U 0,0,00,0,1,0,2,0,0.0,0,O,0,0

%N表T（N，k）由除以N的第k个素数的最大幂指数的反对偶得到。

%F T（n，k）=对数（A060176（n，k））/log（A000040（k））=A054841（n）右侧的第k位。

%e a（12,1）=2，因为4=2^2=p_1^2除以12，但8=2^3不等于。

%e a（12,2）=1，因为3=p_2除以12，但9=3^2不等于。

%e另请参见A249344中的示例，它是该数组的转置。

%e阵列的左上角：

%e 0，0，0。。。

%e 1，0，0，0,0，0。。。

%e 0，1，0，0，0,0，0。。。

%e 2，0，0，0,0，0。。。

%e 0，0，1，0，0。。。

%e。。。

%t t[n_，k_]：=整数指数[n，素数[k]]；

%t表[t[n-k+1，k]，{n，1,15}，{k，n，1，-1}]//Flatten（*Jean-François Alcover_，2019年11月18日*）

%o（方案）

%o（定义（A060175 n）（A249344bi（A004736 n）（A002260 n））

%o（定义（A249344双行列）（let（（p（A000040行）））（let loop（（n col）（i 0）））

%o_Antti Karttunen，2014年10月28日

%o（Python）

%o来自sympy import prime

%o定义a（n，k）：

%o p=素数（n）

%o i=z=0

%o而p**i<=k：

%o如果k%（p**i）==0:z=i

%o i+=1

%o返回z

%o对于范围（1，10）中的n：打印（[a（n-k+1，k）对于范围（1,n+1）中的k）]）#_Indranil Ghosh，2017年6月24日

%o（PARI）a（n，k）=估值（n，质数（k））；\\_米歇尔·马库斯，2017年6月24日

%Y转座：A249344。

%Y列1:A007814。

%Y列2:A007949。

%Y列3:A112765。

%Y第4列：A214411。

%Y参考A002260、A004736、A054841、A060176、A085604、A090622、A115627、A249421、A249422。

%K轻松，不，tabl

%O 1,10号

%2001年3月14日，安利底特律

%E Antti Karttune_于2014年10月28日纠正错误示例并计算更多术语