OFFSET
1,1
COMMENTS
All terms == 5 (mod 6).
Terms are not in numerical order.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 53 = 13*(1+3)+1 where 13 and 53 are prime.
MAPLE
A:= select(t -> isprime(t) and isprime(t*convert(convert(t, base, 10), `+`)+1), [$2..10^4]):
map(t -> t*convert(convert(t, base, 10), `+`)+1, A);
PROG
(PARI) f(p) = p*sumdigits(p)+1;
isok(p) = isprime(f(p));
apply(f, select(isok, primes(300))) \\ Michel Marcus, Dec 19 2020
CROSSREFS
KEYWORD
nonn,look
AUTHOR
J. M. Bergot and Robert Israel, Dec 18 2020
STATUS
approved