A362754 a(1) = 1, a(2) = 6; for n > 2, a(n) is the smallest positive number that has not yet appeared that shares a factor with a(n-1) and also contains as a factor the smallest prime that is not a factor of a(n-1).
1, 6, 10, 12, 15, 18, 20, 24, 30, 14, 21, 28, 36, 40, 42, 35, 50, 45, 48, 60, 56, 54, 70, 63, 66, 55, 22, 33, 44, 72, 75, 78, 65, 26, 39, 52, 84, 80, 90, 98, 96, 100, 102, 85, 34, 51, 68, 108, 105, 110, 99, 88, 114, 95, 38, 57, 76, 120, 112, 126, 130, 117, 104, 132, 135, 138, 115, 46, 69, 92, 144
Offset: 1
Keywords
Examples
a(3) = 10 as a(2) = 6 = 2*3, and 10 is the smallest unused number that shares a factor with 6 while also containing 5 as a prime factor, the smallest prime not a factor of 6. a(4) = 12 as a(3) = 10 = 2*5, and 12 is the smallest unused number that shares a factor with 10 while also containing 3 as a prime factor, the smallest prime not a factor of 10.
Links
- Scott R. Shannon, Table of n, a(n) for n = 1..10000
- Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^12, showing squarefree composites in green, numbers neither squarefree nor prime powers in blue, and numbers in A286708 in large light blue.
- Scott R. Shannon, Image of the first 250000 terms. The green line is a(n) = n.
- Scott R. Shannon, Image of the first 250000 terms in color. Terms with a lowest prime factor 2, 3, 5, 7, 11, >=13 are colored white, red, yellow, green, blue, violet and light gray respectively.
Programs
-
Mathematica
nn = 120; c[] := False; m[] := 2; MapIndexed[Set[{a[First[#2]], c[#1]}, {#1, True}] &, {1, 6}]; j = a[2]; Do[q = 2; While[Divisible[j, q], q = NextPrime[q]]; k = m[q]; While[Or[c[#], PrimePowerQ[#], CoprimeQ[j, k]] &[q k], k++]; k *= q; While[c[m[q] q], m[q]++]; Set[{a[n], c[k], j}, {k, True, k}], {n, 3, nn}]; Array[a, nn] (* Michael De Vlieger, May 02 2023 *)
Comments