%I #12 Sep 08 2022 08:45:43
%S 2,4,9,23,55,126,291,678,1578,3667,8523,19815,46066,107089,248950,
%T 578740,1345409,3127695,7271007,16903042,39294807,91349342,212361454,
%U 493680487,1147667895,2668004163,6202357186,14418731129,33519483178
%N Transform of A056594 by the T_{1,1} transformation (see link)
%H G. C. Greubel, <a href="/A159334/b159334.txt">Table of n, a(n) for n = 0..1000</a>
%H R. Choulet, <a href="http://www.apmep.asso.fr/IMG/pdf/curtz1.pdf">Curtz-like transformation </a>
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,4,-2,1).
%F O.g.f.: -(2-2*x+3*x^2+x^4)/((x^2+1)*(x^3-2*x^2+3*x-1)).
%F for n>=0 a(n+5)=3*a(n+4)-3*a(n+3)+4*a(n+2)-2*a(n+1)+a(n)
%t Join[{2, 4, 9, 23, 55, 126}, LinearRecurrence[{3, -3, 4, -2, 1}, {291, 678, 1578, 3667, 8523}, 45]] (* _G. C. Greubel_, Jun 25 2018 *)
%o (PARI) x='x+O('x^50); Vec(-(2-2*x+3*x^2+x^4)/((x^2+1)*(x^3-2*x^2+3*x-1))) \\ _G. C. Greubel_, Jun 25 2018
%o (Magma) I:=[291, 678, 1578, 3667, 8523]; [2, 4, 9, 23, 55, 126] cat [n le 5 select I[n] else 3*Self(n-1) - 3*Self(n-2) + 4*Self(n-3) -2*Self(n-4) +Self(n-5): n in [1..50]]; // _G. C. Greubel_, Jun 25 2018
%Y Cf. A159328, A159329, A159330, A159331.
%K easy,nonn
%O 0,1
%A _Richard Choulet_, Apr 10 2009