login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A100113
a(n) = n if n <= 2, otherwise (smallest squarefree number m not occurring earlier such that gcd(m, a(n-1)) > 1).
4
1, 2, 6, 3, 15, 5, 10, 14, 7, 21, 30, 22, 11, 33, 39, 13, 26, 34, 17, 51, 42, 35, 55, 65, 70, 38, 19, 57, 66, 46, 23, 69, 78, 58, 29, 87, 93, 31, 62, 74, 37, 111, 102, 82, 41, 123, 105, 77, 91, 119, 85, 95, 110, 86, 43, 129, 114, 94, 47, 141, 138, 106, 53, 159, 165, 115, 130
OFFSET
1,2
COMMENTS
a(A100112(n)) and A100114(A100112(n)) define a pair of inverse permutations of the squarefree numbers: a(A100112(A100114(n))) = A100114(A100112(a(n))) = A005117(n);
A100115(n) = if n is squarefree then a(A100112(n)), otherwise n.
Comments from N. J. A. Sloane, Oct 29 2020 (Start)
An alternative definition is that this is the lexicographically earliest infinite sequence of distinct positive squarefree numbers with the property that gcd(a(n), a(n-1)) > 1 for n >= 3.
Described in this way, this is a squarefree version of the EKG sequence A064413, and it is easy to modify the proof that that sequence is a permutation of the positive integers so as to show that the present sequence is a permutation of the positive squarefree numbers, as claimed in the first comment.
Conjecture: With the three exceptions p = 2, 5, 13, and 31, when a prime p appears it is preceded by 2*p and followed by 3*p.
(End)
LINKS
Eric Weisstein's World of Mathematics, Squarefree
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Nov 07 2004
STATUS
approved