login
A366103
a(n) = n if n is 1 or a prime. If n is composite, p is the greatest prime not exceeding n, and k = n - p, then a(n) is the least novel m*a(k), where m is not a prior term.
1
1, 2, 3, 4, 5, 6, 7, 8, 18, 27, 11, 9, 13, 10, 24, 36, 17, 12, 19, 14, 30, 45, 23, 15, 32, 48, 64, 80, 29, 16, 31, 20, 42, 63, 84, 105, 37, 21, 44, 66, 41, 22, 43, 25, 52, 78, 47, 26, 56, 99, 112, 140, 53, 28, 68, 102, 132, 165, 59, 33, 61, 34, 70, 114, 152, 175
OFFSET
1,2
COMMENTS
There are infinitely many primes, and given any nonprime positive integer n there exists a greatest prime p such that p < n. Therefore the sequence is well defined, infinite, and the "least novel m*a(k)" condition ensures that every number will eventually be included. The sequence is therefore a permutation of A000027, in which (by the definition) the primes are in natural order.
LINKS
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^20, highlighting prime p = a(p) in red.
Michael De Vlieger, Log log scatterplot of a(n), n = 1..1024, showing primes in red, composite prime powers in gold, squarefree composites in green, and numbers neither squarefree nor prime powers in blue, and highlighting powerful numbers that are not prime powers in light blue.
Michael De Vlieger, Color-coded list of terms 30 to a row, 900 total, showing primes in red and squarefree semiprimes in light green, demonstrating a tendency for the former to directly precede the latter. The primes maintain their arrangement therefore outside of 2, 3, and 5, they tend to fall within reduced residues mod 30.
FORMULA
a(n) = n for n prime.
EXAMPLE
a(1) = 1, a(2) = 2 and a(3) = 3 since 1,2 and 3 are all in A008578 (1 and primes).
For n = 4 write 4 = 3+1 then a(4) is the least novel multiple m*a(1) = m, where m is not a term; thus a(4) = 4.
a(5) = 5 since 5 is prime
For n = 6 write 6 = 5+1 then a(6) = least novel multiple m*a(1) = m, where m is not a term, so a(6) = 6.
a(7) = 7 since 7 is prime.
For n = 8 write 8 = 7+1 then a(8) = 8
For n = 9 write 9 = 7+2 then a(9) = least novel m*a(2) = 2*m, where m is not yet a term, thus a(9) = 2*9 = 18.
Sequence can be listed as an irregular table with each row starting with a term in A008578 (1 or prime), in which row lengths equate to the corresponding gaps between consecutive primes; A001223.
1;
2;
3, 4;
5, 6;
7, 8, 18, 27;
11, 9;
13, 10, 24, 36;
17, 12;
19, 14, 30, 45;
23, 15, 32, 48, 64, 80;
...
MATHEMATICA
nn = 1000;
c[_] := False; m[_] := 1;
Monitor[Do[If[CompositeQ[n], k = n - Prime@ PrimePi[n];
While[Or[c[m[k]*k], c[m[k]]], m[k]++]; k *= m[k],
k = n];
Set[{a[n], c[k]}, {k, True}], {n, nn}], n];
Array[a, nn]
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
STATUS
approved