login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A340585 Noncube perfect powers. 2
4, 9, 16, 25, 32, 36, 49, 81, 100, 121, 128, 144, 169, 196, 225, 243, 256, 289, 324, 361, 400, 441, 484, 529, 576, 625, 676, 784, 841, 900, 961, 1024, 1089, 1156, 1225, 1296, 1369, 1444, 1521, 1600, 1681, 1764, 1849, 1936, 2025, 2048, 2116, 2187, 2209, 2304, 2401, 2500 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
This was the original definition of A239870. However, the true definition of that sequence seems to be slightly different.
LINKS
FORMULA
Sum_{n>=1} 1/a(n) = 1 - zeta(3) + Sum_{k>=2} mu(k)*(1-zeta(k)) = 1 - A002117 + A072102 = 0.6724074652... - Amiram Eldar, Jan 12 2021
MAPLE
filter:= proc(n) local g;
g:= igcd(op(ifactors(n)[2][.., 2]));
g > 1 and (g mod 3 <> 0)
end proc:
select(filter, [$1..10000]); # Robert Israel, Jan 12 2021
MATHEMATICA
Select[Range[2, 2500], (g = GCD @@ FactorInteger[#][[;; , 2]]) > 1 && !Divisible[g, 3] &] (* Amiram Eldar, Jan 12 2021 *)
PROG
(PARI) for(n=2, 2500, if( ispower(n) % 3, print1(n, ", ")))
CROSSREFS
Sequence in context: A331219 A337534 A126589 * A010409 A010457 A244833
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Jan 12 2021
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 20:05 EDT 2024. Contains 371254 sequences. (Running on oeis4.)