cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A364246 a(1) = 1. Thereafter, a(n) is the least novel multiple of either prime(k+1) if rad(a(n-1)) = A002110(k), or Product_{prime q; q < gpf(a(n-1)); and q!|a(n-1)} q otherwise.

Original entry on oeis.org

1, 2, 3, 4, 6, 5, 12, 10, 9, 8, 15, 14, 30, 7, 60, 21, 20, 18, 25, 24, 35, 36, 40, 27, 16, 33, 70, 39, 770, 42, 45, 22, 105, 26, 1155, 28, 75, 32, 48, 50, 51, 10010, 54, 55, 84, 65, 462, 80, 57, 170170, 63, 90, 49, 120, 56, 135, 34, 15015, 38, 255255, 44, 210
Offset: 1

Views

Author

David James Sycamore, Jul 15 2023

Keywords

Comments

In other words, if the squarefree kernel of a(n-1) is a primorial term then a(n) is the least novel multiple of the smallest prime which does not divide a(n-1). Otherwise a(n) is the least novel multiple of the product of all primes < gpd(a(n-1)) which do not divide a(n-1). Primes >= 5 arrive late (as least unused term), and a(k) is prime(m) iff a(k-1) is A002110(m-1). The pattern around a prime is P(k), prime(k+1), 2*P(k), m*prime(k+1) for some multiplier m, where P(k) = A002110(k). The sequence is conjectured to be a permutation of the positive integers, with primes in natural order.
A common mode in this sequence is alternation of squarefree semiprime q(j)*q(k), j < k, followed by P(k-1)/q(j). The alternation often occurs in runs such that each iteration increments k. Example: a(241..246): q(2)*q(17) -> P(16)/q(2) -> q(2)*q(18) -> P(17)/q(2) -> q(2)*q(19) -> P(18)/q(2). a(16539..16572) represents a run of 17 alternations. - Michael De Vlieger, Jul 17 2023

Examples

			a(5) = 6 a primorial number so the next term is the smallest prime not dividing 6, thus a(7) = 5.
a(26) = 33 = 3*11 and the product of primes < 11 which do not divide 11 is 2*5*7 = 70, which has not occurred previously, therefore a(27) = 70.
		

Crossrefs

Programs

  • Mathematica
    nn = 120; c[] := False; m[] := 1; a[1] = j = 1; c[1] = True;
    f[x_] := If[# == Prime@ Range[PrimePi@ #[[-1]]], Prime[PrimePi@ #[[-1]] + 1],
      Times @@ Complement[Prime@ Range[PrimePi@ #[[-1]] - 1], #]] &[
      FactorInteger[x][[All, 1]]];
    Do[While[Set[k, f[j]]; c[k m[k]], m[k]++]; k *= m[k];
      Set[{a[n], c[k], j}, {k, True, k}], {n, 2, nn}];
    Array[a, nn] (* Michael De Vlieger, Jul 17 2023 *)

Extensions

More terms from Michael De Vlieger, Jul 17 2023

A364280 Lexicographically earliest sequence of distinct positive integers such that a(n) is the least novel multiple of m, the product of all primes q < gpf(a(n-2)*a(n-1)) which do not divide a(n-2)*a(n-1); a(1) = 1, a(2) = 2.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 10, 9, 8, 11, 105, 12, 13, 385, 18, 14, 15, 16, 17, 15015, 20, 19, 51051, 30, 21, 22, 25, 42, 23, 230945, 84, 24, 35, 26, 33, 70, 27, 28, 40, 36, 29, 37182145, 48, 31, 1078282205, 54, 32, 34, 30030, 37, 6678671, 60060, 38, 51, 5005, 44, 39
Offset: 1

Views

Author

David James Sycamore, Jul 17 2023

Keywords

Comments

It follows from the definition that the sequence is infinite.
Let r(n) = a(n-2)*a(n-1).
If rad(r(n)) is a primorial, then every prime q < gpf(r(n)) divides r(n), so m = 1, the empty product, and a(n) = u, the smallest missing number in the sequence so far.
If rad(r(n)) is not a primorial, then m > 1, and significant spikes can occur in scatterplot when there are multiple primes < gpf(r(n)) which do not divide r(n) (e.g., a(12) = 105, a(15) = 385, a(21) = 15015).
The only way a prime can occur is as u.
The sequence is a permutation of the positive integers since no number appears more than once and m = 1 eventually introduces any number not already placed consequent to terms arising from m > 1.

Examples

			a(4) = 4, a(5) = 5, and 3 is the only prime < 5 which does not divide 20, therefore m = 3 and a(6) = 6 since 3 has occurred once already.
a(10) = 8, a(11) = 11 and the product of all primes < 11 which do not divide 8*11 = 88 is 3*5*7 = 105, which has not occurred earlier, so a(12) = 105.
		

Crossrefs

Programs

  • Mathematica
    nn = 120; c[] := False; m[] := 1; a[1] = i = 1; a[2] = j = 2; c[1] = c[2] = True;
      f[x_] := Times @@ Complement[Prime@ Range[PrimePi@ #[[-1]] - 1], #] &[
        FactorInteger[x][[All, 1]]];
      Do[While[Set[k, f[i j]]; c[k m[k]], m[k]++]; k *= m[k];
        Set[{a[n], c[k], i, j}, {k, True, j, k}], {n, 3, nn}];
    Array[a, nn] (* Michael De Vlieger, Jul 17 2023 *)

Extensions

More terms from Michael De Vlieger, Jul 17 2023
Showing 1-2 of 2 results.