cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A262957 Numerators of the n-th iteration of the alternating continued fraction formed from the positive integers, starting with (1 - ...).

This page as a plain text file.
%I A262957 #80 Mar 04 2025 07:36:56
%S A262957 2,3,19,64,538,2833,29169,210308,2572158,23595915,334778571,
%T A262957 3732092084,60305234822,791741083537,14359827157009,217037153818264,
%U A262957 4366918714540522,74685204276602819,1651116684587556019,31524723785455714840,759659139498065625218,16017463672140861567617
%N A262957 Numerators of the n-th iteration of the alternating continued fraction formed from the positive integers, starting with (1 - ...).
%C A262957 As n->inf, a(n)/A263295(n) converges to 0.57663338973... (A346590); this number has a surprisingly elegant standard continued fraction representation of [0; 1, 1, 2, 1, 3, 4, 1, 5, 6, 1, 7, 8, ...].
%C A262957 From _Robert Israel_, Dec 22 2015: (Start)
%C A262957 a(n) is the numerator of b(n)/c(n) where
%C A262957 b(1) = 2, b(2) = 3, c(1) = 3, c(2) = 5,
%C A262957 b(n+1) = (((-1)^n*(n-1)+n*(n+2))*b(n) - (1+(-1)^n*(n+1))*b(n-1))/(n-(-1)^n),
%C A262957 c(n+1) = (((-1)^n*(n-1)+n*(n+2))*c(n) - (1+(-1)^n*(n+1))*c(n-1))/(n-(-1)^n).
%C A262957 Conjecture: b(n) and c(n) are coprime for all n, so that a(n) = b(n).
%C A262957 I have verified this for n <= 10000. (End)
%H A262957 Robert Israel, <a href="/A262957/b262957.txt">Table of n, a(n) for n = 1..448</a>
%H A262957 Peter Bala, <a href="/A262957/a262957.pdf">A note on A262957 and A263295</a>
%e A262957 (1-1/(2+1)) = 2/3, so a(1) = 2;
%e A262957 (1-1/(2+1/(3-1))) = 3/5, so a(2) = 3;
%e A262957 (1-1/(2+1/(3-1/(4+1)))) = 19/33, so a(3) = 19;
%e A262957 (1-1/(2+1/(3-1/(4+1/(5-1))))) = 64/111, so a(4) = 64.
%p A262957 P[1]:= 2: P[2]:= 3:
%p A262957 Q[1]:= 3; Q[2]:= 5;
%p A262957 for i from 2 to 100 do
%p A262957   P[i+1]:= ((-1)^i*(i-1) + i^2 + 2*i)/(i-(-1)^i)*P[i] + (1 + (i+1)*(-1)^i)/((-1)^i-i)*P[i-1];
%p A262957   Q[i+1]:= ((-1)^i*(i-1) + i^2 + 2*i)/(i-(-1)^i)*Q[i] + (1 + (i+1)*(-1)^i)/((-1)^i-i)*Q[i-1];
%p A262957 od:
%p A262957 seq(numer(P[i]/Q[i]),i=1..100); # _Robert Israel_, Dec 22 2015
%o A262957 (PARI) a(n) = if(n%2==0, s=-1, s=1); t=1; while(n>-1, t=n+1+s/t; n--; s=-s); denominator(t=1/t)
%o A262957 vector(30, n, a(n)) \\ _Mohamed Sabba_, Dec 22 2015
%Y A262957 Same principle as A244279 and A244280 - except here we begin with subtraction rather than addition.
%Y A262957 Cf. A263295 (denominators), A346590.
%K A262957 nonn,frac
%O A262957 1,1
%A A262957 _Mohamed Sabba_, Nov 19 2015
%E A262957 More terms from _Mohamed Sabba_, Dec 22 2015