|
|
A302034
|
|
A028234 analog for a factorization process based on the Ludic sieve (A255127); Discard all instances of the (smallest) Ludic factor A272565(n) from n.
|
|
7
|
|
|
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, 19, 13, 1, 27, 1, 7, 7, 29, 11, 15, 1, 31, 13, 1, 11, 33, 1, 17, 5, 35, 1, 9, 1, 37, 17, 19, 1, 39, 7, 5, 11, 41, 1, 21, 1, 43, 35, 11, 1, 45, 1, 23, 1, 47, 13, 3, 1, 49, 23, 25, 1, 51, 13, 13, 19
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,6
|
|
COMMENTS
|
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 sequence of distinct Ludic numbers (A003309) in ascending order, while applying A302035 to the same terms gives the corresponding "exponents" of these Ludic factors in this nonstandard "Ludic factorization of n", 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 and examples in A302032.
|
|
LINKS
|
|
|
FORMULA
|
|
|
PROG
|
(PARI)
\\ Assuming A269379 and its inverse A269380 have been precomputed, then the following is reasonably fast:
A302034(n) = if(1==n, n, my(k=0); while((n%2), n = A269380(n); k++); n = (n/2^valuation(n, 2)); while(k>0, n = A269379(n); k--); (n));
|
|
CROSSREFS
|
Cf. A302036 (gives the positions of 1's).
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|