login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A081460 Consider the mapping f(r) = (1/2)*(r + N/r) from rationals to rationals where N = 5. Starting with a = 2 and applying the mapping to each new (reduced) rational number gives 2, 9/4, 161/72, 51841/23184, ..., tending to N^(1/2). Sequence gives values of the denominators. 8
1, 4, 72, 23184, 2403763488, 25840354427429161536, 2986152136938872067784669198846010266752, 39878504028822311675150039382403961856254569551519724209276629577579916539865344 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
Related sequence pairs (numerator, denominator) can be obtained by choosing N = 2, 3, 6, etc.
The sequence satisfies the Pell equation A081459(n+1)^2 - 5*a(n+1)^2 = 1. - Vincenzo Librandi, Dec 20 2011
LINKS
FORMULA
a(n) = 2*a(n-1)*A081459(n-1). - Mario Catalani (mario.catalani(AT)unito.it), May 21 2003
a(n) = A000045(A007283(n-1))/2. - Ehren Metcalfe, Oct 07 2017
From Amiram Eldar, Apr 07 2023: (Start)
a(n) = A079613(n-1)/2.
a(n) = Product_{k=1..n-1} L(3*2^(n-1-k)), where L(k) is the k-th Lucas number (A000032). (End)
a(n) = A001076(2^(n-1)). - Robert FERREOL, Apr 18 2023
MATHEMATICA
Table[Fibonacci[2^(n - 1)*3], {n, 1, 8}]/2 (* Amiram Eldar, Apr 07 2023 *)
PROG
(PARI) {r=2; N=5; for(n=1, 8, a=numerator(r); b=denominator(r); print1(b, ", "); r=(1/2)*(r + N/r))}
(Magma) m:=8; f:=[ n eq 1 select 2 else (Self(n-1)+5/Self(n-1))/2: n in [1..m] ]; [ Denominator(f[n]): n in [1..m] ]; // Bruno Berselli, Dec 20 2011
CROSSREFS
Sequence in context: A344693 A172478 A087315 * A327040 A327112 A284673
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Mar 22 2003
EXTENSIONS
Edited and extended by Klaus Brockhaus and Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Apr 06 2003
a(8) corrected by Vincenzo Librandi, Dec 20 2011
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 25 11:36 EDT 2024. Contains 373701 sequences. (Running on oeis4.)