The On-Line Encyclopedia of Integer Sequences (OEIS)

A373545 revision #17

A373545

a(1) = 1, a(2) = 2; for n > 2, a(n) is the smallest unused positive number that shares a factor with a(n-1) if a(n-1) is odd otherwise is coprime to a(n-1) if a(n-1) is even.

1, 2, 3, 6, 5, 10, 7, 14, 9, 12, 11, 22, 13, 26, 15, 18, 17, 34, 19, 38, 21, 24, 23, 46, 25, 20, 27, 30, 29, 58, 31, 62, 33, 36, 35, 28, 37, 74, 39, 42, 41, 82, 43, 86, 45, 40, 47, 94, 49, 56, 51, 48, 53, 106, 55, 44, 57, 54, 59, 118, 61, 122, 63, 60, 67, 134, 65, 50, 69, 66, 71, 142, 73, 146

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

COMMENTS

The terms are concentrated along lines of different gradient; the odd terms all being on the line a(n) = n or very close to it, while the even terms are distributed among the other lines dependent on the number of prime factors they contain. The topmost line contains even semiprimes while the lower lines, those with gradient less than 1, contain powers of 2 multiplied by a larger single prime. However the topmost of these lower lines, which is more diffuse, also contains all other even numbers. Also noticeable is some of the lower lines terminating after which these values appear to move into the line contains all other even numbers.

Many odd terms are fixed points, this not occurring only when such a number would share a factor with the previous even number. This first occurs at a(65) = 67, when 65 cannot be chosen as it would share a factor with a(64) = 60.

The term selection rules would allow for consecutive odd numbers although this never occurs in the terms studied and is unlikely to ever occur. Likewise for the terms studied all primes appear in their natural order. No power of 2, other than 2 itself, can be a term.

From Michael De Vlieger, Jun 11 2024: (Start)

Odd prime p precedes 2*p.

Odd prime p appears for odd n, but additionally, tends to occupy a(n) such that n mod A002110(k) is a reduced residue of the largest A002110(k) < n. Thus, a(n) = p for n <= p. Example, a(65) = 67 is displaced. However, generally primes represent fixed points. (End)

LINKS

Scott R. Shannon, Table of n, a(n) for n = 1..10000

Michael De Vlieger, Plot f(a(n)) in rows of 210 terms, n = 1..16380, where color function f renders primes in red, perfect prime powers in gold, squarefree composites in greens, and numbers neither squarefree nor prime powers in blue or purple, where purple additionally signifies powerful numbers that are not prime powers. Bright green represents a primorial, yellow green an even squarefree semiprime, olive green an odd squarefree semiprime, and dark green other squarefree composites. Shows pattern of primes a(n) = p followed by even squarefree semiprimes 2*p, generally where n mod A002110(k) is a reduced residue. The lower image plots f(n), only showing primes in their place for comparison.

Scott R. Shannon, Image of the first 100000 terms. Numbers with one, two, three, four, or five and more prime factors, counted with multiplicity, are show as red, yellow, green, blue and violet respectively. The white line is a(n) = n.

Scott R. Shannon, Image of the first 100000 terms. As above except here the colors show the number of distinct prime factors.

EXAMPLE

a(7) = 7 as a(6) = 10 is an even number and 7 is the smallest unused positive that is coprime to 10.

MATHEMATICA

kk = 2; nn = 120; c[_] := False; Array[Set[{a[#], c[#]}, {#, True}] &, kk];

j = a[kk]; u = kk + 1;

Do[If[OddQ[j],

If[PrimePowerQ[j],

p = FactorInteger[j][[1, 1]];

k = #1 + Boole[#2 > 0] & @@ QuotientRemainder[u, p];

While[c[k p], k++]; k *= p,

k = u; While[Or[c[k], CoprimeQ[j, k]], k++]],