OFFSET
1,1
COMMENTS
When A224610(i) = 1 then prime(i) is in this sequence.
Subsequence of A177104. - R. J. Mathar, Apr 19 2013
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..10000
MATHEMATICA
Reap[For[p = 2, p < 200000, p = NextPrime[p], If[PrimeQ[q = 2*p^3 - 1] && PrimeQ[r = 2*p*q^2 - 1], Sow[p]]]][[2, 1]] (* Jean-François Alcover, Apr 19 2013 *)
bpQ[n_]:=Module[{c=2n^3-1}, AllTrue[{c, 2n*c^2-1}, PrimeQ]]; Select[ Prime[ Range[ 15000]], bpQ] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Sep 05 2015 *)
PROG
(Magma) [p: p in PrimesUpTo(180000) | IsPrime(q) and IsPrime(2*p*q^2-1) where q is 2*p^3-1 ]; // Bruno Berselli, Apr 19 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Apr 12 2013
STATUS
approved