OFFSET
0,4
REFERENCES
Murat Sahin and Elif Tan, Conditional (strong) divisibility sequences, Fib. Q., 56 (No. 1, 2018), 18-31.
LINKS
Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (-1,8,8)
FORMULA
G.f.: -x*(2*x+1)/((x+1)*(8*x^2-1)). - Alois P. Heinz, Mar 10 2018
From Colin Barker, Mar 11 2018: (Start)
a(n) = (2^(3*n/2) - 1) / 7 for n even.
a(n) = 3*2^((3*(n-1))/2+1)/7 + 1/7 for n odd.
a(n) = -a(n-1) + 8*a(n-2) + 8*a(n-3) for n>2.
(End)
MAPLE
a:= n-> (<<0|1|0>, <0|0|1>, <8|8|-1>>^n. <<0, 1, 1>>)[1, 1]:
seq(a(n), n=0..35); # Alois P. Heinz, Mar 10 2018
MATHEMATICA
Fold[Append[#1, Inner[Times, 2 Boole[OddQ@ #2] + {1, 2}, {#1[[-1]], #1[[-2]]}, Plus]] &, {0, 1}, Range[2, 30]] (* or *)
CoefficientList[Series[-x (2 x + 1)/((x + 1) (8 x^2 - 1)), {x, 0, 30}], x] (* Michael De Vlieger, Mar 10 2018 *)
PROG
(PARI) concat(0, Vec(x*(1 + 2*x) / ((1 + x)*(1 - 8*x^2)) + O(x^40))) \\ Colin Barker, Mar 11 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Mar 10 2018
EXTENSIONS
More terms from Altug Alkan, Mar 10 2018
STATUS
approved