%I #20 Sep 08 2022 08:44:54
%S 1,1,2,3,5,8,13,528,541,1069,1610,2679,4289,6968,283009,289977,572986,
%T 862963,1435949,2298912,3734861,151693352,155428213,307121565,
%U 462549778,769671343,1232221121,2001892464,81307919681,83309812145,164617731826,247927543971
%N Denominators of continued fraction convergents to sqrt(425).
%H Vincenzo Librandi, <a href="/A041809/b041809.txt">Table of n, a(n) for n = 0..200</a>
%H <a href="/index/Rec#order_14">Index entries for linear recurrences with constant coefficients</a>, signature (0, 0, 0, 0, 0, 0, 536, 0, 0, 0, 0, 0, 0, 1).
%F G.f.: -(x^12-x^11+2*x^10-3*x^9+5*x^8-8*x^7+13*x^6+8*x^5+5*x^4+3*x^3+2*x^2+x+1) / (x^14+536*x^7-1). - _Colin Barker_, Nov 25 2013
%F a(n) = 536*a(n-7) + a(n-14) for n>13. - _Vincenzo Librandi_, Dec 24 2013
%t Denominator[Convergents[Sqrt[425], 30]] (* _Vincenzo Librandi_, Dec 24 2013 *)
%t LinearRecurrence[{0,0,0,0,0,0,536,0,0,0,0,0,0,1},{1,1,2,3,5,8,13,528,541,1069,1610,2679,4289,6968},40] (* _Harvey P. Dale_, Dec 02 2014 *)
%o (Magma) I:=[1,1,2,3,5,8,13,528,541,1069,1610,2679,4289, 6968]; [n le 14 select I[n] else 536*Self(n-7)+Self(n-14): n in [1..50]]; // _Vincenzo Librandi_, Dec 24 2013
%Y Cf. A041808, A040404.
%K nonn,frac,easy
%O 0,3
%A _N. J. A. Sloane_.
%E More terms from _Colin Barker_, Nov 25 2013
|