OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
MAPLE
f:= proc(p) local k;
k:= 1:
do
k:= nextprime(k);
if isprime(3*p*k*(3*p*k+1)-1) and
isprime(3*p*k*(3*p*k+1)+1) and
isprime(3*p*k*(3*p*k+3)-1) and
isprime(3*p*k*(3*p*k+3)+1) then return k
fi
od
end proc:
map(f, [seq(ithprime(i), i=1..50)]); # Robert Israel, Jan 04 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Apr 29 2008
EXTENSIONS
Corrected and extended by D. S. McNeil, Dec 10 2010
STATUS
approved