N： =2*10^7:#对于术语<=N
P： =选择（i素数，[2，seq（i，i=3…N/2，2）]）：
nP:=nops（P）：
SP:=0:
对于i从1开始，而P[i]^2<=N do
m： =列表工具：-二进制位置（P，N/P[i]）；
SP:=SP，op（P[i]*P[i.m]）；
日期：
SP:=排序（[SP]）：
SS:=ListTools:-部分和（SP）：
五： =矢量（6）：
SI：=矢量（6）：
二： =矢量（6，1）：
对于从1到6的i，执行SI[i]：=SS[i+1]-SS[1]od：
计数：=1:V[1]：=4:m:=4:im:={1}：
当计数<6 do时
对于im-do中的j
II[j]：=II[j]+1；
SI[j]：=不锈钢[II[j]+j]-不锈钢[II[j]]；
od；
m： =最小值（SI）；
im:=选择（j->SI[j]=m，{$1..20}）；
对于从1到20的k，而{$1..k}子集im-do
如果V[k]=0，则V[k]：=m；计数：=计数+1 fi
od；
日期：
转换（V，列表）；
