%I #16 Sep 08 2022 08:44:54
%S 1,1,5,31,129,160,3969,4129,20485,127039,528641,655680,16264961,
%T 16920641,83947525,520605791,2166370689,2686976480,66653806209,
%U 69340782689,344016936965,2133442404479,8877786554881,11011228959360,273147281579521,284158510538881
%N Denominators of continued fraction convergents to sqrt(164).
%H Vincenzo Librandi, <a href="/A041303/b041303.txt">Table of n, a(n) for n = 0..200</a>
%H <a href="/index/Rec">Index entries for linear recurrences with constant coefficients</a>, signature (0,0,0,0,0,4098,0,0,0,0,0,-1).
%F G.f.: -(x^10-x^9+5*x^8-31*x^7+129*x^6-160*x^5-129*x^4-31*x^3-5*x^2-x-1) / ((x^6-64*x^3-1)*(x^6+64*x^3-1)). - _Colin Barker_, Nov 15 2013
%F a(n) = 4098*a(n-6) - a(n-12). - _Vincenzo Librandi_, Dec 15 2013
%t Denominator[Convergents[Sqrt[164], 30]] (* _Vincenzo Librandi_, Dec 15 2013 *)
%t LinearRecurrence[{0,0,0,0,0,4098,0,0,0,0,0,-1},{1,1,5,31,129,160,3969,4129,20485,127039,528641,655680},30] (* _Harvey P. Dale_, Nov 25 2015 *)
%o (Magma) I:=[1,1,5,31,129,160,3969,4129,20485,127039, 528641,655680]; [n le 12 select I[n] else 4098*Self(n-6)-Self(n-12): n in [1..40]]; // _Vincenzo Librandi_, Dec 15 2013
%Y Cf. A041302, A040151.
%K nonn,frac,easy
%O 0,3
%A _N. J. A. Sloane_.
%E More terms from _Colin Barker_, Nov 15 2013