%I #26 Jul 06 2024 01:31:10
%S 4,6,12,16,34,48,68,384,1080,5892,9816,34008,50034
%N Numbers k such that 6^k - 5^k - 4^k - 3^k - 2^k - 1 is prime.
%C a(11) > 7500.
%C All terms are even. - _Jon Perry_, Apr 12 2014
%e 6^4 - 5^4 - 4^4 - 3^4 - 2^4 - 1 = 317 is prime. Thus, 4 is a term.
%o (PARI) for(n=1,7500,if(ispseudoprime(6^n-5^n-4^n-3^n-2^n-1),print(n)))
%o (Python)
%o from sympy import isprime
%o def afind(limit, k0=1):
%o pow6, pow5, pow4, pow3, pow2 = 6**k0, 5**k0, 4**k0, 3**k0, 2**k0
%o for k in range(k0, limit+1):
%o if isprime(pow6 - pow5 - pow4 - pow3 - pow2 - 1): print(k, end=", ")
%o pow6 *= 6; pow5 *= 5; pow4 *= 4; pow3 *= 3; pow2 *= 2
%o afind(1100) # _Michael S. Branicky_, Aug 21 2021
%Y Cf. A147977.
%Y Cf. A000043, A240503.
%K nonn,more
%O 1,1
%A _Derek Orr_, Apr 06 2014
%E a(11) from _Michael S. Branicky_, Aug 21 2021
%E a(12)-a(13) from _Michael S. Branicky_, Jul 05 2024
|