%I#29 2024年2月7日01:16:30
%S 1,1,1,1,2,2,1,1,1,1,1,2,2,2,2,2,2,4,3,3,6,6,1,1,1,2,2,1,1,1,1,
%温度2,2,2,2,4,4,3,3,2,6,6,2,2,2,2,4,2,2,4,1,8,6,6,
%U 6,6,12,12,3,3,3-3,6,6,3,3A,3,6,16,66,6
%N阶乘基中N的非零位数的乘积。
%C a(0)=1,因为空乘积总是给出1。
%H Antti Karttunen,n的表,n=0..5040的a(n)</a>
%F对于所有n,a(A227157(n))=A208575(A227151(n。
%t a[n_]:=模[{k=n,m=2,r,p=1},而[{k,r}=商余数[k,m];k!=0||r!=0,如果[r>0,p*=r];m++];p] ;数组[a,100,0](*_Amiram Eldar_,2024年2月7日*)
%o(MIT/GNU方案)
%o(定义(A227153 n)(应用*(删除匹配项(n->factbase n)零?)))
%o(定义(n->factbase n)(let loop((n n)(fex(if(零?n)(列表0)(列表)))(i 2))(cond((零?n)fex)(else(loop(floor->exact(/n i))))
%o(Python)
%o来自操作员导入mul
%o定义A(n,p=2):
%o如果n<p其他A(n//p,p+1)*10+n%p,则返回n
%o定义a(n):
%o如果n<2则返回1,否则reduce(mul,[int(i)for i in str(A(n-1))if i!=“0”])
%o打印([a(n)代表范围(201)中的n)]#_Indranil Ghosh,2017年6月19日
%Y A227157给出了与A208575相等的位置。
%Y参见A007623、A227154、A227191。
%K nonn,基础
%0、5
%2013年7月4日,安蒂·卡图内
|