cp's OEIS Frontend

A215466 Expansion of x*(1-4*x+x^2) / ( (x^2-7*x+1)*(x^2-3*x+1) ).

%I A215466 #57 Nov 10 2024 13:18:37
%S A215466 0,1,6,38,252,1705,11628,79547,544824,3733234,25585230,175356611,
%T A215466 1201893336,8237850373,56462937882,387002396990,2652553009008,
%U A215466 18180866487757,124613506702404,854113665498719,5854182112700460
%N A215466 Expansion of x*(1-4*x+x^2) / ( (x^2-7*x+1)*(x^2-3*x+1) ).
%C A215466 From _Peter Bala_, Aug 05 2019: (Start)
%C A215466 Let  U(n;P,Q), where P and Q are integer parameters, denote the Lucas sequence of the first kind. Then, excluding the case P = -1, the sequence ( U(n;P,1) + U(2*n;P,1) )/(P + 1) is a fourth-order linear divisibility sequence with o.g.f. x*(1 - 2*(P - 1)*x + x^2)/((1 - P*x + x^2)*(1 - (P^2 - 2)*x + x^2)). This is the case P = 3. See A000027 (P = 2), A165998 (P = -2) and A238536 (P = -3).
%C A215466 More generally, the sequence U(n;P,1) + U(2*n;P,1) + ... + U(k*n;P,1) is a linear divisibility sequence of order 2*k. As an example, see A273625 (P = 3, k = 3 and then sequence normalized with initial term 1). (End)
%H A215466 Vincenzo Librandi, <a href="/A215466/b215466.txt">Table of n, a(n) for n = 0..1000</a>
%H A215466 Wikipedia, <a href="https://en.wikipedia.org/wiki/Lucas_sequence">Lucas sequence</a>
%H A215466 Peter Bala, <a href="/A238600/a238600_1.pdf">Divisibility sequences from strong divisibility sequences</a>
%H A215466 E. L. Roettger and H. C. Williams, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL24/Roettger/roettger12.html">Appearance of Primes in Fourth-Order Odd Divisibility Sequences</a>, J. Int. Seq., Vol. 24 (2021), Article 21.7.5.
%H A215466 H. C. Williams and R. K. Guy, <a href="http://dx.doi.org/10.1142/S1793042111004587">Some fourth-order linear divisibility sequences</a>, Intl. J. Number Theory 7 (5) (2011) 1255-1277.
%H A215466 H. C. Williams and R. K. Guy, <a href="http://www.emis.de/journals/INTEGERS/papers/p33/p33.Abstract.html">Odd and even linear divisibility sequences of order 4</a>, INTEGERS, 2015, #A33.
%H A215466 <a href="/index/Di#divseq">Index to divisibility sequences</a>
%H A215466 <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (10,-23,10,-1).
%F A215466 a(n) = L(n)*F(3n)/4 if n even, = F(n)*L(3n)/4 if n odd, where L=A000032, F=A000045.
%F A215466 a(n) = 3*A004187(n)/4 + A001906(n)/4.
%F A215466 a(n) = 10*a(n-1) - 23*a(n-2) + 10*a(n-3) - a(n-4), a(0)=0, a(1)=1, a(2)=6, a(3)=38. - _Harvey P. Dale_, Nov 02 2015
%F A215466 a(n) = (1/4)*(Fibonacci(2*n) + Fibonacci(4*n)) = (1/4)*(A001906(n) + A033888(n)). - _Peter Bala_, Aug 05 2019
%F A215466 E.g.f.: exp(5*x/2)*(cosh(x)+exp(x)*cosh(sqrt(5)*x))*sinh(sqrt(5)*x/2)/sqrt(5). - _Stefano Spezia_, Aug 17 2019
%F A215466 a(n) = -a(-n) for all n in Z. - _Michael Somos_, Dec 29 2022
%p A215466 A215466 := proc(n)
%p A215466     if type(n,'even') then
%p A215466         A000032(n)*combinat[fibonacci](3*n)/4 ;
%p A215466     else
%p A215466         combinat[fibonacci](n)*A000032(3*n)/4 ;
%p A215466     end if;
%p A215466 end proc:
%t A215466 CoefficientList[Series[x*(1 - 4*x + x^2)/((x^2 - 7*x + 1)*(x^2 - 3*x + 1)), {x, 0, 40}], x] (* _Vincenzo Librandi_, Dec 23 2012 *)
%t A215466 LinearRecurrence[{10,-23,10,-1},{0,1,6,38},30] (* _Harvey P. Dale_, Nov 02 2015 *)
%o A215466 (Magma) I:=[0,1,6,38]; [n le 4 select I[n] else 10*Self(n-1)-23*Self(n-2)+10*Self(n-3)-Self(n-4): n in [1..30]]; // _Vincenzo Librandi_, Dec 23 2012
%o A215466 (Magma) /* By definition: */ m:=20; R<x>:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!((1-4*x+x^2)/((x^2-7*x+1)*(x^2-3*x+1)))); // _Bruno Berselli_, Dec 24 2012
%o A215466 (PARI) a(n)=([0,1,0,0; 0,0,1,0; 0,0,0,1; -1,10,-23,10]^n*[0;1;6;38])[1,1] \\ _Charles R Greathouse IV_, Nov 13 2015
%o A215466 (PARI) {a(n) = my(w=quadgen(5)^(2*n)); imag(w^2+w)/4}; /* _Michael Somos_, Dec 29 2022 */
%Y A215466 Cf. A000032, A000045, A001906, A033888, A165998, A215465, A238536, A273625.
%K A215466 nonn,easy
%O A215466 0,3
%A A215466 _R. J. Mathar_, Aug 11 2012