|
|
A224627
|
|
Prime numbers p such that 2*p^3-1, 2*p*q^2-1, 2*p*r^2-1, and 2*p*s^2-1 are prime numbers.
|
|
0
|
|
|
19460899, 86276401, 87980803, 167646631, 300722029, 343507111, 479516311, 906597943, 998757829, 1031308249, 1112697199, 1311383431, 1962194053
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
MATHEMATICA
|
Reap[ For[p = 2, p < 2*10^9, p = NextPrime[p], If[PrimeQ[q = 2*p^3 - 1] && PrimeQ[r = 2*p*q^2 - 1] && PrimeQ[s = 2*p*r^2 - 1] && PrimeQ[2*p*s^2 - 1], Print[p]; Sow[p]] ]][[2, 1]] (* Jean-François Alcover, Apr 22 2013 *)
|
|
PROG
|
(PFGW & SCRIPTIFY)
SCRIPT
DIM n, 1
DIM q
DIMS t
OPENFILEOUT myf, a(n).txt
LABEL a
SET n, n+1
SETS t, %d\,; p(n)
SET q, 2*p(n)^3-1
PRP q, t
IF ISPRP THEN GOTO b
GOTO a
LABEL b
SET q, 2*p(n)*q^2-1
PRP q, t
IF ISPRP THEN GOTO c
GOTO a
LABEL c
SET q, 2*p(n)*q^2-1
PRP q, t
IF ISPRP THEN GOTO d
GOTO a
LABEL d
SET q, 2*p(n)*q^2-1
PRP q, t
IF ISPRP THEN GOTO e
GOTO a
LABEL e
WRITE t, myf
GOTO a
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|