%I #26 Sep 08 2022 08:44:53
%S 1,1,3,4,59,63,185,248,3657,3905,11467,15372,226675,242047,710769,
%T 952816,14050193,15003009,44056211,59059220,870885291,929944511,
%U 2730774313,3660718824,53980837849,57641556673,169263951195,226905507868,3345941061347,3572846569215
%N Denominators of continued fraction convergents to sqrt(60).
%C Interspersion of 4 linear recurrences with constant coefficients. - _Gerry Martens_, Jun 10 2015
%H Vincenzo Librandi, <a href="/A041105/b041105.txt">Table of n, a(n) for n = 0..200</a>
%H <a href="/index/Rec#order_08">Index entries for linear recurrences with constant coefficients</a>, signature (0,0,0,62,0,0,0,-1).
%F G.f.: -(x^2-x-1)*(x^4+4*x^2+1) / ((x^4-8*x^2+1)*(x^4+8*x^2+1)). - _Colin Barker_, Nov 12 2013
%F a(n) = 62*a(n-4) - a(n-8). - _Vincenzo Librandi_, Dec 11 2013
%p numtheory:-cfrac(sqrt(60),100,'con','den'):
%p den[1..-2]; # _Robert Israel_, Jun 09 2015
%t Denominator[Convergents[Sqrt[60], 30]] (* _Vincenzo Librandi_, Dec 11 2013 *)
%t d0 := LinearRecurrence[{62, -1}, {1, 59}, 20]
%t d1 := LinearRecurrence[{62, -1}, {1, 63}, 20] (* A258684 *)
%t d2 := LinearRecurrence[{62, -1}, {3, 185}, 20]
%t d3 := LinearRecurrence[{62, -1}, {4, 248}, 20]
%t Flatten[MapIndexed[{d0[[#]] , d1[[#]], d2[[#]] , d3[[#]]} &,
%t Range[10]]] (* _Gerry Martens_, Jun 09 2015 *)
%t LinearRecurrence[{0, 0, 0, 62, 0, 0, 0, -1},{1, 1, 3, 4, 59, 63, 185, 248},30] (* _Ray Chandler_, Aug 03 2015 *)
%o (Magma) I:=[1, 1, 3, 4, 59, 63, 185, 248]; [n le 8 select I[n] else 62*Self(n-4)-Self(n-8): n in [1..40]]; // _Vincenzo Librandi_, Dec 11 2013
%Y Cf. A041104, A040052, A020817, A010513.
%Y Cf. A258684.
%K nonn,cofr,easy,frac
%O 0,3
%A _N. J. A. Sloane_
%E More terms from _Colin Barker_, Nov 12 2013
|