login
A338977
Primes p*A007953(p)+1 for p in A338976.
1
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
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
Sequence in context: A140811 A247657 A241099 * A090686 A082277 A289154
KEYWORD
nonn,look
AUTHOR
J. M. Bergot and Robert Israel, Dec 18 2020
STATUS
approved