%I #34 Jul 14 2015 16:51:32
%S 4,9,76,161,1364,2889,24476,51841,439204,930249,7881196,16692641,
%T 141422324,299537289,2537720636,5374978561,45537549124,96450076809,
%U 817138163596,1730726404001,14662949395604
%N Numerators of continued fraction convergents to sqrt(20).
%H Vincenzo Librandi, <a href="/A041030/b041030.txt">Table of n, a(n) for n = 0..200</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (0,18,0,-1).
%F a(2n-1) = ceiling(1/(4/(Fibonacci(6n)*sqrt(5)-Lucas(6n)+2)-2)), a(2n) = ceiling(1/(1-2/(Fibonacci(6n+3)*sqrt(5)-Lucas(6n+3)+2))-2). - _Thomas Baruchel_
%F G.f.: (4+9*x+4*x^2-x^3)/(1-18*x^2+x^4).
%F From _Gerry Martens_, Jul 11 2015: (Start)
%F Interspersion of 2 sequences [a0(n),a1(n)] for n>0 :
%F a0(n) = -((2+sqrt(5))/(9+4*sqrt(5))^n)+(-2+sqrt(5))*(9+4*sqrt(5))^n.
%F a1(n) = (1/(9+4*sqrt(5))^n+(9+4*sqrt(5))^n)/2. (End)
%t Table[Numerator[FromContinuedFraction[ContinuedFraction[Sqrt[20],n]]],{n,1,50}] (* _Vladimir Joseph Stephan Orlovsky_, Mar 17 2011 *)
%t Numerator[Convergents[Sqrt[20], 30]] (* _Vincenzo Librandi_, Oct 28 2013 *)
%t a0[n_] := -((2+Sqrt[5])/(9+4*Sqrt[5])^n)+(-2+Sqrt[5])*(9+4*Sqrt[5])^n //Simplify
%t a1[n_] := (1/(9+4*Sqrt[5])^n+(9+4*Sqrt[5])^n)/2 // Simplify
%t Flatten[MapIndexed[{a0[#], a1[#]} &,Range[20]]] (* _Gerry Martens_, Jul 11 2015 *)
%Y Cf. A010476, A041031, A087953.
%K nonn,cofr,frac,easy
%O 0,1
%A _N. J. A. Sloane_