OFFSET
1,2
COMMENTS
Is this (apart from the first 2 entries) the same as A006885? - R. J. Mathar, Aug 16 2020
Yes, a(n) = A006885(n) for n >= 3. - A.H.M. Smeets, Nov 13 2025
LINKS
A.H.M. Smeets, Table of n, a(n) for n = 1..97
PROG
(PARI) firstMiss(A) = { my(i); if(#A == 0 || A[1] > 0, return(0)); for(i = 1, A[#A] + 1, if(!setsearch(A, i), return(i))); };
iter(A) = { my(a = firstMiss(A)); while(!setsearch(A, a), A = setunion(A, Set([a])); a = if(a % 2, 3*a+1, a/2)); A; };
makeVec(m) = { my(v = [], A = Set([0]), i); for(i = 1, m, if (length(v) == 0 || v[length(v)] != A[#A], v = concat(v, A[#A])); if (i < m, A = iter(A))); v; };
makeVec(10000)
CROSSREFS
KEYWORD
nonn
AUTHOR
Markus Sigg, Aug 10 2020
EXTENSIONS
a(16)-a(22) from Jinyuan Wang, Aug 13 2020
a(23)-a(27) from A006885 - A.H.M. Smeets, Nov 15 2025
STATUS
approved
