%I#25 2021年5月7日10:42:25
%S 2,6210357357064845819350301990,
%电话286956625737391526978466144187168109530,
%电话:1038300112150956151877699324649731518883355380534272386781875587619359740733888844803014212990
%N a(0)=2,对于N>0,a(N)=a(N-1)*A242378(N,a(N-1。
%C不规则表格A053632的乘法编码(A007188和A260443样式)。
%H Indranil Ghosh,n的表格,n=0..9的a(n)</a>
%F a(0)=2,对于n>0,a(n)=a(n-1)*A242378(n,a(n-1))。
%F其他身份。对于所有n>=0:
%F A001222(a(n))=A000079(n)。
%F A048675(a(n))=A028362(1+n)。
%F A248663(a(n))=A068052(n)。
%F A007913(a(n))=A285102(n)。
%o(PARI)
%o A003961(n)={my(f=因子(n));对于(i=1,#f~,f[i,1]=下一素数(f[i、1]+1));因子回退(f);};
%o A242378(k,n)={而(k>0,n=A003961(n);k=k-1);n;};
%o A285101(n)={if(0==n,2,A285101(n-1)*A242378(n,A285101(n-1));};
%o(方案)(定义(A285101n)(如果(零?n)2(*(A285101(-n 1)))(A242378bi n(A285109(-n))));;关于A242378bi,请参见A242378。
%o(Python)
%o来自sympy import factorint,prime,primepi
%o来自操作员导入mul
%o从functools导入reduce
%o定义a003961(n):
%o f=因子(n)
%o如果n==1,则返回1,否则减少(mul,[prime(primepi(i)+1)**f[i]for i in f])
%o定义a242378(k,n):
%o当k>0时:
%o n=a003961(n)
%o k-=1
%o返回n
%o l=[2]
%o表示范围(1,7)内的n:
%o x=l[n-1]
%o l.append(x*a242378(n,x))
%o打印(l)#_Indranil Ghosh,2017年6月27日
%Y参见A001222、A003961、A007913、A028362、A048675、A053632、A068052、A242378、A248663、A285102。
%Y另请参见A007188、A260443。
%K nonn公司
%0、1
%2017年4月15日,安蒂·卡图宁
|