OFFSET
1,5
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..7725
MAPLE
A210651 := proc(n)
for k from 1 by 2 do
if isprime(3*k*2^n-1) then
return k;
end if;
end do:
end proc: # R. J. Mathar, Mar 29 2012
PROG
PFGW64 and SCRIPTIFY from Primeform group
Command: pfgw64 -f in.txt
in.txt file :
SCRIPT
DIM nn, 0
DIM kk
DIMS tt
OPENFILEOUT myfile, a(n).txt
LABEL loopn
SET nn, nn+1
IF nn>7725 THEN END
SET kk, -1
LABEL loopk
SET kk, kk+2
SETS tt, %d, %d\,; nn; kk
PRP 3*kk*2^nn-1, tt
IF ISPRP THEN GOTO a
IF ISPRIME THEN GOTO a
GOTO loopk
LABEL a
WRITE myfile, tt
GOTO loopn
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Mar 27 2012
STATUS
approved