SELECT 1
   USE  d:\最密三生素数文件\六生成元17.DBF ALIAS 六生成元17
   SELECT 2
   USE  d:\最密三生素数文件\六生成元19.DBF ALIAS 六生成元19
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 19
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 385
           && @ 4,8 say 孪素
           X=scy17+510510*(N-1)
           MZ=X%19
                IF  MZ=0 OR MZ=4 OR MZ=6 OR MZ=10 OR MZ=12 OR MZ=16
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE scy19 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）的生成元17，制作生成元19（上一个台阶）

SELECT 1
   USE  d:\最密三生素数文件\六生成元19.DBF ALIAS 六生成元19
   SELECT 2
   USE  d:\最密三生素数文件\六生成元23.DBF ALIAS 六生成元23
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 23
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 5005
           && @ 4,8 say 孪素
           X=scy19+9699690*(N-1)
           MZ=X%23
                IF  MZ=0 OR MZ=4 OR MZ=6 OR MZ=10 OR MZ=12 OR MZ=16
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE scy23 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）的生成元19，制作生成元23（上一个台阶）

SELECT 1
   USE  d:\最密三生素数文件\六生成元23.DBF ALIAS 六生成元23
   SELECT 2
   USE  d:\最密三生素数文件\六生成元29.DBF ALIAS 六生成元29
    kssj=SECONDS()                      &&取出开始时间
    For N=1 to 29
    @ 3,6 say N
       SELECT 1
       go 1
       for m=1 to 85085
           && @ 4,8 say 孪素
           X=scy23+223092870*(N-1)
           MZ=X%29
                IF  MZ=0 OR MZ=4 OR MZ=6 OR MZ=10 OR MZ=12 OR MZ=16
                  else
                  SELECT   2      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE scy29 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）的生成元23，制作生成元29（上一个台阶），最后生成元29有19569550个元素

SELECT 1
  USE d:\最密三生素数文件\素数表.DBF ALIAS 素数表
  SELECT 2
  USE d:\最密三生素数文件\六生成元29.DBF ALIAS 六生成元29
  SELECT 3
  USE d:\最密三生素数文件\对称6素.DBF ALIAS 对称6素表
     kssj=SECONDS()                      &&取出开始时间
  For i=1 to 2
         @ 3,6 say i
         SELECT 2
       go 1
       SELECT 2
    for m=1 to 1956955
     @ 6,12 say m
       k=scy29
       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=4  OR qmz=6 OR qmz=10 OR qmz=12 OR qmz=16
                        EXIT         
                ENDIF
                SKIP            &&素数表指针向下移动一个
        ENDFOR
        IF j>sss-10
                SELECT  3  &&打开盛放素数的表
                APPEND BLANK                        &&增加一条空记录
                REPLACE dc6s 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,"运行时间提示")
用最密对称6生成元29制作6生素数

SELECT 1
  USE d:\最密三生素数文件\素数表.DBF ALIAS 素数表
  SELECT 2
  USE d:\最密三生素数文件\六生成元29.DBF ALIAS 六生成元29
  SELECT 3
  USE d:\最密三生素数文件\对称6素.DBF ALIAS 对称6素表
     kssj=SECONDS()                      &&取出开始时间
  For i=15 to 16
         @ 3,6 say i
         SELECT 2
       go 1
       SELECT 2
    for m=1 to 1956955
     @ 6,12 say m
       k=scy29
       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=4  OR qmz=6 OR qmz=10 OR qmz=12 OR qmz=16
                        EXIT         
                ENDIF
                SKIP            &&素数表指针向下移动一个
        ENDFOR
        IF j>sss-10
                SELECT  3  &&打开盛放素数的表
                APPEND BLANK                        &&增加一条空记录
                REPLACE dc6s 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,"运行时间提示")
用生成元制作最密6生素数