A153160 a(n) = A007916(n)^6.

64, 729, 15625, 46656, 117649, 1000000, 1771561, 2985984, 4826809, 7529536, 11390625, 24137569, 34012224, 47045881, 64000000, 85766121, 113379904, 148035889, 191102976, 308915776, 481890304, 594823321, 729000000, 887503681

Offset

1,1

Links

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

Mathematica

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

Prog

(Python)
from sympy import mobius, integer_nthroot
def A153160(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**6 # Chai Wah Wu, Nov 21 2024

Crossrefs

Cf. A007916, A153147, A153157, A153158, A153159, A113851.

Sequence in context: A351604 A343286 A240845 * A113851 A046455 A092758

Adjacent sequences: A153157 A153158 A153159 * A153161 A153162 A153163

Keyword

nonn

Author

Vladimir Joseph Stephan Orlovsky, Dec 19 2008

Extensions

Edited and extended by Ray Chandler, Dec 22 2008

Status

approved