OFFSET
1,1
COMMENTS
Primes p such that ((1+sqrt(2))^p - (1-sqrt(2))^p) / (2*sqrt(2)) is a composite number.
PROG
(PARI) a000129(n) = ([2, 1; 1, 0]^n)[2, 1];
forprime(p=2, 1e3, if(!isprime(a000129(p)), print1(p, ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Mar 16 2016
STATUS
approved