A257405 For n=1 or prime, a(n)=n; otherwise, a(n) is the smallest number not already used which shares a factor with a(n-1).

1, 2, 3, 6, 5, 10, 7, 14, 4, 8, 11, 22, 13, 26, 12, 9, 17, 34, 19, 38, 16, 18, 23, 46, 20, 15, 21, 24, 29, 58, 31, 62, 28, 30, 25, 35, 37, 74, 32, 36, 41, 82, 43, 86, 40, 42, 47, 94, 44, 33, 27, 39, 53, 106, 48, 45, 50, 52, 59, 118, 61, 122, 54, 51, 57, 60, 67

Offset

1,2

Comments

Since limsup(prime(n)) is infinity, it is easy to see that the sequence contains all even numbers. So it is natural to conjecture that the sequence is a permutation of the positive integers.

Links

Peter J. C. Moses, Table of n, a(n) for n = 1..1000

Mathematica

seq={1}; Do[cmplSeq=Complement[Range[2*Max[seq]], seq];
If[PrimeQ[n], AppendTo[seq, n], AppendTo[seq, Min[Select[cmplSeq, GCD[Last[seq], #]>1&]]]], {n, 2, 100}]; seq (* Ivan N. Ianakiev, Apr 25 2015 *)

Crossrefs

Cf. A000040, A064413, A257311, A257312, A257313, A257314, A257315.

Sequence in context: A080029 A091239 A145110 * A373545 A043547 A338702

Adjacent sequences: A257402 A257403 A257404 * A257406 A257407 A257408

Keyword

nonn

Author

Vladimir Shevelev, Apr 22 2015

Extensions

More terms from Peter J. C. Moses, Apr 22 2015

Status

approved