login
A356014
Consider the exponents in the prime factorization of n, and replace each run of k consecutive e's by a unique e; the resulting list corresponds to the exponents in the prime factorization of a(n).
3
1, 2, 3, 4, 3, 2, 3, 8, 9, 10, 3, 12, 3, 10, 3, 16, 3, 18, 3, 20, 21, 10, 3, 24, 9, 10, 27, 20, 3, 2, 3, 32, 21, 10, 3, 4, 3, 10, 21, 40, 3, 10, 3, 20, 45, 10, 3, 48, 9, 50, 21, 20, 3, 54, 21, 40, 21, 10, 3, 12, 3, 10, 63, 64, 21, 10, 3, 20, 21, 10, 3, 72, 3
OFFSET
1,2
COMMENTS
We ignore the exponents (all 0's) for the prime numbers beyond the greatest prime factor of n.
This sequence operates on prime exponents as A090079 and A337864 operate on binary and decimal digits, respectively.
FORMULA
a(a(n)) = a(n).
a(n^k) = a(n)^k for any k >= 0.
a(n) = A319521(A356008(n)).
A007814(a(n)) = A007814(n).
a(n) = 3 iff n belongs to A294674 \ {1}.
a(n) = 4 iff n belongs to A061742 \ {1}.
a(n) = 8 iff n belongs to A115964.
EXAMPLE
For n = 99:
- 99 = 11^1 * 7^0 * 5^0 * 3^2 * 2^0,
- the list of exponents is: 1 0 0 2 0,
- compressing consecutive values, we obtain: 1 0 2 0,
- so a(99) = 7^1 * 5^0 * 3^2 * 2^0 = 63.
PROG
(PARI) a(n) = { my (v=1, e=-1, k=0); forprime (p=2, oo, if (n==1, return (v), if (e!=e=valuation(n, p), v*=prime(k++)^e); n/=p^e)) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Jul 23 2022
STATUS
approved