A283549 Composite numbers k such that tau(k^(k-1)) is a prime.
4, 9, 16, 27, 49, 64, 121, 125, 169, 289, 1681, 1849, 2401, 3481, 4913, 5329, 11881, 12769, 16129, 18769, 24649, 29791, 32041, 32761, 38809, 39601, 44521, 63001, 69169, 76729, 78125, 79507, 85849, 96721, 124609, 130321, 134689, 143641, 167281, 175561, 187489, 237169, 316969, 326041, 332929, 380689, 383161, 434281, 491401
Offset: 1
Keywords
Links
- Robert G. Wilson v, Table of n, a(n) for n = 1..10556
Programs
-
Mathematica
k = 1; lst = {}; While[k < 100001, If [ !PrimeQ@ k && PrimeQ[ DivisorSigma[0, k^(k -1)]], AppendTo[lst, k]]; k++]; lst (* or *) mx = 10^6; Union@ Flatten@ Reap[ Do[ Sow@ Select[ Prime[ Range[ PrimePi[ mx^(1/e) ]]]^e, PrimeQ[1 + e (#-1)] &], {e, 2, Log2[mx]}]][[2, 1]] (* Giovanni Resta, Mar 10 2017 *) -
PARI
is(n)=!isprime(n) && ispseudoprime(numdiv(n^(n-1))) \\ Charles R Greathouse IV, Mar 10 2017
Comments