%I#28 2023年10月23日16:15:32
%S 2,5,9,13,17,23,27,32,38,44,48,56,60,66,74,80,84,93,97105113119,
%电话12313313914515316116517718118819620221022226232240,
%电话:25025426627027829029630031231832733534347359367377385391395411415421433441449461465473481
%N A099777的部分和,其中A099776(N)给出了第N个偶数的除数。
%H Antti Karttunen,n表,n=1..10082的a(n)</a>
%H A.Karttunen,<A href=“https://oeis.org/plot2a?name1=A263086&;名称2=A263085&;tform1=未转换&;tform2=未转换&;移位=0&;radiop1=比率&;drawlines=true“>使用OEIS Plot2脚本绘制的比率a(n)/A263085(n)</a>
%F a(1)=2;对于n>1,a(n)=A000005(2*n)+a(n-1)[其中A000005。
%F其他身份。对于所有n>=1:
%F a(n)=A263084(n)+A263085(n)。
%F a(n)~n/2*(3*log(n)+log(2)+6*gamma-3),其中gamma是Euler-Marcheroni常数A001620.-_瓦茨拉夫·科特索维奇,2019年2月13日
%F From _Ridouane Oudra,2019年8月24日:(开始)
%F a(n)=和{k=1..n}A000005(2*k)
%F a(n)=A006218(n)+A060831(n)。(结束)
%p与(数字理论):seq(相加(tau(2*k),k=1..n),n=1..60);#_Ridouane Oudra,2019年8月24日
%t累加[DivisorSigma[0,2范围@69]](*_Michael De Vlieger_,2015年10月13日*)
%o(方案,带有备忘录-宏定义)
%o(定义(A263086 n)(如果(=1 n)(A099777 n)(+(A09977 n)(A2630086(-n 1))))
%o(PARI)a(n)=总和(k=1,n,numdiv(2*k));\\_米歇尔·马库斯,2019年8月25日
%o(Python)
%o来自数学导入isqrt
%o定义A263086(n):返回(t:=isqrt(m:=n>>1)
%Y参见A000005、A006218、A099777、A263084、A2630.85、A060831。
%Y另请参阅A262518、A262519。
%K nonn公司
%O 1,1号机组
%2015年10月12日,安蒂·卡图内
|