A175723 a(1)=a(2)=1; thereafter a(n) = gpf(a(n-1)+a(n-2)), where gpf = "greatest prime factor".

1, 1, 2, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5, 2, 7, 3, 5

Offset

1,3

Comments

Rapidly enters a loop with period 3,5,2,7.

More generally, if a(1) and a(2) are distinct positive numbers with a(1)+a(2) >= 2, the sequence eventually enters the cycle {7,3,5,2} [Back and Caragiu].

Links

Table of n, a(n) for n=1..121.

G. Back and M. Caragiu, The greatest prime factor and recurrent sequences, Fib. Q., 48 (2010), 358-362.

Mathematica

nxt[{a_, b_}]:={b, FactorInteger[a+b][[-1, 1]]}; Transpose[NestList[nxt, {1, 1}, 120]][[1]] (* or *) PadRight[{1, 1, 2}, 130, {5, 2, 7, 3}] (* Harvey P. Dale, Feb 24 2015 *)

Crossrefs

Cf. A000045, A006530, A020639.

Similar or related sequences: A177904, A177923, A178094, A178095, A178174, A178179, A180101, A180107, A221183.

Sequence in context: A178094 A364874 A122556 * A084346 A165911 A354764

Adjacent sequences: A175720 A175721 A175722 * A175724 A175725 A175726

Keyword

nonn

Author

N. J. A. Sloane, Dec 16 2010

Status

approved