OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..2300
EXAMPLE
MAPLE
spf:= n -> add(t[1]*t[2], t=ifactors(n)[2]):
filter:= proc(p) local v;
if not isprime(p) then return false fi;
v:= p mod spf(p-1);
isprime(v) and p mod spf(p+1) = v
end proc:
select(filter, [seq(i, i=3..10^5, 2)]);
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Nov 26 2020
STATUS
approved