login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Lucky numbers that are pure (perfect) powers with exponent greater than 2.
1

%I #18 Nov 24 2021 04:11:56

%S 729,2187,29791,50625,50653,83521,130321,132651,166375,185193,226981,

%T 300763,707281,804357,823543,1295029,1860867,1874161,2685619,3176523,

%U 4782969,5545233,5764801,7189057,9393931,12977875,13845841,16194277

%N Lucky numbers that are pure (perfect) powers with exponent greater than 2.

%C Intersection of A000959 and A076467.

%H Kevin P. Thompson, <a href="/A139308/b139308.txt">Table of n, a(n) for n = 1..113</a> (terms < 10^10 with terms 1..94 from Giovanni Resta)

%e a(1) = 729 = 3^6 and 729 is lucky, i.e., A000959(118) = 729.

%t t = 2 Range@ 10000000 - 1; f[n_] := Block[{k = t[[n]]}, t = Delete[t, Table[{k}, {k, k, Length@t, k}]]]; Do[f@n, {n, 2, 1000000}]; fQ[n_] := ! IntegerQ@ Log[2, n] && GCD @@ Last /@ FactorInteger@n > 2; Select[t, fQ] (* _Robert G. Wilson v_, Oct 16 2010 *)

%Y Cf. A000959, A076467, A031162.

%K nonn

%O 1,1

%A _Walter Kehowski_, Jun 07 2008

%E a(16)-a(28) from _Robert G. Wilson v_, Oct 16 2010