|
|
|
|
5, 23, 53, 137, 191, 827, 569, 1553, 857, 1091, 2087, 2339, 1811, 3347, 2711, 5867, 3449, 10979, 7127, 5231, 8753, 7883, 6131, 6311, 11057, 11633, 8111, 13649, 18899, 19661, 10091, 10631, 17393, 11171, 9209, 16703, 17207, 19727, 19979, 29741, 31181, 34739, 12809, 38699, 25523, 26783, 39461, 28547
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
All terms == 5 (mod 6).
Terms are not in numerical order.
|
|
LINKS
|
|
|
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
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|