|
%I #32 Dec 31 2023 00:27:51
%S 1,0,4,4,20,36,116,260,724,1764,4660,11716,30356,77220,198644,507524,
%T 1302100,3332196,8540596,21869380,56031764,143509284,367636340,
%U 941673476,2412218836,6178912740,15827788084,40543439044,103854591380
%N Expansion of (1-x)/(1 - x - 4*x^2).
%C First differences of A006131.
%C This sequence {a(n)} appears in the formula for powers of c = (1 + sqrt(17))/2 = A222132, the fundamental (integer) algebraic number of Q(sqrt(17)): c^n = a(n) + A006131(n-1)*c. This is also valid for positive powers of 1/c = (-1 + sqrt(17)) /8. See the formula below and in A006131 in terms of Chebyshev or Fibonacci polynomials. - _Wolfdieter Lang_, Nov 27 2023
%H G. C. Greubel, <a href="/A052923/b052923.txt">Table of n, a(n) for n = 0..1000</a>
%H INRIA Algorithms Project, <a href="http://ecs.inria.fr/services/structure?nbr=908">Encyclopedia of Combinatorial Structures 908</a>
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (1,4).
%H <a href="/index/Ch#Cheby">Index entries for sequences related to Chebyshev polynomials.</a>
%F G.f.: (1-x)/(1 - x - 4*x^2).
%F a(n) = a(n-1) + 4*a(n-2), with a(0)=1, a(1)=0.
%F a(n) = Sum_{alpha=RootOf(-1+z+4*z^2)} (1/17)*(-1+9*alpha)*alpha^(-1-n).
%F If p[1]=0, and p[i]=4, ( i>1), and if A is Hessenberg matrix of order n defined by: A[i,j]=p[j-i+1], (i<=j), A[i,j]=-1, (i=j+1), and A[i,j]=0 otherwise. Then, for n>=1, a(n)=det A. - _Milan Janjic_, Apr 29 2010
%F From _Wolfdieter Lang_, Nov 27 2023: (Start)
%F a(n) = 4*A006131(n-2), with A006131(-2) = 1/4 and A006131(-1) = 0.
%F a(n) = -(-2*i)^n*S(n-2, i/2), with i = sqrt(-1), and the S-Chebyshev polynomials (see A049310). S(-n, x) = -S(n-2, x). The Fibonacci polynomials are F(n, x) = (-i)^(n-1)*S(n-1, i*x). (End)
%p spec := [S,{S=Sequence(Prod(Sequence(Z),Z,Union(Z,Z,Z,Z)))},unlabeled]: seq(combstruct[count](spec,size=n), n=0..20);
%p seq(coeff(series((1-x)/(1 -x -4*x^2), x, n+1), x, n), n = 0..30); # _G. C. Greubel_, Oct 16 2019
%t LinearRecurrence[{1,4}, {1,0}, 30] (* _G. C. Greubel_, Oct 16 2019 *)
%o (PARI) my(x='x+O('x^30)); Vec((1-x)/(1 -x -4*x^2)) \\ _G. C. Greubel_, Oct 16 2019
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 30); Coefficients(R!( (1-x)/(1 -x -4*x^2) )); // _G. C. Greubel_, Oct 16 2019
%o (Sage)
%o def A052923_list(prec):
%o P.<x> = PowerSeriesRing(ZZ, prec)
%o return P((1-x)/(1 -x -4*x^2)).list()
%o A052923_list(30) # _G. C. Greubel_, Oct 16 2019
%o (GAP) a:=[1,0];; for n in [3..30] do a[n]:=a[n-1]+4*a[n-2]; od; a; # _G. C. Greubel_, Oct 16 2019
%o a := n -> -(2*I)^n*ChebyshevU(n-2, -I/4):
%o seq(simplify(a(n)), n = 0..28); # _Peter Luschny_, Dec 03 2023
%Y Cf. A006131, A026581, A049310, A222132.
%K easy,nonn
%O 0,3
%A encyclopedia(AT)pommard.inria.fr, Jan 25 2000
%E More terms from _James A. Sellers_, Jun 06 2000
|
