|
%I #12 Jun 13 2015 00:53:11
%S 1,4,16,64,320,1344,5440,21824,103744,447808,1840448,7427392,33985856,
%T 148624704,619779392,2521191744,11221570880,49269495104,207933019456,
%U 853358105920,3726080251200,16339070997824,69569923978560
%N G.f.: (1+3*x+12*x^2+48*x^3)/(1-x-256*x^4)
%H <a href="/index/Rec">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,0,256).
%F a(n)= +a(n-1) + 256*a(n-4).
%t Clear[M, v, t, n];
%t M = {{0, t, 0, 0}, {0, 0, t, 0}, {0, 0, 0, t}, {t, 0, 0, 1/t}};
%t v[0] = {1, 1, 1, 1};
%t v[n_] := v[n] = M.v[n - 1];
%t CharacteristicPolynomial[M, x];
%t t = 2;
%t a = Table[t^n*v[n][[1]], {n, 0, 30}]
%o (PARI) Vec((1+3*x+12*x^2+48*x^3)/(1-x-256*x^4)+O(x^44))
%K nonn,easy
%O 0,2
%A _Roger L. Bagula_, Mar 26 2009
|
