|
%I #24 Sep 08 2022 08:45:06
%S 3,4,8,22,64,184,524,1488,4224,11992,34048,96672,274480,779328,
%T 2212736,6282592,17838080,50647424,143802560,408296704,1159271424,
%U 3291504000,9345523712,26534621696,75339399936,213910160384,607352285184
%N Binomial transform of generalized Lucas numbers S(n) = S(n-1) + S(n-2) + S(n-3), S(0)=3, S(1)=1, S(2)=3.
%C Limit as n-> infinity of a(n)/a(n-1) is 1+c, where c = 1.83928675...
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, San Diego, 1995.
%H G. C. Greubel, <a href="/A073313/b073313.txt">Table of n, a(n) for n = 0..1000</a>
%H M. Bernstein and N. J. A. Sloane, <a href="http://arXiv.org/abs/math.CO/0205301">Some canonical sequences of integers</a>, Linear Alg. Applications, 226-228 (1995), 57-72; erratum 320 (2000), 210. [Link to arXiv version]
%H M. Bernstein and N. J. A. Sloane, <a href="/A003633/a003633_1.pdf">Some canonical sequences of integers</a>, Linear Alg. Applications, 226-228 (1995), 57-72; erratum 320 (2000), 210. [Link to Lin. Alg. Applic. version together with omitted figures]
%H H. Prodinger, <a href="https://www.fq.math.ca/Scanned/32-5/prodinger.pdf">Some information about the binomial transform.</a>, The Fibonacci Quarterly, 32, 1994, 412-415.
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (4,-4,2).
%F a(n) is the trace of the n-th power of 3 X 3 matrix: first row (2, 1, 0), second row (1, 1, 1), third row (1, 0, 1). It satisfies recurrence a(n) = 4*a(n-1) - 4*a(n-2) + 2*a(n-3), a(0)=3, a(1)=4, a(2)=8.
%F G.f.: (3 - 8*x + 4*x^2)/(1 - 4*x + 4*x^2 - 2*x^3).
%t f[n_]:= f[n]=4*f[n-1]-4*f[n-2]+2*f[n-3]; f[0]=3; f[1]=4; f[2]=8; Table[f[n], {n, 0, 30}]
%t LinearRecurrence[{4,-4,2},{3,4,8},30] (* _Harvey P. Dale_, May 08 2015 *)
%o (PARI) my(x='x+O('x^30)); Vec((3-8*x+4*x^2)/(1-4*x+4*x^2-2*x^3)) \\ _G. C. Greubel_, Apr 22 2019
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 30); Coefficients(R!( (3-8*x+4*x^2)/(1-4*x+4*x^2-2*x^3) )); // _G. C. Greubel_, Apr 22 2019
%o (Sage) ((3-8*x+4*x^2)/(1-4*x+4*x^2-2*x^3)).series(x, 30).coefficients(x, sparse=False) # _G. C. Greubel_, Apr 22 2019
%Y Cf. A001644.
%K easy,nonn
%O 0,1
%A Mario Catalani (mario.catalani(AT)unito.it), Jul 26 2002
|
