login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A158761
a(n)=a(n-2)+64*a(n-3) with a(0)=1, a(1)=4, a(3)=16.
0
1, 4, 16, 68, 272, 1092, 4624, 18500, 74512, 314436, 1258512, 5083204, 21382416, 85627972, 346707472, 1454102596, 5826897680, 23643380804, 98889463824, 396564832324, 1612065835280, 6725490517060, 26992215104016
OFFSET
0,2
COMMENTS
The sequence is also 2^n times the sum of the three top-row entries
of the n-th power of the matrix (0,2,0; 0,0,2;, 2,1/8,0) .
FORMULA
G.f.: -(1+4*x+15*x^2)/(-1+x^2+64*x^3). [From R. J. Mathar, Mar 26 2009]
MATHEMATICA
Clear[M, v, t, n];
M = {{0, t, 0}, {0, 0, t}, {t, 1/t^3, 0}};
v[0] = {1, 1, 1};
v[n_] := v[n] = M.v[n - 1];
CharacteristicPolynomial[M, x];
t = 2; a = Table[t^n*v[n][[1]], {n, 0, 30}]
CROSSREFS
Sequence in context: A259815 A283036 A307051 * A179611 A290912 A089979
KEYWORD
nonn,easy
AUTHOR
Roger L. Bagula, Mar 25 2009
EXTENSIONS
Definition replaced by recurrence - The Associate Editors of the OEIS - Oct 05 2009
STATUS
approved