%I#12 2018年10月9日07:28:17
%S 0,1,2,5,9,15,27,42,65,991482143084356058391145154820802769,
%电话:36594812627881451051813506172572196127821350954411755243,
%电话:6892885735106285131357161893198944243817298060363446
%A048272和A022567的N卷积。
%C也是A015723和A000009的卷积。
%H Alois P.Heinz,n的表格,n=0..10000时的a(n)</a>
%H Mircea Merca,<a href=“http://dx.doi.org/10.1016/j.jnt.2015年8月14日“>最近对正整数除数的卷积的组合解释,《数论杂志》,第160卷,2016年3月,第60-75页,推论3.5。
%F a(n)=总和{k=1..n}A048272(k)*A022567(n-k)=总和_{k=0..n}A015723(k)*A000009(n-k)。
%F a(n)~3^(1/4)*log(2)*exp(Pi*sqrt(2*n/3))/(2^(7/4)*Pi*n^(1/4))_Vaclav Kotesovec_,2018年10月9日
%p b:=proc(n)选项记忆`如果`(n=0,1,加(加(d*[0,1][1+
%p irem(d,2)],d=数值[除数](j))*b(n-j),j=1..n)/n)
%p端:
%p g:=proc(n,i)选项记忆`如果`(i*(i+1)/2<n,0,`如果`(n=0,[1,0],
%p加((1->[l[1],l[2]+l[1]*j])(g(n-i*j,i-1)),j=0..分钟(n/i,1))
%p端:
%p a:=n->添加(b(n-j)*g(j$2)[2],j=0..n):
%p序列(a(n),n=0..60);#_Alois P.Heinz,2016年6月18日
%t表[Sum[Count[#,_?OddQ]-Count[#,_?EvenQ]&@Divisors@k SeriesCoefficient[q,q^2]^-2,{q,0,#}]&[n-k],{k,n}],{n,0,40}]
%Y参见A048272、A022567、A015723、A000009。
%K nonn公司
%0、3
%A R.J.Mathar_,2016年6月18日
|