login
A372974
a(1) = 1, a(2) = 2; for n > 2, a(n) is the smallest unused positive number such that a(n) is coprime to a(n-1) and omega(a(n)) does not equal omega(a(n-1)).
3
1, 2, 15, 4, 21, 5, 6, 7, 10, 3, 14, 9, 20, 11, 12, 13, 18, 17, 22, 19, 24, 23, 26, 25, 28, 27, 34, 29, 30, 31, 33, 8, 35, 16, 39, 32, 45, 37, 36, 41, 38, 43, 40, 47, 42, 53, 44, 49, 46, 59, 48, 61, 50, 67, 51, 64, 55, 71, 52, 73, 54, 79, 56, 81, 58, 83, 57, 70, 69, 89, 60, 77, 78, 85, 66, 65, 84
OFFSET
1,2
COMMENTS
The fixed points show an unusual pattern; they begin 1, 2, 4, 69, 190, 438, 545, 725, 732, 909 and it appears, based on a graph of the sequence (see the attached image of the first 5000 terms) there may be no more. However more exist at 324388, 330574, 333069, 333531,..., 369752. Then once again there is a large gap until 2704713, 2726054, 2760963, ... . It is unclear what causes this behavior.
The sequence is conjectured to be a permutation of the positive integers.
LINKS
Scott R. Shannon, Image of the first 5000 terms. The green line is a(n) = n.
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^16, showing primes in red, perfect prime powers in gold, squarefree composites in green, and numbers neither squarefree nor prime powers in blue or purple, with purple indicating powerful numbers that are not prime powers.
EXAMPLE
a(3) = 15 as a(2) = 2 and omega(2) = A001221(2) = 1, and 15 is coprime to 2 while omega(15) = A001221(15) = 2 which does not equal 1. No smaller number satisfies both of these requirements.
MATHEMATICA
nn = 120; c[_] := False; Array[Set[{a[#], c[#]}, {#, True}] &, 2]; j = 2; u = 3;
Do[k = u;
While[Or[GCD[j, k] > 1, PrimeNu[k] == #, c[k]] &[PrimeNu[j]], k++];
Set[{a[n], c[k], j}, {k, True, k}];
If[k == u, While[c[u], u++]], {n, 3, nn}];
Array[a, nn] (* Michael De Vlieger, May 28 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Scott R. Shannon, May 26 2024
STATUS
approved