|
%I #16 Feb 05 2017 09:34:17
%S 1,59,2416,6230,47680,120824,798976,2017760,12928000,32622464,
%T 207425536,523312640,3321118720,8378415104,53147140096,134076293120,
%U 850391203840,2145307295744,13606407110656,34325263155200,217703105167360,549205596176384,3483252048265216
%N a(n) = 20*a(n-2) - 64*a(n-4).
%H Colin Barker, <a href="/A119886/b119886.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (0,20,0,-64).
%F G.f.: -(576*x^4-5050*x^3-2396*x^2-59*x-1) / ((2*x-1)*(2*x+1)*(4*x-1)*(4*x+1)). - _Colin Barker_, Nov 17 2012
%F a(n) = 2^(n-4)*(-3266 + 585*(-2)^n + 258*(-1)^n + 2583*2^n) for n>0. - _Colin Barker_, Feb 05 2017
%t M = {{0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0}, {1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0}, {0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0}, {1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0}, {0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0}, {0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0}, {0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0}, {1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1}, {1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1}, {0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0}, {0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0}, {0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0}, {0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1}, {0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0}, {0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1}, {0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0}} v[1] = Table[Fibonacci[n], {n, 1, 16}] v[n_] := v[n] = M.v[n - 1] a = Table[Floor[v[n][[1]]], {n, 1, 50}]
%t (* Second program: *)
%t A = SparseArray[{{1, 8} -> 1, Band[{1, 4}] -> 1, Band[{1, 2}, {3, 4}] -> 1, Band[{5, 6}, {7, 8}] -> 1}, {8, 8}]; M = ArrayFlatten[{{A+Transpose[A], IdentityMatrix[8]}, {IdentityMatrix[8], A+Transpose[A]}}]; v[1] = Array[ Fibonacci, 16]; v[n_] := v[n] = M.v[n-1]; A119886 = Array[v, 50][[All, 1]] (* _Jean-François Alcover_, Feb 05 2017 *)
%o (PARI) Vec(-(576*x^4-5050*x^3-2396*x^2-59*x-1) / ((2*x-1)*(2*x+1)*(4*x-1)*(4*x+1)) + O(x^30)) \\ _Colin Barker_, Feb 05 2017
%Y Cf. A060595, A060638, A002409, A099193, A000937, A099195.
%K nonn,easy
%O 0,2
%A _Roger L. Bagula_, Aug 09 2006
%E New name from _Joerg Arndt_, Feb 05 2017
|
