login
A098581
Expansion of (1+2*x+4*x^2)/(1-x-8*x^4).
1
1, 3, 7, 7, 15, 39, 95, 151, 271, 583, 1343, 2551, 4719, 9383, 20127, 40535, 78287, 153351, 314367, 638647, 1264943, 2491751, 5006687, 10115863, 20235407, 40169415, 80222911, 161149815, 323033071, 644388391, 1286171679, 2575370199
OFFSET
0,2
FORMULA
a(n) = a(n-1) + 8*a(n-4).
a(n) = Sum_{k=0..n} binomial(n-k, floor(k/3)) * 2^k.
MATHEMATICA
CoefficientList[Series[(1+2x+4x^2)/(1-x-8x^4), {x, 0, 40}], x] (* or *) LinearRecurrence[{1, 0, 0, 8}, {1, 3, 7, 7}, 40] (* Harvey P. Dale, Feb 04 2015 *)
PROG
(PARI) x='x+O('x^30); Vec((1+2*x+4*x^2)/(1-x-8*x^4)) \\ G. C. Greubel, Feb 03 2018
(Magma) I:=[1, 3, 7, 7]; [n le 4 select I[n] else Self(n-1) + 8*Self(n-4): n in [1..30]]; // G. C. Greubel, Feb 03 2018
CROSSREFS
Sequence in context: A081218 A240422 A130003 * A238997 A240260 A240427
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Sep 16 2004
STATUS
approved