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