OFFSET
1,2
COMMENTS
A057130 gives the product of prime numbers (-1 mod 6) in the order of occurrence.
EXAMPLE
PROG
(PFGW Script)
SCRIPT
DIM i, 0
DIM j
DIM k
DIM n, 1
OPENFILEOUT myf, a(n).txt
OPENFILEIN maf, a007528.txt
LABEL a
SET i, i+1
IF i>100 THEN END
GETNEXT j, maf
SET n, n*j
SET k, 0
LABEL b
SET k, k+1
PRP k*6+1, k
IF ISPRP THEN GOTO c
GOTO b
LABEL c
PRP k*6*n-1, k
IF ISPRP THEN GOTO d
GOTO b
LABEL d
WRITE myf, k
GOTO a
(PARI) lista(nn) = {my(pp = 1); forprime (p = 1, nn, if (Mod(p, 6) == -1, pp *= p; my(k=1); while (!isprime(6*k*pp-1) || !isprime(6*k*pp+1), k++); print1(k, ", "); ); ); } \\ Michel Marcus, Nov 25 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Nov 24 2019
STATUS
approved