OFFSET
1,1
COMMENTS
a(n) = prime(m) if A033955(m) is prime.
LINKS
Robert Israel, Table of n, a(n) for n = 1..4000
EXAMPLE
a(3) = 23 is a term because (23 mod 2) + ... + (23 mod 19) = 1+2+3+2+1+10+6+4 = 29 is prime.
MAPLE
f:= proc(n) local i, p;
p:= ithprime(n);
add(p mod ithprime(i), i=1..n-1)
end proc:
map(ithprime, select(t -> isprime(f(t)), [$1..2000]));
PROG
(PARI) isok(p) = if (isprime(p), my(s=0); forprime(q=2, precprime(p-1), s += p % q); isprime(s); ); \\ Michel Marcus, Jan 11 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Jan 11 2021
STATUS
approved