OFFSET
1,2
COMMENTS
If j=a(n-1) is in A046363 then a(n) = A001414(j) = prime q, and a(n+1) is the least novel multiple of q. Otherwise a(n) is the least novel multiple of the smallest prime factor of j. After a prime term a(n) the sequence produces a string of terms each divisible by the smallest prime factor of a(n+1) until arriving at a term in A046363, whereupon a new prime appears and the process repeats.
Conjectured to be a permutation of the positive integers A000027, in which the primes do not appear in order (prime order starts:2,5,7,3,13,11,19,17,31,23,43..).
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^12, showing primes in red, perfect prime powers in gold, squarefree composites in green, and numbers neither squarefree nor prime powers in blue and purple, with purple additionally representing powerful numbers that are not squarefree.
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^20, with the same color function as immediately above. Note various trajectories of primes.
EXAMPLE
a(2)=2 is not a term in A046360, and has smallest prime factor 2, so a(3) = 4, the least novel multiple of 2. Likewise a(4)=6 since a(3)=4 is not in A046360 and the smallest prime factor of 4 is 2.
a(6)=10 since 5 is the smallest prime factor of 5, and 10 is the smallest novel multiple of 5.
If a(n-1) = prime p, a(n) is the least novel multiple of p, for example a(12) = 3 and since a(4) = 6 it follows that a(13) = 9. Likewise a(19) = 13, and since no prior term is divisible by 13, a(20) = 36.
MATHEMATICA
nn = 120; c[_] := False; m[_] := 1;
a[1] = 1; j = a[2] = 2; c[1] = c[2] = True; m[1] = m[2] = 2;
f[x_] := f[x] = Total@ Flatten[ConstantArray[#1, #2] & @@@ FactorInteger[x]];
Do[(If[PrimeQ[#2], k = #2, k = #1]; While[c[k*m[k]], m[k]++]; k *= m[k]) & @@
{FactorInteger[j][[1, 1]], f[j]};
Set[{a[i], c[k], j}, {k, True, k}], {i, 3, nn}];
Array[a, nn] (* Michael De Vlieger, Sep 28 2024 *)
CROSSREFS
KEYWORD
nonn,new
AUTHOR
David James Sycamore, Sep 27 2024
STATUS
approved