|
PROG
|
(PARI) for(n=1, 7500, if(ispseudoprime(6^n-5^n-4^n-3^n-2^n-1), print(n)))
(Python)
from sympy import isprime
def afind(limit, k0=1):
pow6, pow5, pow4, pow3, pow2 = 6**k0, 5**k0, 4**k0, 3**k0, 2**k0
for k in range(k0, limit+1):
if isprime(pow6 - pow5 - pow4 - pow3 - pow2 - 1): print(k, end=", ")
pow6 *= 6; pow5 *= 5; pow4 *= 4; pow3 *= 3; pow2 *= 2
|