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-6 of 6 results.

A302031 An omega (A001221) analog based on the Ludic sieve (A255127): a(1) = 0; for n > 1, a(n) = 1 + a(A302034(n)).

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 2, 3, 1, 2, 2, 1, 2, 3, 1, 2, 2, 3, 1, 2, 1, 2, 2, 2, 1, 3, 2, 2, 2, 2, 1, 2, 1, 2, 3, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 3, 2, 2, 2
Offset: 1

Views

Author

Antti Karttunen, Apr 02 2018

Keywords

Crossrefs

Cf. A302036 (positions of terms < 2).
Differs from similar A302041 for the first time at n=59, where a(59) = 2, while A302041(59) = 1.

Programs

Formula

a(1) = 0; for n > 1, a(n) = 1 + a(A302034(n)).
a(n) = A001221(A302026(n)).
a(n) = A069010(A269388(n)).

A302044 A028234 analog for factorization process based on the sieve of Eratosthenes (A083221).

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 1, 5, 1, 3, 1, 7, 5, 1, 1, 9, 1, 5, 1, 11, 1, 3, 1, 13, 7, 7, 1, 15, 1, 1, 5, 17, 7, 9, 1, 19, 11, 5, 1, 21, 1, 11, 1, 23, 1, 3, 1, 25, 25, 13, 1, 27, 1, 7, 7, 29, 1, 15, 1, 31, 13, 1, 11, 33, 1, 17, 5, 35, 1, 9, 1, 37, 17, 19, 11, 39, 1, 5, 11, 41, 1, 21, 7, 43, 35, 11, 1, 45, 1, 23, 1, 47, 13, 3, 1, 49, 19, 25, 1, 51, 1, 13, 25
Offset: 1

Views

Author

Antti Karttunen, Mar 31 2018

Keywords

Comments

Iterating n, a(n), a(a(n)), a(a(a(n))), ..., until 1 is reached, and taking the smallest prime factor (A020639) of each term gives a sequence of distinct primes in ascending order, while applying A302045 to the same terms gives the corresponding exponents (multiplicities) of those primes. Permutation pair A250245/A250246 maps between this non-standard prime factorization and the ordinary factorization of n. See also comments and examples in A302042.

Crossrefs

Programs

  • PARI
    \\ Assuming A250469 and its inverse A268674 have been precomputed, then the following is fast enough:
    A302044(n) = if(1==n,n,my(k=0); while((n%2), n = A268674(n); k++); n = (n/2^valuation(n, 2)); while(k>0, n = A250469(n); k--); (n));
    
  • PARI
    up_to = 65537;
    ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
    A000265(n) = (n/2^valuation(n, 2));
    A020639(n) = if(n>1, if(n>n=factor(n, 0)[1, 1], n, factor(n)[1, 1]), 1); \\ From A020639
    v078898 = ordinal_transform(vector(up_to,n,A020639(n)));
    A078898(n) = v078898[n];
    A302044(n) = { my(c = A000265(A078898(n))); if(1==c,1,my(p = prime(-1+primepi(A020639(n))+primepi(A020639(c))), d = A078898(c), k=0); while(d, k++; if((1==k)||(A020639(k)>=p),d -= 1)); (k*p)); };

Formula

For n > 1, a(n) = A250469^(r)(A000265(A078898(n))), where r = A055396(n)-1 and A250469^(r)(n) stands for applying r times the map x -> A250469(x), starting from x = n.
a(n) = A250245(A028234(A250246(n))).

A302032 Discard the least ludic factor of n: a(n) = A255127(A260738(c) + r - 1, A260739(c)), where r = A260738(n), c = A260739(n) are the row and the column index of n in the table A255127; a(n) = 1 if c = 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 3, 5, 1, 6, 1, 7, 5, 8, 1, 9, 5, 10, 9, 11, 1, 12, 1, 13, 7, 14, 1, 15, 7, 16, 15, 17, 7, 18, 1, 19, 11, 20, 1, 21, 1, 22, 21, 23, 1, 24, 19, 25, 19, 26, 1, 27, 11, 28, 27, 29, 11, 30, 1, 31, 13, 32, 11, 33, 1, 34, 33, 35, 1, 36, 13, 37, 17, 38, 1, 39, 35, 40, 39, 41, 1, 42, 31, 43, 35, 44, 1, 45, 1, 46, 45, 47, 13, 48, 1, 49, 23, 50
Offset: 1

Views

Author

Antti Karttunen, Mar 31 2018

Keywords

Comments

Original definition: A032742 analog for a nonstandard factorization process based on the Ludic sieve (A255127); Discard a single instance of the Ludic factor A272565(n) from n.
Like [A020639(n), A032742(n)] or [A020639(n), A302042(n)], also ordered pair [A272565(n), a(n)] is unique for each n. Iterating n, a(n), a(a(n)), a(a(a(n))), ..., until 1 is reached, and taking the Ludic factor (A272565) of each term gives a multiset of Ludic numbers (A003309) in ascending order, unique for each natural number n >= 1. Permutation pair A302025/A302026 maps between this "Ludic factorization" and the ordinary prime factorization of n. See also comments in A302034.
The definition of "discard the least ludic factor" is based on the table A255127 of the ludic sieve, where row r lists the (r+1)-th ludic number k = A003309(r+1), determined at the r-th step of the sieve, followed by the numbers crossed out at this step, namely, every k-th of the numbers remaining so far after k. If the number n is in row r = A260738(n), column c = A260739(n) of that table, then its least ludic factor is A272565(n) = A003309(r+1), the 1st entry of the r-th row. To discard that factor means to consider the number which is r-1 rows below the number c in that table, whence a(n) = A255127(A260738(c)+r-1, A260739(c)) - unless n is a ludic number, in which case a(n) = 1. - M. F. Hasler, Nov 06 2024

Examples

			Frem _M. F. Hasler_, Nov 06 2024: (Start)
For ludic numbers 1, 2, 3, 5, 7, 11, 13, 17, 23, 25, 29, 37, ..., a(n) = 1.
For n = 4, an even number, we have r = A260738(4) = 1: It is listed in row 1 of the table A255127, which lists all numbers that were crossed out at the first step: namely, the ludic number k = 2 and every other larger number. Also, in this row 1, the number 4 is in column c = A260739(4) = 2. Therefore, we apply r-1 = 0 times the map A269379 to c = 2, whence a(4) = 2.
The number n = 6 is also even and therefore listed in row r = 1, now in column c = 3, whence a(6) = 3. Similarly, a(8) = 4 and a(2k) = k for all k >= 1.
The number n = 9 was crossed out at the 2nd step (so r = A260738(9) = 2), when k = 3 was added to the ludic numbers and every 3rd remaining number crossed out; 9 was the first of these (after k = 3) so it is in column c = A260739(9) = 2. Now we have to apply r-1 = 1 times the map A269379 to c. That map yields the number which is located just below the argument (here c = 2) in the table A255127. Since 2 is a ludic number, in the first column, we get the next larger ludic number, 3, whence a(9) = 3.
The number 15 was the (c = 3)rd number to be crossed out at the (r = 2)nd step. Hence a(15) = A269379^{r-1} (c) = A269379(3) = 5 (again, the next larger ludic number).
The number 19 was the (c = 2)nd number to be crossed out at the (r = 3)rd step (when k = 5, its least ludic factor, was added to the list of ludic numbers). Hence a(19) = A269379^2(2) = A269379(3) = 5 again (skipping twice to the next larger ludic number).
(End)
To illustrate how this sequence allows one to compute the complete "ludic factorization" of a number, we consider n = 100.
For n = 100, its Ludic factor A272565(100) is 2, and we have seen that a(n) = 100/2 = 50.
For n = 50, its Ludic factor A272565(50) is 2 again, and again a(50) = 50/2 = 25.
Since n = 25 = A003309(1+9) is a ludic number, it equals its Ludic factor A272565(25) = 25. Because it appeared at the A260738(25) = 9th step, we apply A269379 eight times to the column index A260739(25) = 1, a fixed point, so a(25) = A269379^8(1) = 1.
Collecting the Ludic factors given by A272565 we get the multiset of factors: [2, 2, 25] = [A003309(1+1), A003309(1+1), A003309(1+9)]. By definition, A302026(100) = prime(1)*prime(1)*prime(9) = 2*2*23 = 92, the product of the corresponding primes.
If we start from n = 100, iterating the map n -> A302034(n) [instead of A302032] and apply A272565 to each term obtained we get just a single instance of each Ludic factor: [2, 25]. Then by applying A302035 to the same terms we get the corresponding exponents (multiplicities) of those factors: [2, 1].
		

Crossrefs

Cf. the following analogs A302031 (omega), A302037 (bigomega).
Cf. also A032742, A302042.

Programs

  • PARI
    \\ Assuming A269379 and its inverse A269380 have been precomputed, then the following is reasonably fast:
    A302032(n) = if(1==n,n,my(k=0); while((n%2), n = A269380(n); k++); n = n/2; while(k>0, n = A269379(n); k--); (n))

Formula

For n > 1, a(n) = A269379^r'(A260739(n)), where r' = A260738(n)-1 and A269379^r'(n) stands for applying r' times the map x -> A269379(x), starting from x = n.
a(n) = A302025(A032742(A302026(n))).
From M. F. Hasler, Nov 06 2024: (Start)
a(n) = 1 if n is a ludic number, i.e., in A003309. Otherwise:
a(n) = A255127(A260738(c) + r - 1, A260739(c)), with r = A260738(n), c = A260739(n).
In particular, a(2n) = n for all n. (End)

A302026 Permutation of natural numbers mapping "Ludic factorization" to ordinary factorization: a(1) = 1, a(2n) = 2*a(n), a(A269379(n)) = A003961(a(n)).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 25, 20, 27, 22, 19, 24, 23, 26, 21, 28, 29, 30, 49, 32, 45, 34, 35, 36, 31, 50, 33, 40, 37, 54, 41, 44, 81, 38, 43, 48, 125, 46, 75, 52, 47, 42, 121, 56, 63, 58, 77, 60, 53, 98, 39, 64, 55, 90, 59, 68, 135, 70, 61, 72, 169, 62, 51, 100, 67, 66, 175, 80, 99, 74, 71, 108, 343, 82, 105, 88
Offset: 1

Views

Author

Antti Karttunen, Apr 03 2018

Keywords

Comments

See comments and examples in A302032 to see how Ludic factorization proceeds.

Crossrefs

Cf. A302025 (inverse permutation).
Cf. A005940, A250246, A269172, A269388 (similar or related permutations).

Programs

Formula

a(1) = 1, a(2n) = 2*a(n), a(2n+1) = A003961(a(A269380(2n+1))).
a(n) = A250246(A269172(n)).
a(n) = A005940(1+A269388(n)).
Other identities. For all n >= 1:
A001221(a(n)) = A302031(n).
A001222(a(n)) = A302037(n).

A302036 Ludic powers: numbers k such that A302031(k) < 2; numbers k such that A260739(k) is a power of 2.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 21, 23, 25, 29, 31, 32, 37, 41, 43, 45, 47, 49, 53, 55, 61, 64, 67, 71, 73, 77, 83, 85, 89, 91, 93, 97, 101, 107, 109, 115, 119, 121, 127, 128, 131, 143, 145, 149, 151, 157, 161, 167, 173, 175, 179, 181, 189, 191, 193, 197, 205, 209, 211, 221, 223, 227, 229, 233, 235, 239, 247, 253, 256, 257
Offset: 1

Views

Author

Antti Karttunen, Apr 02 2018

Keywords

Comments

An analog of A000961 for factorization process based on the Ludic sieve (A255127).
Numbers k for which A302031(k) < 2, or equally, for which A302034(k) = 1, or equally, for which A209229(A260739(k)) = 1.

Crossrefs

Cf. A000079, A003309, A254100 (subsequences).
Cf. also A000961, A302038, A302040.

Programs

  • PARI
    for(n=1,257,if(A302031(n)<2,print1(n,","))); \\ See also code in A302031.

A302035 a(1) = 0, for n > 1, a(n) = A001511(A260739(n)); Number of instances of (the smallest) Ludic factor A272565(n) in n.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 4, 1, 1, 2, 2, 3, 1, 1, 3, 1, 1, 1, 2, 1, 1, 2, 5, 2, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 4, 1, 1, 4, 3, 1, 1, 2, 1, 1, 2, 3, 2, 1, 1, 2, 1, 1, 1, 6, 1, 1, 1, 2, 3, 1, 1, 3, 2, 1, 1, 2, 1, 1, 2, 4, 2, 1, 1, 2, 3, 1, 1, 3, 1, 1, 1, 2, 5, 1, 1, 5, 1, 1, 1, 2, 2, 1, 1, 3, 2
Offset: 1

Views

Author

Antti Karttunen, Apr 01 2018

Keywords

Comments

An A067029 analog for "Ludic factorization": iterating the map n -> A302034(n) until 1 is reached, and taking the Ludic factor (A272565) of each term gives a sequence of distinct Ludic numbers (A003309) in ascending order, while applying this function (A302035) to those terms gives the corresponding "exponents" of those Ludic factors, that is, the count of consecutive occurrences of each when iterating the map n -> A302032(n), which gives the same factors with repetitions. Permutation pair A302025/A302026 maps between the Ludic factorization and the ordinary prime factorization of n. See also comments and examples in A302032.

Crossrefs

Formula

a(1) = 0; for n > 1, a(n) = A001511(A260739(n)).
For n > 1, a(n) = A302025(A067029(A302026(n))).
Showing 1-6 of 6 results.