A153159 a(n) = A007916(n)^5.

32, 243, 3125, 7776, 16807, 100000, 161051, 248832, 371293, 537824, 759375, 1419857, 1889568, 2476099, 3200000, 4084101, 5153632, 6436343, 7962624, 11881376, 17210368, 20511149, 24300000, 28629151, 39135393, 45435424

Offset

1,1

Links

Table of n, a(n) for n=1..26.

Mathematica

Select[Range[2, 100], GCD@@Last/@FactorInteger@#==1&]^5

Prog

(Python)
from sympy import mobius, integer_nthroot
def A153159(n):
    def f(x): return int(n+1-sum(mobius(k)*(integer_nthroot(x, k)[0]-1) for k in range(2, x.bit_length())))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m**5 # Chai Wah Wu, Nov 21 2024

Crossrefs

Cf. A007916, A153147, A153157, A153158, A153160, A113850.

Sequence in context: A351603 A343325 A343285 * A113850 A046454 A050997

Adjacent sequences: A153156 A153157 A153158 * A153160 A153161 A153162

Keyword

nonn

Author

Vladimir Joseph Stephan Orlovsky, Dec 19 2008

Extensions

Edited and extended by Ray Chandler, Dec 22 2008

Status

approved