|
%I
%S 25,2,3,2,2,1,1,4,25,4,1,1,2,2,3,2,50,2,3,2,2,1,1,4,25,4,1,1,2,2,3,2,
%T 50,2,3,2,2,1,1,4,25,4,1,1,2,2,3,2,50,2,3,2,2,1,1,4,25,4,1,1,2,2,3,2,
%U 50,2,3,2,2,1,1,4,25,4,1,1,2,2,3,2,50,2,3,2,2,1,1,4
%N Continued fraction for sqrt(647).
%H <a href="/index/Con#confC">Index entries for continued fractions for constants</a>
%H <a href="/index/Rec#order_16">Index entries for linear recurrences with constant coefficients</a>, signature (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1).
%F a(n)=(1/128)*{-377*(n mod 16)+15*[(n+1) mod 16]-[(n+2) mod 16]+7*[(n+3) mod 16]-[(n+4) mod 16]+7*[(n+5) mod 16]+31*[(n+6) mod 16]+175*[(n+7) mod 16]-161*[(n+8) mod 16]-17*[(n+9) mod 16]+7*[(n+10) mod 16]+15*[(n+11) mod 16]+7*[(n+12) mod 16]+15*[(n+13) mod 16]-[(n+14) mod 16]+391*[(n+15) mod 16]}-25*[C(2*n,n) mod 2], with n>=0 [From _Paolo P. Lava_, May 15 2009]
%p with(numtheory): Digits := 300: convert(evalf(sqrt(647)),confrac);
%K nonn,cofr,easy
%O 0,1
%A _N. J. A. Sloane_.
|
