OFFSET
1,1
EXAMPLE
a(1) = 6 because 6*6 + 1 = 37 can be written as 5*7 + 2.
MAPLE
alist := proc(upto) local L, q, p, n, r; L := []; q := 2;
for n from 1 to upto do
p := q; q := nextprime(p); r := p * q + 1 ;
if modp(r, 6) = 0 and isprime(r + 1) then
L := [op(L), iquo(r, 6)] fi od;
L end: alist(350); # Peter Luschny, Jun 20 2021
MATHEMATICA
(Select[6Range[10^6]+1, PrimeQ[#] && MatchQ[FactorInteger[#-2], {{p_, 1}, {q_, 1}} /; q == NextPrime[p]]&]-1)/6 (* Jean-François Alcover, Jul 07 2021 *)
PROG
(PARI) a342564(plim)={my(p1=5); forprime(p2=7, plim, my(p=p1*p2+2); if(isprime(p), print1((p-1)/6, ", ")); p1=p2)};
a342564(2400)
(Python)
from primesieve.numpy import n_primes
from numbthy import isprime
primesarray = numpy.array(n_primes(10000, 1))
for i in range (0, 9999):
totest = int(primesarray[i] * primesarray[i+1] + 2)
if (isprime(totest)) and (((totest-1)%6) == 0):
print((totest-1)//6) # Karl-Heinz Hofmann, Jun 20 2021
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Hugo Pfoertner, Jun 20 2021
STATUS
approved