登录
方阵转置A285321型.

%I#8 2019年12月7日12:18:29

%S 2,4,3,8,9,6,16,27,12,5,32,81,18,25,10,64243,2425,20,15128729,36,

%电话625,40,45,302562187,483125,50,75,60,75126561,5415625,80135,90,

%U 49,14102419683,7278125100225120343,28,21号

%N方阵A285321的转置。

%C见A285321。

%o(方案)(定义(A285322 n)(A285312bi(A004736 n)(A002260 n));;对于A285321bi,请参见A285321。

%o(Python)

%o来自操作员导入mul

%o来自sympy import prime,primefactors

%o def a019565(n):如果n>0,则返回reduce(mul,(prime(i+1)for i,v in enumerate(bin(n)[:1:-1])if v=='1'))_

%o def a007947(n):如果n<2,则返回1,否则减少(mul,素数(n))

%o定义a065642(n):

%o如果n==1:返回1

%o r=a007947(n)

%o n=n+r

%哦,虽然是007947(n)=回复:

%o n+=r

%o返回n

%o定义A(n,k):如果n==1,则返回a019565(k),否则返回a065642(A(n-1,k))

%o表示范围(1,11)中的n:打印[A(n-k+1,k)表示范围(1,n+1)中的k]#_Indranil Ghosh_,2017年4月19日

%Y转座:A285321。

%K nonn,表

%O 1,1号机组

%A _Antti Karttunen,2017年4月17日