%I#6 2022年9月14日20:40:30

%S 0,0,1,3,4,0,0,1,3,4,16,8,9,11,12,14,0,13,46,8,11,9,12,14,16,17,19，

%电话21,22,24,26,30,33,36,8,9,11,12,14,16,17,19,21,222,24,25,27,29,30,32，

%U 33,36,40,42,46,50,52,56,58,62,66,68,72,76,22,24,25,27,29,30,32,33,35,37,38,40,42，43,45

%N a（N）=s（1）*t（N）+s（2）*ts（k）*t（n-k+1），其中k=楼层（n/2），s=A000201（下部Wythoff层序），t=A023533。

%H G.C.Greubel，n表，n=1..5000的a（n）</a>

%tb[j_]：=b[j]=和[KroneckerDelta[j，二项式[m+2,3]]，{m，0,15}]；

%t A025120[n_]：=A025120[n]=总和[下限[（n-j+2）*GoldenRatio]*b[j]，{j，下限[（n+4）/2]，n+1}]；

%t表[A025120[n]，{n，100}]（*_G.C.格鲁贝尔，2022年9月14日*）

%o（岩浆）

%o A023533:=func<n|二项式（Floor（（6*n-1）^（1/3））+2，3）ne n select 0 else 1>；

%o A025120:=func<n|（&+[楼层（k*（1+Sqrt（5））/2）*A023533（n+2-k）：[1..楼层（（n+1）/2）]]）中的k）>；

%o[A025120（n）：[1..100]]中的n；//_G.C.Greubel，2022年9月14日

%o（SageMath）

%o@CachedFunction

%o定义b（j）：（0..13）中m的返回和（布尔（j==二项式（m+2.3））

%o@CachedFunction

%o定义A025120（n）：返回总和（下限（（n-j+2）*golden_ratio）*b（j）for j in（（（n+4）//2）。。n+1））

%o[A025120（n）代表n in（1..100）]#_G.C.格鲁贝尔，2022年9月14日

%Y参考A000201、A023533。

%K nonn公司

%O 1,4型

%百灵鸟金伯利_