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!)
A221183 a(0)=0, a(1)=1; thereafter a(n) = gpf(2*a(n-1)+a(n-2)), where gpf = "greatest prime factor" (A006530). 1
0, 1, 2, 5, 3, 11, 5, 7, 19, 5, 29, 7, 43, 31, 7, 5, 17, 13, 43, 11, 13, 37, 29, 19, 67, 17, 101, 73, 19, 37, 31, 11, 53, 13, 79, 19, 13, 5, 23, 17, 19, 11, 41, 31, 103, 79, 29, 137, 101, 113, 109, 331, 257, 13, 283, 193, 223, 71, 73, 31, 5, 41, 29, 11, 17, 5, 3, 11, 5, 7, 19, 5, 29, 7, 43, 31, 7, 5, 17 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
Rapidly enters a loop of length 62: [5, 3, 11, 5, 7, 19, 5, 29, 7, 43, 31, 7, 5, 17, 13, 43, 11, 13, 37, 29, 19, 67, 17, 101, 73, 19, 37, 31, 11, 53, 13, 79, 19, 13, 5, 23, 17, 19, 11, 41, 31, 103, 79, 29, 137, 101, 113, 109, 331, 257, 13, 283, 193, 223, 71, 73, 31, 5, 41, 29, 11, 17].
LINKS
Greg Back and Mihai Caragiu, The Greatest Prime Factor and Recurrent Sequences, Fibonacci Quart. 48 (2010), no. 4, 358-362. [Discusses similar sequences]
PROG
(PARI) gpf(n) = if (n==1, 1, vecmax(factor(n)[, 1]));
lista(nn) = {print1(x=0, ", "); print1(y=1, ", "); for (n=2, nn, z = gpf(x+2*y); print1(z, ", "); x = y; y = z; ); } \\ Michel Marcus, Mar 01 2016
CROSSREFS
Sequence in context: A091809 A110315 A275677 * A178174 A094744 A229608
KEYWORD
nonn
AUTHOR
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 March 29 06:15 EDT 2024. Contains 371265 sequences. (Running on oeis4.)