login
Odd powers that are not prime powers.
1

%I #24 Feb 16 2025 08:33:18

%S 225,441,1089,1225,1521,2025,2601,3025,3249,3375,3969,4225,4761,5625,

%T 5929,7225,7569,8281,8649,9025,9261,9801,11025,12321,13225,13689,

%U 14161,15129,16641,17689,18225,19881,20449,21025,21609,23409,24025,25281,25921,27225

%N Odd powers that are not prime powers.

%C Numbers in A075109 but not in A000961.

%C Also odd perfect powers having no primitive root (intersection of A075109 and A175594).

%H Arkadiusz Wesolowski, <a href="/A216419/b216419.txt">Table of n, a(n) for n = 1..10000</a>

%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/PrimitiveRoot.html">Primitive Root</a>.

%F Sum_{n>=1} 1/a(n) = 1/2 + Sum_{k>=2} mu(k)*(1-zeta(k)*(2^k-1)/2^k) - Sum_{p prime} 1/(p*(p-1)) = 0.0158808884... - _Amiram Eldar_, Dec 21 2020

%e 81 = 9^2 as well as 81 = 3^4, therefore 81 is not a term.

%e 225 can be expressed so in one way as (3*5)^2, therefore 225 is a term.

%t nn = 27500; lst = Union[Flatten[Table[n^i, {i, Prime[Range[PrimePi[Log[2, nn]]]]}, {n, 2, nn^(1/i)}]]]; Select[lst, OddQ[#] && ! IntegerQ@PrimitiveRoot[#] &] (* Most of the code is from _T. D. Noe_ *)

%o (Magma) [n : n in [3..27225 by 2] | IsPower(n) and EulerPhi(n) ne CarmichaelLambda(n)]; // _Arkadiusz Wesolowski_, Nov 09 2013

%Y Cf. A008683, A075109, A136141, A175594.

%K nonn,changed

%O 1,1

%A _Arkadiusz Wesolowski_, Sep 06 2012