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)).
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
Keywords
Examples
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.
Links
- Scott R. Shannon, Table of n, a(n) for n = 1..10000
- 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.
Programs
-
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 *)
Comments