OFFSET
1,2
COMMENTS
A non divisor prime of a(n-1) is any prime p < Gpf(a(n-1)) which does not divide a(n-1). A007947(a(n-1)) is in A002110 iff a(n-1) is a term in A055932. Sequence is conjectured to be a permutation of the natural numbers (A000027) with primes in order.
Scatterplot shows trajectories of numbers whose smallest prime factor is prime p, e.g., for p = 5, numbers in A084967, p = 7, those in A084968, p = 11 those in A084969, etc. - Michael De Vlieger, Oct 09 2024
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^20.
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^14, 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 the latter additionally representing powerful numbers that are not prime powers.
EXAMPLE
a(1) = 1 = A002110(0) so a(2) = 2 (smallest prime not already a term).
a(2) = 2 = A002110(1) so a(3) = 3.
a(3) = 3 not a term in A002110 so a(4) is least novel multiple of 2, the least non divisor prime of 3. Therefore a(4) = 4 since 2 has occurred earlier.
a(39) = 42, not a term in A002110 so a(40) = 25, the least novel multiple of 5, the smallest non divisor prime of 42.
MATHEMATICA
nn = 120; c[_] := False; m[_] := 1; f[x_] := FactorInteger[x][[All, 1]];
Array[Set[{a[#], c[#], m[#]}, {#, True, 2}] &, 2]; j = 2; v = 3;
Do[If[Or[IntegerQ@ Log2[j], And[EvenQ[j], Union@ Differences@ PrimePi[#] == {1}]],
k = v; While[c[k*m[k]], m[k]++]; k *= m[k],
k = 2; While[Divisible[j, k], k = NextPrime[k]];
While[c[k*m[k]], m[k]++]; k *= m[k]] &[f[j]];
Set[{a[n], c[k], j}, {k, True, k}];
If[k == v, While[c[v], v = NextPrime[v]]], {n, 3, nn}];
Array[a, nn] (* Michael De Vlieger, Oct 09 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
David James Sycamore, Oct 06 2024
EXTENSIONS
More terms from Michael De Vlieger, Oct 09 2024
STATUS
approved