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.

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

This page as a plain text file.
%I A153159 #10 Nov 22 2024 01:46:49
%S A153159 32,243,3125,7776,16807,100000,161051,248832,371293,537824,759375,
%T A153159 1419857,1889568,2476099,3200000,4084101,5153632,6436343,7962624,
%U A153159 11881376,17210368,20511149,24300000,28629151,39135393,45435424
%N A153159 a(n) = A007916(n)^5.
%t A153159 Select[Range[2,100],GCD@@Last/@FactorInteger@#==1&]^5
%o A153159 (Python)
%o A153159 from sympy import mobius, integer_nthroot
%o A153159 def A153159(n):
%o A153159     def f(x): return int(n+1-sum(mobius(k)*(integer_nthroot(x, k)[0]-1) for k in range(2, x.bit_length())))
%o A153159     m, k = n, f(n)
%o A153159     while m != k: m, k = k, f(k)
%o A153159     return m**5 # _Chai Wah Wu_, Nov 21 2024
%Y A153159 Cf. A007916, A153147, A153157, A153158, A153160, A113850.
%K A153159 nonn
%O A153159 1,1
%A A153159 _Vladimir Joseph Stephan Orlovsky_, Dec 19 2008
%E A153159 Edited and extended by _Ray Chandler_, Dec 22 2008