OFFSET
1,1
COMMENTS
I suppose that there are infinitely many twin primes in the sequence, e.g., (239,241), (809,811), (1877,1879).
EXAMPLE
a(3)=83 since 83+(8+3)+(8^2+3^2) = 83+11+73 = 167 is a prime.
PROG
(PARI) isok(n) = {if (isprime(n), digs = digits(n, 10); isprime(n + sum(i=1, #digs, digs[i]) + sum(i=1, #digs, digs[i]^2)); , 0; ); } \\ Michel Marcus, Jul 18 2013
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Carmine Suriano, Jul 21 2010
STATUS
approved