%I#33 2024年1月8日01:37:26
%S 0,0,0,1,1,1,1,1,2,1,2,1,1,3,1,2,2,2,1,4,1,2,3,1,4,1,3,2,2,
%温度2,5,1,2,2,4,1,5,1,3,3,2,1,6,2,3,3,1,5,5,2,2,1,1,7,1,2,4,4,2,5,1.3,
%U 2,5,1,8,1,2,3,3,35,1,6,3,2,1,8,2,2,5,1.8,3,2,2,2,2,8,1,4,5,1,5,5,2,1,8,1,5
%N大于5的N的除数。
%H Seiichi Manyama,n的表格,n的a(n)=1..10000</a>
%H<a href=“/index/Di#divisors”>与数字除数相关的序列的索引条目</a>。
%F G.F.:和{k>=1}x^(6*k)/(1-x^k)。
%F L.g.F.:-log(产品{k>=6}(1-x^k)^(1/k))。
%F G.F.:和{k>=6}x^k/(1-x^k)_Seiichi Manyama,2023年1月7日
%F求和{k=1..n}a(k)~n*(log(n)+2*gamma-197/60),其中gamma是欧拉常数(A001620)_Amiram Eldar,2024年1月8日
%t表[DivisorSum[n,1&,#>5&],{n,1,110}]
%t nmax=110;系数列表[级数[Sum[x^(6k)/(1-x^k),{k,1,nmax}],{x,0,nmax{],x]//下降[#,1]&
%t nmax=110;系数列表[系列[-Log[产品[(1-x^k)^(1/k),{k,6,nmax}]],{x,0,nmax}],x]范围[0,nmax]//Drop[#,1]&
%o(PARI)a(n)=汇总(n,d,d>5);\\_米歇尔·马库斯,2021年4月22日
%o(PARI)我的(N=100,x='x+o('x^N));concat([0,0,0,0,0],Vec(总和(k=6,N,x^k/(1-x^k)))\\_Seiichi Manyama_,2023年1月7日
%A135539的Y列k=6。
%Y参考A000005、A001620、A023645、A032741、A321014、A338648、A338650、A338651、A33865、A33863。
%K nonn,简单
%O 1,12号
%A _Ilya Gutkovskiy_,2021年4月22日
%E a(1)-a(5)由_David a.Corneth于2022年6月13日编制
|