OFFSET
1,1
COMMENTS
No terms end in a 3, since if p == 3 (mod 10), then p^10 - p^9 - p^8 - p^7 - p^6 - p^5 - p^4 - p^3 - p^2 - p - 1 == 5 (mod 10) and is therefore not prime. - Michel Marcus, Jun 25 2014
LINKS
K. D. Bajpai, Table of n, a(n) for n = 1..11323
MATHEMATICA
Select[Prime[Range[2100]], PrimeQ[#^10-Total[#^Range[9]]-1]&] (* Harvey P. Dale, Sep 08 2019 *)
PROG
(Python)
import sympy
from sympy import isprime
{print(n, end=', ') for n in range(5*10**4) if isprime(n**10-n**9-n**8-n**7-n**6-n**5-n**4-n**3-n**2-n-1) and isprime(n)}
(PARI) for(n=1, 5*10^4, if(ispseudoprime(n)&&ispseudoprime(n^10-sum(i=0, 9, n^i)), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jun 04 2014
STATUS
approved