SELECT 1
   USE D:\最密三生素数文件\素数表.DBF ALIAS 素数表
   SELECT 2
   USE D:\最密三生素数文件\最密9生素数.DBF ALIAS 密9中
   kssj=SECONDS()
   SELECT 1
   GO 1
   For i=1 to 50876228
    @ 5,12 say i
   SELECT 1 &&打开盛放素数的表
   jl=recno()
   PUBLIC A
   A=素数
   SKIP 4
   PUBLIC B
   B=素数
   SKIP 4
   PUBLIC C
   C=素数
   cj1=B-A &&求前后两个素数的差值
   cj2=C-B &&求前后两个素数的差值
     IF cj1=12 AND cj2=18
     SELECT 2
     APPEND BLANK  &&增加一条空记录
     REPLACE jss WITH A+15 &&将A+4付给变量szx
     ENDIF
     SELECT 1
     GO jl+1
     ENDFOR
   =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
用素数表筛选最密9生素数中项0,2,6,8,12,18,20,26,30

  SELECT 1
   USE  d:\最密三生素数文件\对称8素13.DBF ALIAS 对称8素13
   SELECT 2
    USE  d:\最密三生素数文件\对称8素末.DBF ALIAS 对称8素11
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 13
    @ 3,6 say N
       SELECT 2
       go 1
       for m=1 to 8
           && @ 4,8 say 孪素
           X=dc8s+2310*(N-1)
           MZ=X%13
                IF  MZ=0 OR MZ=1 OR MZ=2 OR MZ=6 OR MZ=7 OR MZ=11 OR MZ=12
                  else
                  SELECT   1      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE dc8s13 WITH X     &&将N值付给素
                ENDIF
                SELECT 2
                SKIP
       endfor
    ENDFOR
   =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
制作最密对称8生元（从11升级到13）

SELECT 1
   USE  d:\最密三生素数文件\对称8素13.DBF ALIAS 对称8素13
   SELECT 2
   USE  d:\最密三生素数文件\对称8素17.DBF ALIAS 对称8素17
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 17
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 48
           && @ 4,8 say 孪素
           X=dc8s13+30030*(N-1)
           MZ=X%17
                IF  MZ=0 OR MZ=2 OR MZ=3 OR MZ=6 OR MZ=7 OR MZ=9 OR MZ=12 OR MZ=14
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE dc8s17 WITH X     &&将N值付给素
                ENDIF
                SELECT 2
                SKIP
       endfor
    ENDFOR
   =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
制作最密对称8生元（从13升级到17）

SELECT 1
   USE  d:\最密三生素数文件\对称8素17.DBF ALIAS 对称8素17
   SELECT 2
   USE  d:\最密三生素数文件\对称8素19.DBF ALIAS 对称8素19
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 19
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 432
           && @ 4,8 say 孪素
           X=dc8s17+510510*(N-1)
           MZ=X%19
                IF  MZ=0 OR MZ=1 OR MZ=2 OR MZ=5 OR MZ=6 OR MZ=7 OR MZ=12 OR MZ=14
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE dc8s19 WITH X     &&将N值付给素
                ENDIF
                SELECT 1
                SKIP
       endfor
    ENDFOR
   =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
制作最密对称8生元（从17升级到19）

SELECT 1
   USE  d:\最密三生素数文件\对称8素23.DBF ALIAS 对称8素23
   SELECT 2
   USE  d:\最密三生素数文件\对称8素29.DBF ALIAS 对称8素29
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 29
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 71280
           && @ 4,8 say 孪素
           X=dc8s23+223092870*(N-1)
           MZ=X%29
                IF  MZ=0 OR MZ=2 OR MZ=6 OR MZ=12 OR MZ=14 OR MZ=20 OR MZ=24 OR MZ=26
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE dc8s29 WITH X     &&将N值付给素
                ENDIF
                SELECT 1
                SKIP
       endfor
    ENDFOR
   =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
制作最密对称8生元（从23升级到29）

CLEAR ALL
SELECT 1
USE d:\最密三生素数文件\素数表.DBF ALIAS 素数表
SELECT 2
USE d:\最密三生素数文件\对称8素29.DBF ALIAS 对称8素
SELECT 3
USE d:\最密三生素数文件\对称8素表.DBF ALIAS 对称8素表
     kssj=SECONDS()                      &&取出开始时间
  For i=1 to 2
         @ 3,6 say i
         SELECT 2
       go 1
       SELECT 2
    for m=1 to 1496880
       k=dc8s29
       N=k+6469693230*(i-1)
        Kf=INT(SQRT(N))       &&求出奇数的开方根
        SELECT  1             &&打开素数表
        LOCATE FOR 素数>=kf      &&根据开方根,查找最大素数
        DO CASE
                CASE EOF()          &&如果超出素数表最后一条记录
                        GO BOTTOM
                CASE 素数>kf     &&如果找到的素数比开方根值大,就将记录指针退回一个,保证使用的最大素数在开方根内
                        SKIP -1
        ENDCASE
        SSS=RECNO()    &&记录下最大素数（开方根内）记录位置（如果找到的素数等于开方根，则直接读取这个素数的位置）
        GO 11           &&从第一条记录开始读取素数(3)
        FOR j=1 TO SSS-10      &&内循环开始。这个循环实质上是从小到大顺序,依次读取素数。循环值是记录序号
                qmz=MOD(N,素数)   &&以读取的素数为条件,对被判断值求模 【qmz=MOD(bpz,素数)】文本号括着的原在此行语句体部分.特别需要注意的是在用标点符号时，一定改成英文模式及半角状态。
                IF qmz=0 OR qmz=2  OR MZ=6 OR MZ=12 OR MZ=14 OR MZ=20 OR MZ=24 OR MZ=26
                        EXIT         
                ENDIF
                SKIP            &&素数表指针向下移动一个
        ENDFOR
        IF j>sss-10
                SELECT  3  &&打开盛放素数的表
                APPEND BLANK                        &&增加一条空记录
                REPLACE dc8s WITH N             &&将N值付给素数
        ENDIF
        SELECT 2
        skip
      ENDFOR
   ENDFOR
        =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
用生成元制作最密8生素数

CLEAR ALL
SELECT 1
USE d:\最密三生素数文件\素数表.DBF ALIAS 素数表
SELECT 2
USE d:\最密三生素数文件\对称8素29.DBF ALIAS 对称8素
SELECT 3
USE d:\最密三生素数文件\对称8素表.DBF ALIAS 对称8素表
     kssj=SECONDS()                      &&取出开始时间
  For i=1 to 2
         @ 3,6 say i
         SELECT 2
       go 1
       SELECT 2
    for m=1 to 1496880
     @ 6,12 say m
       k=dc8s29
       N=k+6469693230*(i-1)
        Kf=INT(SQRT(N))       &&求出奇数的开方根
        SELECT  1             &&打开素数表
        LOCATE FOR 素数>=kf      &&根据开方根,查找最大素数
        DO CASE
                CASE EOF()          &&如果超出素数表最后一条记录
                        GO BOTTOM
                CASE 素数>kf     &&如果找到的素数比开方根值大,就将记录指针退回一个,保证使用的最大素数在开方根内
                        SKIP -1
        ENDCASE
        SSS=RECNO()    &&记录下最大素数（开方根内）记录位置（如果找到的素数等于开方根，则直接读取这个素数的位置）
        GO 11           &&从第一条记录开始读取素数(3)
        FOR j=1 TO SSS-10      &&内循环开始。这个循环实质上是从小到大顺序,依次读取素数。循环值是记录序号
                qmz=MOD(N,素数)   &&以读取的素数为条件,对被判断值求模 【qmz=MOD(bpz,素数)】文本号括着的原在此行语句体部分.特别需要注意的是在用标点符号时，一定改成英文模式及半角状态。
                IF qmz=0 OR qmz=2  OR qmz=6 OR qmz=12 OR qmz=14 OR qmz=20 OR qmz=24 OR qmz=26
                        EXIT         
                ENDIF
                SKIP            &&素数表指针向下移动一个
        ENDFOR
        IF j>sss-10
                SELECT  3  &&打开盛放素数的表
                APPEND BLANK                        &&增加一条空记录
                REPLACE dc8s WITH N             &&将N值付给素数
        ENDIF
        SELECT 2
        skip
      ENDFOR
   ENDFOR
        =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
上楼有点瑕疵

SELECT 1
  USE D:\最密三生素数文件\素数表.dbf ALIAS 素数表
  kssj=SECONDS()  &&取出开始时间
  s=1.000000000000000000
  SELECT 1
    GO 10
    A=素数
    s=s*(A*(A-16)/(A-8)^2)
        SKIP 1
        FOR  j=1 TO 50876228
        SELECT 1
        B=素数
        s=s*(B*(B-16)/(B-8)^2)
        SELECT 1  
        SKIP
        ENDFOR

 ?s
 =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
对称最密8生素数的中项和合成，求其系数（公式中，公共系数）

SELECT 1
   USE  d:\最密三生素数文件\六生成元11.DBF ALIAS 六生成元11
   SELECT 2
   USE  d:\最密三生素数文件\六生成元13.DBF ALIAS 六生成元13
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 13
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 5
           && @ 4,8 say 孪素
           X=scy11+2310*(N-1)
           MZ=X%13
                IF  MZ=0 OR MZ=4 OR MZ=6 OR MZ=10 OR MZ=12 OR MZ=3
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE scy13 WITH X     &&将N值付给素
                ENDIF
                SELECT 1
                SKIP
       endfor
    ENDFOR
   =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
用最密6生素数（0,4,2,4,2,4）的生成元11，制作生成元13（上一个台阶）

SELECT 1
   USE  d:\最密三生素数文件\六生成元13.DBF ALIAS 六生成元13
   SELECT 2
   USE  d:\最密三生素数文件\六生成元17.DBF ALIAS 六生成元17
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 17
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 35
           && @ 4,8 say 孪素
           X=scy13+30030*(N-1)
           MZ=X%17
                IF  MZ=0 OR MZ=4 OR MZ=6 OR MZ=10 OR MZ=12 OR MZ=16
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE scy17 WITH X     &&将N值付给素
                ENDIF
                SELECT 1
                SKIP
       endfor
    ENDFOR
   =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
用最密6生素数（0,4,2,4,2,4）的生成元13，制作生成元17（上一个台阶）