%I#12 2023年9月9日11:34:01
%S 1,2,2,3,2,4,2,3,3,4,2,6,2,44,4,4,1,6,6,4,2,6,3,6,2,8,2,4,1,4,4,
%温度4,9,2,4,4,6,2,8,2,6,6,4,2,8,3,6,6,1,6,5,4,4,12,2,4,5,8,2,6,
%U 4,8,2,9,2,4,6,6,4,8,1,8,4,4,2,12,4,4
%可被N整除的最小指数奇数的指数奇数除数。
%C First在n=32时与A049599和A282446不同,在n=64时与A353898不同。
%H Amiram Eldar,n的表,n=1..10000的a(n)</a>
%H Vaclav Kotesovec,图表-渐近比率(100000项)。
%F a(n)=A322483(A356191(n))。
%F与a的乘积(p^e)=上限(e+3)/2)。
%F Dirichlet g.F.:zeta(s)*zeta(2*s)*Product_{p prime}(1+1/p^s-1/p^(3*s))。
%F From _Vaclav Kotesovec_,2023年9月9日:(开始)
%设F(s)=Product_{pprime}(1-1/p^(2*s)-1/p^(3*s)+1/p^(4*s))。
%F Dirichlet g.F.:zeta(s)^2*zeta(2*s)*F(s)。
%F和{k=1..n}a(k)~(Pi^2*F(1)*n/6)*(log(n)+2*gamma-1+12*zeta'(2)/Pi^2+F'(1)/F(1)),其中
%F F(1)=乘积{p素数}(1-1/p^2-1/p^3+1/p^4)=0.535896153828337999808502631318545950648222374514145271510108346133288119。。。,
%F F’(1)=F(1)*Sum_{p prime}(-4+3*p+2*p^2)*log(p)/(1-p-p^2+p^4)=F(1)*1.4525924794445159559037143958285473414824651141192913672347667991。。。
%F和gamma是Euler-Marcheroni常数A001620。(结束)
%t f[p_,e_]:=上限[(e+3)/2];a[1]=1;a[n_]:=倍@@f@@FactorInteger[n];数组[a,100]
%o(PARI)a(n)=vecprod(应用(x->细胞((x+3)/2),系数(n)[,2]));
%Y参见A322483、A356191。
%Y参见A049599、A282446、A353898。
%K nonn,简单,多
%O 1,2号机组
%A阿米拉姆·埃尔达尔,2023年9月8日
|