A101504 a(n+1)=a(n)/x+x, with x = Max(a(k): 0<=k<n and a(k)|a(n)), a(0)=a(1)=1.

1, 1, 2, 3, 4, 4, 5, 6, 5, 6, 7, 8, 6, 7, 8, 9, 6, 7, 8, 9, 10, 7, 8, 9, 10, 11, 12, 8, 9, 10, 11, 12, 13, 14, 9, 10, 11, 12, 13, 14, 15, 8, 9, 10, 11, 12, 13, 14, 15, 16, 10, 11, 12, 13, 14, 15, 16, 17, 18, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21

Offset

0,3

Comments

a(n)>a(n+1) iff a(n) is composite and a(n)>Max(a(k):0<=k<n);

a(A101505(n)) = n and a(m) < n for m < A101505(n).

Links

Ivan Neretin, Table of n, a(n) for n = 0..10000

Mathematica

Nest[Append[#, (x = Max@Intersection[Most@#, Divisors[#[[-1]]]]) + #[[-1]]/x] &, {1, 1}, 77] (* Ivan Neretin, Aug 24 2016 *)

Crossrefs

Sequence in context: A092032 A058222 A064064 * A125568 A248110 A108872

Adjacent sequences: A101501 A101502 A101503 * A101505 A101506 A101507

Keyword

nonn

Author

Reinhard Zumkeller, Dec 05 2004

Status

approved