OFFSET
1,1
COMMENTS
Primes not in A059801. - Robert Israel, Nov 03 2024
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= p -> isprime(p) and not isprime(4^p-3^p):
select(filter, [seq(i, i=3..1000, 2)]); # Robert Israel, Nov 03 2024
MATHEMATICA
Select[Prime[Range[70]], CompositeQ[4^#-3^#]&] (* Harvey P. Dale, Mar 14 2025 *)
PROG
(PARI) apmb(a, b, n) = { forprime(x=2, n, y=a^x-b^x; if(!ispseudoprime(y), print1(x", "); ) ) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, Oct 26 2003
EXTENSIONS
Offset corrected by Mohammed Yaseen, Jul 18 2022
STATUS
approved
