%I#12 2022年11月30日08:25:02
%S 1,0,1,1,2,5,12,28,68181531167154911862765299237880903907,
%电话:358061914729777762639952274442521123673024455729809348,
%电话:27292248240133614280479671803467718030415533464970976743183094283634259901080027574354621878244653093329527120022
%N G.f.A(x)满足:A(x”)=1+x^2+x^3*A(x/(1-x))/(1-x)。
%H G.C.Greubel,n表,n=0..695的a(n)</a>
%F a(0)=1,a(1)=0,a(2)=1;a(n)=和{k=0..n-3}二项式(n-3,k)*a(k)。
%t nmax=31;A[_]=0;做[A[x_]=1+x^2+x^3 A[x/(1-x)]/(1-x)+O[x]^(nmax+1)//正常,nmax+1];系数列表[A[x],x]
%ta[0]=1;a[1]=0;a[2]=1;a[n]:=a[n]=和[二项式[n-3,k]a[k],{k,0,n-3}];表[a[n],{n,0,31}]
%o(岩浆)
%o函数a(n)
%o如果n lt 3,则返回(1+(-1)^n)/2;
%o否则返回(&+[二项式(n-3,j)*a(j):[0..n-3]]中的j);
%o结束条件:;返回a;
%o端函数;
%o[0..35]]中的[a(n):n;//_G.C.Greubel,2022年11月30日
%o(SageMath)
%o@CachedFunction
%o定义a(n):#a=A346051
%o如果(n<3):返回(1,0,1)[n]
%o else:返回和(二项式(n-3,k)*a(k),k在范围(n-2)内)
%o[a(n)表示n在(51)范围内]#_G.C.Greubel_,2022年11月30日
%Y参见A000994、A000995、A000996、A00099.7、A0009998。
%Y参考A007476、A210540、A346050、A346052。
%K nonn公司
%0、6
%A _Ilya Gutkovskiy,2021年7月2日
|