OFFSET
1,1
LINKS
MAPLE
b:= proc(n, t) option remember; `if`(n<t, n,
b(`if`(n::even, n/2, 3*n+1), t))
end:
a:= n-> b((2*n+1)$2):
seq(a(n), n=1..80); # Alois P. Heinz, Jan 22 2022
MATHEMATICA
b[n_, t_] := b[n, t] = If[n<t, n, b[If[EvenQ[n], n/2, 3n+1], t]];
a[n_] := b[2n+1, 2n+1];
Table[a[n], {n, 1, 80}] (* Jean-François Alcover, Apr 25 2022, after Alois P. Heinz *)
PROG
(PARI) a(n) = my(N=2*n+1, m=N); while(m >= N, m = if (m%2, 3*m+1, m/2)); m; \\ Michel Marcus, Jan 22 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Apr 12 2001
EXTENSIONS
More terms from Larry Reeves (larryr(AT)acm.org), Apr 13 2001
STATUS
approved