OFFSET
1,1
COMMENTS
All terms are even.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 22 is a member because A337183(22) = 37 and 37+22 = 59 are both prime.
MAPLE
f:= proc(n) local F, b, i;
F:= sort(map(t -> t[1]$t[2], ifactors(n)[2]), `>`);
b:= convert(F, `+`);
(add(F[i]*b^(i-1), i=1..nops(F)));
end proc:
filter:= proc(n) local v;
v:= f(n); isprime(v) and isprime(v+n)
end proc:
select(filter, 2*[$1..1000]);
MATHEMATICA
Select[Range[2, 1690], AllTrue[If[PrimeQ@ #, #, FromDigits[#, Total[#]] &@ Flatten[ConstantArray @@@ FactorInteger[#]]] + {0, #}, PrimeQ] &] (* Michael De Vlieger, Jan 29 2021 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
J. M. Bergot and Robert Israel, Jan 29 2021
STATUS
approved