A050138 a(1)=2, a(2)=6. For n >= 2, a(n) = floor(a(n-1)/2) if this is not among 0,a(1),...,a(n-1); otherwise a(n) = 3*n.

2, 6, 3, 1, 15, 7, 21, 10, 5, 30, 33, 16, 8, 4, 45, 22, 11, 54, 27, 13, 63, 31, 69, 34, 17, 78, 39, 19, 9, 90, 93, 46, 23, 102, 51, 25, 12, 114, 57, 28, 14, 126, 129, 64, 32, 138, 141, 70, 35, 150, 75, 37, 18, 162, 81, 40, 20, 174, 87, 43

Offset

1,1

Comments

Does this sequence contain every positive integer exactly once?

Inverse: 4, 1, 3, 14, 9, 2, 6, 13, 29, 8, 17, 37, 20, 41, 5, 12, 25, 53, ..., . - Robert G. Wilson v, Apr 09 2018

Links

Robert Israel, Table of n, a(n) for n = 1..10000

Maple

S:= {0, 2, 6}: A[1]:= 2: A[2]:= 6:
for n from 3 to 100 do
  t:= floor(A[n-1]/2);
  if member(t, S) then t:= 3*n fi;
  A[n]:= t;
  S:= S union {t};
od:
seq(A[n], n=1..100); # Robert Israel, Apr 09 2018

Mathematica

f[s_] := Block[{b = Floor[s[[-1]]/2], l = Length@ s}, Append[s, If[MemberQ[s, b], 3l, b]]]; s = {0, 2, 6}; Nest[f, s, 57] (* Robert G. Wilson v, Apr 09 2018 *)

Crossrefs

Cf. A050000, A050128, A050135, A050136, A050136, A050138.

Sequence in context: A021042 A136694 A164104 * A108443 A201761 A011042

Adjacent sequences: A050135 A050136 A050137 * A050139 A050140 A050141

Keyword

nonn,look

Author

Clark Kimberling

Extensions

Name corrected by Robert Israel, Apr 09 2018

Status

approved