OFFSET
1,2
COMMENTS
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
Wikipedia, Collatz Conjecture
MAPLE
collatz:= proc(n) option remember; `if`(n=1, 0,
1 + collatz(`if`(n::even, n/2, 3*n+1)))
end:
b:= proc() 0 end:
g:= proc(n) option remember; local t;
`if`(n=1, 0, g(n-1));
t:= collatz(n); b(t):= b(t)+1
end:
a:= proc(n) option remember; local k; for k
from 1+a(n-1) while g(k)>1 do od; k
end: a(0):=0:
seq(a(n), n=1..100);
MATHEMATICA
collatz[n_] := collatz[n] = If[n==1, 0,
1+collatz[If[EvenQ[n], n/2, 3n+1]]];
b[_] = 0;
g[n_] := g[n] = Module[{t}, If[n==1, 0, g[n-1]];
t = collatz[n]; b[t] = b[t]+1];
a[n_] := a[n] = Module[{k}, For[k = 1+a[n-1],
g[k] > 1, k++]; k]; a[0] = 0;
Array[a, 100] (* Jean-François Alcover, Jan 30 2021, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Alois P. Heinz, Jan 27 2021
STATUS
approved