%I #36 Sep 08 2022 08:45:10
%S 0,1,3,15,63,279,1215,5319,23247,101655,444447,1943271,8496495,
%T 37149111,162426303,710173575,3105078543,13576277079,59359302495,
%U 259535569959,1134762524847,4961500994295,21693078131967,94848240361671
%N Expansion of x/(1 - 3*x - 6*x^2).
%C Binomial transform of A015443. A row of array A083857.
%C Pisano period lengths: 1, 1, 1, 1, 12, 1, 8, 1, 1, 12, 110, 1, 168, 8, 12, 2, 16, 1, 360, 12, ... - _R. J. Mathar_, Aug 10 2012
%H G. C. Greubel, <a href="/A083858/b083858.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (3,6).
%F a(n) = 3*a(n-1) + 6*a(n-2), a(0)=0, a(1)=1.
%F a(n) = (3*sqrt(33)/2 + 21/2)^(n/2)/sqrt(33) - (21/2 - 3*sqrt(33)/2)^(n/2)*(-1)^n/sqrt(33).
%F G.f.: Q(0)*x/2, where Q(k) = 1 + 1/(1 - x*(6*k+3 + 6*x )/( x*(6*k+6 + 6*x ) + 1/Q(k+1) )); (continued fraction). - _Sergei N. Gladkovskii_, Sep 21 2013
%F a(n) = B(n, k + 2^(n-1)) - B(n,k) where B(n,k) is formed by the family of recursions b(n) = 3*(b(n-1) + b(n-2))/2, with b(0) = 1 and b(1) = k, as explained further in A249861. - _Richard R. Forberg_, Nov 04 2014
%t a[n_]:=(MatrixPower[{{1,2},{1,-4}},n].{{1},{1}})[[2,1]]; Table[Abs[a[n]],{n,-1,40}] (* _Vladimir Joseph Stephan Orlovsky_, Feb 19 2010 *)
%t LinearRecurrence[{3,6}, {0,1}, 30] (* _G. C. Greubel_, Jan 16 2018 *)
%o (Sage) [lucas_number1(n,3,-6) for n in range(0, 24)] # _Zerinvary Lajos_, Apr 22 2009
%o (PARI) x='x+O('x^30); concat([0], Vec(x/(1-3*x-6*x^2))) \\ _G. C. Greubel_, Jan 16 2018
%o (Magma) I:=[0,1]; [n le 2 select I[n] else 3*Self(n-1) + 6*Self(n-2): n in [1..30]]; // _G. C. Greubel_, Jan 16 2018
%Y Cf. A015523, A015524.
%K easy,nonn
%O 0,3
%A _Paul Barry_, May 06 2003