%I#26 2022年9月8日08:44:59
%S 1,0,1,3,2,7,13,18,41,711222394217621725434642849515389,
%电话280825117793047169610308847562197102417018648413395711,
%电话:618449811261551205077893734691468008809123848199225535258
%N(1-x)/(1-x-x^2-2*x^3+2*x^4)的展开。
%H G.C.Greubel,n表,n=0..1000时的a(n)</a>
%H INRIA算法项目,<a href=“http://ecs.inria.fr/services/structure?nbr=482“>组合结构百科全书482</a>
%H<a href=“/index/Rec#order_04”>具有常系数的线性重复出现的索引条目,签名(1,1,2,-2)。
%财务报表:(1-x)/(1-x-2*x^3+2*x^4-x^2)。
%F a(n)=a(n-1)+a(n-2)+2*a(n-3)-2*a(n-4),其中a(0)=1,a(1)=0,a(2)=1,a(3)=3。
%F a(n)=Sum_{alpha=RootOf(1-x-2*x^3+2*x^4-x^2)}(-1/353*(-18-106*alpha+33*alpha_2+28*alpha^3)*alpha(-1-n))。
%p规范:=[S,{S=序列(Prod(Z,Z,Union(Z,Z,Sequence(Z))))},未标记]:seq(combstruct[计数](规范,大小=n),n=0..40);
%t系数列表[级数[(1-x)/(1-x-x^2-2x^3+2x^4),{x,0,40}],x](*或*)线性递归[{1,1,2,-2},{1,0,1,3},40](*H arvey P.Dale_,2017年7月2日*)
%o(PARI)我的(x='x+o('x^40));Vec((1-x)/(1-x-2*x^3+2*x^4-x^2))\\_G.C.格鲁贝尔,2019年5月8日
%o(Magma)R<x>:=PowerSeriesRing(整数(),40);系数(R!((1-x)/(1-x-2*x^3+2*x^4-x^2));//_G.C.Greubel,2019年5月8日
%o(鼠尾草)((1-x)/(1-x-2*x^3+2*x^4-x^2)).系列(x,40).系数(x,稀疏=假)#_G.C.格鲁贝尔,2019年5月8日
%o(间隙)a:=[1,0,1,3];;对于[5..40]中的n,做a[n]:=a[n-1]+a[n-2]+2*a[n-3]-2*a[n-4];od;a、 #个_G.C.Greubel_,2019年5月8日
%K容易,不是
%0、4
%百科全书(AT)pommard.inia.fr,2000年1月25日
%E来自James A.Sellers_的更多条款,2000年6月6日
|