OFFSET
1,1
COMMENTS
All terms are odd primes.
a(13) > 10^5.
MATHEMATICA
Select[Prime[Range[1, 100000]], PrimeQ[(15^# + 2^#)/17]&]
PROG
(PARI) forprime(p=3, 10^6, if(ispseudoprime((15^p + 2^p)/17), print1(p, ", ") ) ); \\ Joerg Arndt, Jul 29 2013
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Robert Price, May 07 2013
EXTENSIONS
Removed incorrect first term of "2".
STATUS
approved