%I #38 Oct 14 2022 08:55:28
%S 1,4,24,160,1088,7424,50688,346112,2363392,16138240,110198784,
%T 752484352,5138284544,35086401536,239584935936,1635988275200,
%U 11171226714112,76281907511296,520885446377472,3556828310929408,24287542916415488
%N a(n) = 8*a(n-1) - 8*a(n-2), a(0)=1, a(1)=4.
%C Binomial transform of A001541.
%C Let A be the unit-primitive matrix (see [Jeffery]) A = A_(8,3) = [0,0,0,1; 0,0,2,0; 0,2,0,1; 2,0,2,0]. Then A084130(n) = (1/4)*Trace(A^(2*n)). (Cf. A006012, A001333.) - _L. Edson Jeffery_, Apr 04 2011
%C a(n) is also the rational part of the Q(sqrt*(2)) integer giving the length L(n) of a variant of the Lévy C-curve, given by _Kival Ngaokrajan_, at iteration step n. See A057084. - _Wolfdieter Lang_, Dec 18 2014
%H G. C. Greubel, <a href="/A084130/b084130.txt">Table of n, a(n) for n = 0..1000</a>
%H L. E. Jeffery, <a href="/wiki/User:L._Edson_Jeffery/Unit-Primitive_Matrices">Unit-primitive matrices</a>
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (8,-8).
%F a(n) = (4+sqrt(8))^n/2 + (4-sqrt(8))^n/2.
%F G.f.: (1-4*x)/(1-8*x+8*x^2).
%F E.g.f.: exp(4*x)*cosh(sqrt(8)*x).
%F a(n) = A057084(n) - 4*A057084(n-1). - _R. J. Mathar_, Nov 10 2013
%F From _G. C. Greubel_, Oct 13 2022: (Start)
%F a(2*n) = 2^(3*n-1)*A002203(2*n).
%F a(2*n+1) = 2^(3*n+2)*A000129(2*n+1). (End)
%t LinearRecurrence[{8,-8},{1,4},30] (* _Harvey P. Dale_, Sep 25 2014 *)
%o (PARI) {a(n)= if(n<0, 0, real((4+ 2*quadgen(8))^n))}
%o (Magma) [n le 2 select 4^(n-1) else 8*(Self(n-1) -Self(n-2)): n in [1..41]]; // _G. C. Greubel_, Oct 13 2022
%o (SageMath)
%o A084130=BinaryRecurrenceSequence(8,-8,1,4)
%o [A084130(n) for n in range(41)] # _G. C. Greubel_, Oct 13 2022
%Y Cf. A000129, A001333, A001541, A002203, A006012, A057084, A084130, A084131.
%K easy,nonn
%O 0,2
%A _Paul Barry_, May 16 2003
|