login
a(n+3) = 5*a(n+2) - 2*a(n+1) - 8*a(n), a(0) = 1, a(1) = 0, a(2) = 3.
3

%I #8 Sep 01 2017 03:05:23

%S 1,0,3,7,29,107,421,1659,6597,26299,105029,419771,1678405,6712251,

%T 26846277,107379643,429507653,1718008763,6871991365,27487878075,

%U 109951337541,439805000635,1759219303493,7036875815867,28147500467269

%N a(n+3) = 5*a(n+2) - 2*a(n+1) - 8*a(n), a(0) = 1, a(1) = 0, a(2) = 3.

%C A Jacobsthal related sequence (A001045). This sequence was calculated using the same rules given for A108618; the "initial seed" is the floretion given in the program code, below.

%H G. C. Greubel, <a href="/A110613/b110613.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 (5,-2,-8).

%F G.f.: (1-5*x+5*x^2)/((4*x-1)*(2*x-1)*(x+1)).

%F Program "Superseeker" finds:

%F a(n) + a(n+1) = A007582(n) = A007581(n+1) - A007581(n).

%F a(n+2) - a(n) = A049775(n).

%F a(n) + 2*a(n+1) + a(n+2) = A087440(n+1).

%p seriestolist(series((1-5*x+5*x^2)/((4*x-1)*(2*x-1)*(x+1)), x=0,25)); -or- Floretion Algebra Multiplication Program, FAMP Code: 2tessumseq[(.5'i - .5'k - .5i' + .5k' - .5'ij' - .5'ji' - .5'jk' - .5'kj')('i + j' + 'ij' + 'ji')] Sumtype is set to:sum[Y[15]] = sum(*) (from 3rd term, disregarding signs)

%t LinearRecurrence[{5, -2, -8}, {1, 0, 3}, 50] (* _G. C. Greubel_, Sep 01 2017 *)

%o (PARI) x='x+O('x^50); Vec((1-5*x+5*x^2)/((4*x-1)*(2*x-1)*(x+1))) \\ _G. C. Greubel_, Sep 01 2017

%Y Cf. A001045, A007582, A007581, A049775, A087440, A108618, A110614.

%K easy,nonn

%O 0,3

%A _Creighton Dement_, Jul 31 2005