OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
17 = 16 + 1 = 11[16] (in base 16), and 17 + 1 + 1 = 19 is again prime.
MAPLE
digsum:= (n, b) -> convert(convert(n, base, b), `+`):
select(p -> isprime(p) and isprime(p+digsum(p, 16)), [2, seq(i, i=3..1000, 2)]); # Robert Israel, Nov 07 2018
PROG
(PARI) forprime(p=1, 1999, isprime(p+sumdigits(p, 16))&&print1(p", "))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
M. F. Hasler, Nov 06 2018
STATUS
approved