The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A353917 a(1) = 4. Let j = a(n-1) and let m = omega(gcd(j, k)) with gcd(j, k) > 1. For n > 1, a(n) = least k such that min(omega(j), omega(k)) = m and m < max(omega(j), omega(k)), but neither j | k nor k | j. 3
4, 6, 8, 10, 16, 12, 9, 15, 25, 20, 30, 18, 27, 21, 49, 14, 32, 22, 64, 24, 42, 28, 70, 40, 60, 36, 66, 44, 110, 50, 90, 48, 78, 52, 128, 26, 169, 39, 81, 33, 121, 55, 125, 35, 343, 56, 84, 54, 102, 68, 170, 80, 120, 45, 105, 63, 168, 72, 114, 76, 190, 100, 130 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
The sequence exhibits phases involving alternating composite prime powers and squarefree semiprimes. These manifest in log-log scatterplot in a caustic fashion, where the composite prime power is very much larger than the squarefree semiprime for sufficiently large n.
Let P = the set of distinct prime divisors of j = a(n-1), and let Q = the set of distinct prime divisors of k = a(n). Let g = gcd(j, k) > 1 and let G = {P intersect Q}. Noncoprime j and k implies |G| > 0. This sequence is such that |G| > 0, |P| > |G|, and |Q| == |G|, or vice versa, yet neither j | k nor k | j.
Theorem: primes are prohibited. Proof: since we have gcd(j, k) > 1 and do not allow divisibility, and since primes must either divide or be coprime to another number m, primes do not appear in this sequence.
Theorem: squarefree semiprimes j = pq are followed by k = p^2 or k = q^2. Proof: since omega(j) = |P| = 2 and is squarefree, we have 2 cases pertaining to successor k, both with gcd(j, k) > 1.
1.) |P| == |G| implies |Q| > |G| and |Q| > |P|.
2.) |Q| == |G| implies |P| > |G| and |P| > |Q|.
The first case implies some prime r | k yet gcd(j, r) = 1. But this would require j | k, which is prohibited. The second case suggests either p | k or q | k, but so as to satisfy non-divisibility axiom, we are forced into either k = p^e, e > 1, or k = q^m, m > 1.
Corollary: powers of the same prime appear in natural order in this sequence.
There is a weaker alternation between numbers in A120944 and A350352 as n is sufficiently large. This alternation exhibits prime power factor features akin to the composite prime power-squarefree semiprime alternation.
Conjecture: permutation of composite numbers.
LINKS
Michael De Vlieger, Annotated log log scatterplot of a(n), n = 1..1192, showing records in red and local minima in blue.
Michael De Vlieger, Log-log scatterplot of a(n), n = 1..2^16, showing records in red, local minima in blue, accentuating composite prime powers in green and squarefree semiprimes in gold.
Michael De Vlieger, Prime power factor study of a(n), n = 1..10^4 where n increases from left to right, and pi(p) from bottom to top in the top image. Multiplicity e of p^e is represented by black for e = 1, red for e = 2, and a color function to the maximum multiplicity in the range thereafter. On the lower bar we show composite prime powers in yellow, squarefree semiprimes in orange, numbers in A350352 in green, and those in A126706 in blue.
MATHEMATICA
nn = 2^7; c[_] = 0; j = a[1] = 4; c[4] = 1; u = 6; Do[Set[k, u]; Set[m, PrimeNu[j]]; While[Nand[c[k] == 0, ! Divisible[#2, #1] & @@ Sort[{j, k}], And[#2 > #3, #1 == #3] & @@ Append[Sort[{m, PrimeNu[k]}], PrimeNu[GCD[j, k]]]], k++]; Set[{a[i], c[k]}, {k, i}]; j = k; If[k == u, While[Nand[c[u] == 0, CompositeQ@ u], u++]], {i, 2, nn}]; Array[a, nn]
CROSSREFS
Sequence in context: A288814 A062711 A280739 * A117347 A088011 A289424
KEYWORD
nonn
AUTHOR
Michael De Vlieger, May 10 2022
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 16 13:17 EDT 2024. Contains 372552 sequences. (Running on oeis4.)