OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..150
EXAMPLE
MAPLE
revdigs:= proc(n) local L, k;
L:= convert(n, base, 10);
add(L[-k]*10^(k-1), k=1..nops(L))
end proc:
filter:= proc(n) local a, b;
a:= convert(map(convert, ifactors(n)[2], `*`), `+`);
if not isprime(a) then return false fi;
b:= convert(map(convert, ifactors(revdigs(n))[2], `*`), `+`);
b = a+2 and isprime(b)
end proc:
select(filter, [$1 .. 10^7]);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
J. M. Bergot and Robert Israel, Aug 12 2020
STATUS
approved