OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
isA227907 := proc(n)
local q;
if isprime(n) then
q := nextprime(n) ;
isprime(4*n-q) ;
else
false;
end if;
end proc:
for n from 1 to 3000 do
if isA227907(n) then
printf("%d, ", n) ;
end if;
end do: # R. J. Mathar, Oct 14 2013
PROG
(PARI) isokp(p) = isprime(p) && isprime(4*p - nextprime(p+1)); \\ Michel Marcus, Oct 13 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Irina Gerasimova, Oct 12 2013
STATUS
approved