登录
对于n<4,在Euler变换下左移4位,a(0)=0,a(n)=1。

%I#19 2019年5月8日14:34:31

%S 0,1,1,1,1,2,3,5,7,12,18,30,47,7812520934157194615922663,

%电话4503759412898218913733463691109039186816320913551829950842,

%电话:16401492833664901658849001914720477255535254440163877232183134457819

%N在Euler变换下左移4位,a(0)=0,a(N)=1表示N<4。

%H Alois P.Heinz,n的表格,n=0..1000时的a(n)</a>

%H N.J.A.Sloane,转换</a>

%F a(n)~c*d^n/n^(3/2),其中d=1.8065918193702780027972…和c=1.041173202249532389463_瓦茨拉夫·科特索维奇,2014年6月23日

%F G.F.:x+x^2+x^3+x^4/产品{n>=1}(1-x^n)^a(n).-_伊利亚·古特科夫斯基,2019年5月8日

%p(数字理论):

%p b:=proc(n)选项记忆`如果`(n=0,1,

%p(加(加(d*a(d),d=除数(j))*b(n-j),j=1..n))/n)

%p端:

%p a:=n->`如果`(n<4,符号(n),b(n-4)):

%p序列(a(n),n=0..45);

%tb[n_]:=b[n]=如果[n==0,1,(和[Sum[d*a[d],{d,除数[j]}]*b[n-j],{j,1,n}])/n];a[n_]:=如果[n<4,符号[n],b[n-4];表[a[n],{n,0,41}](*_Jean-François Alcover_,2013年8月1日,摘自_Alois P.Heinz_*)

%A144018的Y列k=4。

%Y参考A316076。

%K nonn,本征

%0、7

%A _Alois P.Heinz,2012年10月18日