OFFSET
1,2
COMMENTS
Suggested by the formula a(a(n)) = 2*n + 3 for the analogous sequence A079000.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
MATHEMATICA
nn = 2^8; a[1] = 2; u = 4; v = {1}; w = {2}; p = 3;
Do[If[MemberQ[w, n], k = p;
w = Append[DeleteCases[w, n], p]; p = NextPrime[p],
If[Length[v] == 0,
k = u; AppendTo[w, u],
k = First[v]; v = Rest[v]]];
Set[{a[n]}, {k}];
If[n + 1 >= u, u++; While[PrimeQ[u], u++]], {n, 2, nn}];
{1}~Join~TakeWhile[Array[a[a[#]] &, nn], IntegerQ] (* Michael De Vlieger, Dec 17 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Dec 14 2024
STATUS
approved