%I M3414 N1380#33 2021年12月20日20:18:42
%S 4,11,26,52,9817128946773711311704251536615246743010396,
%电话:14405197602684362694858364614853991121701462261930362,
%电话:24609931621405556517224657012831320104805513166111648486205732425597193175309
%N二部分分区的数量。
%D M.S.Cheema和H.Gupta,高斯整数分区表。印度国家科学院,《数学表》,第1卷,新德里,1956年,第11页。
%D N.J.A.Sloane,《整数序列手册》,学术出版社,1973年(包括该序列)。
%D N.J.A.Sloane和Simon Plouffe,《整数序列百科全书》,学术出版社,1995年(包括该序列)。
%H Alois P.Heinz,n的表格,n=0..300时的a(n)</a>
%H M.S.Cheema和H.Gupta,《高斯整数分区表》。印度国家科学院,数学表,第1卷,新德里,1956年
%F a(n)=a(n-1)+A000041(n)+A00.0070(n)+A000291(n),对于n>0-阿诺德,2007年12月10日
%F From _Vaclav Kotesovec_,2017年1月7日:(开始)
%财务报表:(4-x-3*x^2+x^3)/((1-x)^3*(1+x))*产品{k>=1}1/(1-x^k)。
%F a(n)~exp(Pi*sqrt(2*n/3))*3*sqert(n)/(2*sqort(2)*Pi^3)。
%F(完)
%p(数字理论):
%p b:=proc(n,k)选项记忆;
%p`if`(n>k,0,1)+`if`,
%p加法(`if`(d>k,0,b(n/d,d)),d=除数(n)减去{1,n})
%p端:
%pa:=n->b((45*2^n)$2):
%p序列(a(n),n=0..50);#_Alois P.Heinz,2013年5月26日
%tb[n_,k_]:=b[n,k]=如果[n>k,0,1]+如果[PrimeQ[n],0,Sum[If[d>k,0,b[n/d,d]],{d,DeleteCase[Divisors[n]、1|n]}]];a[n]:=b[45*2^n,45*2*n];表[a[n],{n,0,50}](*_Jean-François Alcover_,2014年3月20日,在_Alois P.Heinz_*之后)
%t nmax=100;系数列表[系列[(4-x-3*x^2+x^3)/(1-x)^3*(1+x))/产品[1-x^k,{k,1,nmax}],{x,0,nmax{],x](*Vaclav Kotesovec_,2017年1月7日*)
%Y参考A082775,A129306。
%K nonn公司
%0、1
%A _N.J.A.斯隆_
%E延伸超过a(25),作者:_Alois P.Heinz,2013年5月26日
|