|
composite_small:=proc（n:：integer）
说明“确定n是否具有小于100的素因子”；
如果igcd（2305567963945518424753102147331756070，n）=1，则返回false，否则返回true-end if；
结束进程：
#长度为9的主星座模式
p:=[0,4,6,10,16,18,24,28,30]；
#使用isprime（m*n+o+p）
o：=[1273、2263、2683、4003、4633、4993、5893、6883、6943、8623、9613、10243、11563、12823、14863、15133、15553、17863、18433、19753、21163、21793、22483、23053、23113、24103、25783、27733、28723、29983]：
使用（ArrayTools）：
os：=尺寸（o，2）：
米：=30030:
活套停止：=10^11:
循环开始：=0:
打印（13）；
对于从loopstart到loopstop的n do
对于从1到os的do
计数器：=0；wc:=0；wd:=0；
而`和`（wd>-10，wd<9）确实
wd:=wd+1；
如果composite_small（m*n+o[a]+p[wd]）=false，则wd:=wd+1 else wd:=-10结束if
结束do；
如果wd>=9，则
而`和`（计数器>=0，wc<9）则
wc：=wc+1；如果是i素数（m*n+o[a]+p[wc]），则计数器：=计数器+1 else计数器：=-1 end if；
结束do；
结束条件：；
如果计数器=9，则打印（m*n+o[a]）结束；
结束do：
结束do：
| 