OFFSET
0,1
COMMENTS
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (7, -14, 8).
FORMULA
a(n) = 4^(n+8)+12576771*2^(n-4)-1 for n > 8.
G.f.: (753-228*x+1020*x^2+97128*x^3+517164*x^4+9133668*x^5-31930858*x^6 +497474728*x^7-28023638*x^8-3587820988*x^9+3014752848*x^10+125798400*x^11) / ((1-x)*(1-2*x)*(1-4*x)).
G.f. for the sequence starting at a(9): 3*x^9*(5860775285-17850630302*x +11989855016*x^2)/((1-x)*(1-2*x)*(1-4*x)).
MATHEMATICA
CoefficientList[Series[(753 - 228 x + 1020 x^2 + 97128 x^3 + 517164 x^4 + 9133668 x^5 - 31930858 x^6 + 497474728 x^7 - 28023638 x^8 - 3587820988 x^9 + 3014752848 x^10 + 125798400 x^11)/((1 - x) (1 - 2 x) (1 - 4 x)), {x, 0, 40}], x] (* Vincenzo Librandi, Sep 24 2013 *)
PROG
(PARI) {m=18; v=concat([753, 5043, 25779, 213003, 1687623, 18171219, 73344977, 769993485, 4480470831, 17582325855, 69524390079], vector(m-11)); for(n=12, m, v[n]=6*v[n-1]-8*v[n-2]-3); v}
(Magma) [753, 5043, 25779, 213003, 1687623, 18171219, 73344977, 769993485, 4480470831] cat [4^(n+8)+12576771*2^(n-4)-1: n in [9..25]]; // Vincenzo Librandi, Sep 24 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Klaus Brockhaus, May 07 2010
STATUS
approved