OFFSET
0,5
COMMENTS
FORMULA
For n != 1, a(n) = A377435(n) - 1.
EXAMPLE
The perfect-powers in each prescribed range (rows):
.
.
.
9
25 27
36 49
81 100 121 125
144 169 196 216 225 243
289 324 343 361 400 441 484
529 576 625 676 729 784 841 900 961 1000
The binary expansions for n >= 3 (columns):
1001 11001 100100 1010001 10010000 100100001
11011 110001 1100100 10101001 101000100
1111001 11000100 101010111
1111101 11011000 101101001
11100001 110010000
11110011 110111001
111100100
MATHEMATICA
perpowQ[n_]:=n==1||GCD@@FactorInteger[n][[All, 2]]>1;
Table[Length[Select[Range[2^n+1, 2^(n+1)-1], perpowQ]], {n, 0, 15}]
PROG
(Python)
from sympy import mobius, integer_nthroot
def A377467(n):
def f(x): return int(1-sum(mobius(k)*(integer_nthroot(x, k)[0]-1) for k in range(2, x.bit_length())))
return f((1<<n+1)-1)-f((1<<n)) # Chai Wah Wu, Nov 05 2024
CROSSREFS
The version for squarefree numbers is A077643.
The version for prime-powers is A244508.
Including powers of 2 in the range gives A377435.
The version for non-perfect-powers is A377701.
The union of all numbers counted is A377702.
A081676 gives the greatest perfect-power <= n.
A131605 lists perfect-powers that are not prime-powers.
A377468 gives the least perfect-power > n.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Nov 04 2024
EXTENSIONS
a(26)-a(46) from Chai Wah Wu, Nov 05 2024
STATUS
approved