login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Consider the mapping f(a/b) = (a - b)/(ab). Taking a = 2 and b = 1 to start with and carrying out this mapping repeatedly on each new (reduced) rational number gives the following sequence 2/1,1/2,-1/2,-3/-2,-1/6,... Sequence contains the denominators.
2

%I #27 Oct 11 2021 18:45:12

%S 1,2,2,-2,6,-6,42,-42,1806,-1806,3263442,-3263442,10650056950806,

%T -10650056950806,113423713055421844361000442,

%U -113423713055421844361000442,12864938683278671740537145998360961546653259485195806

%N Consider the mapping f(a/b) = (a - b)/(ab). Taking a = 2 and b = 1 to start with and carrying out this mapping repeatedly on each new (reduced) rational number gives the following sequence 2/1,1/2,-1/2,-3/-2,-1/6,... Sequence contains the denominators.

%C The mapping f(a/b) = (a + b)/(a - b). Taking a = 2 and b = 1 to start with and carrying out this mapping repeatedly on each new (reduced) rational number gives the periodic sequence 2/1,3/1,2/1,3/1,...

%H Seiichi Manyama, <a href="/A081478/b081478.txt">Table of n, a(n) for n = 1..26</a>

%F a(2n-1) = A007018(n-1), a(2n) = -A007018(n-1) for n >= 2. - _Jianing Song_, Oct 10 2021

%t Last /@ NestList[{(#1 - #2), #1 #2} & @@ # &, {2, 1}, 16] (* _Michael De Vlieger_, Sep 04 2016 *)

%o (Sage)

%o # Variant with first four terms slightly different. Absolute values.

%o def A081478_abs():

%o x, y = 1, 2

%o yield x

%o while True:

%o yield x

%o x, y = x * y, x//y + 1

%o a = A081478_abs(); print([next(a) for i in range(17)]) # _Peter Luschny_, Dec 17 2015

%Y A003687 gives the numerators.

%Y Cf. A007018.

%K sign,frac

%O 1,2

%A _Amarnath Murthy_, Mar 24 2003

%E More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Apr 06 2003