%I #19 Mar 11 2018 15:34:00
%S 0,1,1,7,9,55,73,439,585,3511,4681,28087,37449,224695,299593,1797559,
%T 2396745,14380471,19173961,115043767,153391689,920350135,1227133513,
%U 7362801079,9817068105,58902408631,78536544841,471219269047,628292358729,3769754152375,5026338869833
%N a(n) = a(n-1) + 2*a(n-2) if n even, or 3*a(n-1) + 4*a(n-2) if n odd, starting with 0, 1.
%D Murat Sahin and Elif Tan, Conditional (strong) divisibility sequences, Fib. Q., 56 (No. 1, 2018), 18-31.
%H Colin Barker, <a href="/A299913/b299913.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (-1,8,8)
%F G.f.: -x*(2*x+1)/((x+1)*(8*x^2-1)). - _Alois P. Heinz_, Mar 10 2018
%F From _Colin Barker_, Mar 11 2018: (Start)
%F a(n) = (2^(3*n/2) - 1) / 7 for n even.
%F a(n) = 3*2^((3*(n-1))/2+1)/7 + 1/7 for n odd.
%F a(n) = -a(n-1) + 8*a(n-2) + 8*a(n-3) for n>2.
%F (End)
%p a:= n-> (<<0|1|0>, <0|0|1>, <8|8|-1>>^n. <<0, 1, 1>>)[1,1]:
%p seq(a(n), n=0..35); # _Alois P. Heinz_, Mar 10 2018
%t Fold[Append[#1, Inner[Times, 2 Boole[OddQ@ #2] + {1, 2}, {#1[[-1]], #1[[-2]]}, Plus]] &, {0, 1}, Range[2, 30]] (* or *)
%t CoefficientList[Series[-x (2 x + 1)/((x + 1) (8 x^2 - 1)), {x, 0, 30}], x] (* _Michael De Vlieger_, Mar 10 2018 *)
%o (PARI) concat(0, Vec(x*(1 + 2*x) / ((1 + x)*(1 - 8*x^2)) + O(x^40))) \\ _Colin Barker_, Mar 11 2018
%Y Bisections give A023001, A083068.
%Y Cf. A299914, A299915, A299916.
%K nonn,easy
%O 0,4
%A _N. J. A. Sloane_, Mar 10 2018
%E More terms from _Altug Alkan_, Mar 10 2018
|