キュバン素数

number-pedia.com

キュバン素数を求めるプログラム。十進BASIC使用。最初の1000個を求めています。

LET s=0
LET t=0
FOR a=1 TO 10000
LET b=3*a^2+3*a+1
FOR c=2 TO INT(SQR(b))
IF MOD(b,c)=0 THEN EXIT FOR
IF MOD(b,c)<>0 THEN LET s=s+1
NEXT C
IF s=INT(SQR(b))-1 THEN LET t=t+1
IF s=INT(SQR(b))-1 THEN PRINT t; b;"prime"
IF t=1000 THEN stop
LET s=0
NEXT A

END