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

?s
=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 8
  A=素数
  s=s*(A^11*(A-12)/(A-1)^12)
        SKIP 1
        FOR  j=1 TO 50876228
        SELECT 1
        B=素数
        s=s*(B^11*(B-12)/(B-1)^12)
        SELECT 1  
        SKIP
        ENDFOR

?s
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
s=0.289234633707104

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

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

SELECT 1
   USE D:\最密三生素数文件\素数表.DBF ALIAS 素数表
   SELECT 2
   USE D:\最密三生素数文件\最密4生素数.DBF ALIAS 密4中
   kssj=SECONDS()
   SELECT 1
   GO 1
   For i=1 to 50876228
    @ 5,12 say i
   SELECT 1 &&打开盛放素数的表
   jl=recno()
   PUBLIC A
   A=素数
   SKIP 1
   PUBLIC B
   B=素数
   SKIP 1
   PUBLIC C
   C=素数
   SKIP 1
   PUBLIC D
   D=素数
   cj1=B-A &&求前后两个素数的差值
   cj2=C-B &&求前后两个素数的差值
   cj3=D-C &&求前后两个素数的差值
     IF cj1=2 AND cj2=4 AND cj3=2
     SELECT 2
     APPEND BLANK  &&增加一条空记录
     REPLACE szx WITH A+4 &&将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,"运行时间提示")

SELECT 1
   USE d:\最密三生素数文件\最密4生素数.DBF ALIAS 四中项
   SELECT 2
   USE d:\最密三生素数文件\统计最密4中项.DBF ALIAS 四中统计
   
   kssj=SECONDS()
   FOR i=1 TO 10
      @ 3,6 say i
        fw=10^i
        SELECT 1
        COUNT ALL FOR szx<fw TO tj
         SELECT 2
         APPEND BLANK     &&增加一条空记录
         REPLACE n WITH i
         REPLACE 范围 WITH fw
         REPLACE 总量 WITH tj
   ENDFOR
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")

SELECT 1
   USE D:\最密三生素数文件\素数表.DBF ALIAS 素数表
   SELECT 2
   USE D:\最密三生素数文件\最密5生素数.DBF ALIAS 密5中
   kssj=SECONDS()
   SELECT 1
   GO 1
   For i=1 to 50876228
    @ 5,12 say i
   SELECT 1 &&打开盛放素数的表
   jl=recno()
   PUBLIC A
   A=素数
   SKIP 1
   PUBLIC B
   B=素数
   SKIP 1
   PUBLIC C
   C=素数
   SKIP 1
   PUBLIC D
   D=素数
   SKIP 1
   PUBLIC E
   E=素数
   cj1=B-A &&求前后两个素数的差值
   cj2=C-B &&求前后两个素数的差值
   cj3=D-C &&求前后两个素数的差值
   cj4=E-D &&求前后两个素数的差值
     IF cj1=2 AND cj2=4 AND cj3=2 AND cj4=4
     SELECT 2
     APPEND BLANK  &&增加一条空记录
     REPLACE wss WITH A+6 &&将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,"运行时间提示")
用素数表制作最密5生素数中项

SELECT 1
   USE d:\最密三生素数文件\最密5生素数.DBF ALIAS 五中项
   SELECT 2
   USE d:\最密三生素数文件\统计最密5中项.DBF ALIAS 五中统计
   
   kssj=SECONDS()
   FOR i=1 TO 10
      @ 3,6 say i
        fw=10^i
        SELECT 1
        COUNT ALL FOR wss<fw TO tj
         SELECT 2
         APPEND BLANK     &&增加一条空记录
         REPLACE n WITH i
         REPLACE 范围 WITH fw
         REPLACE 总量 WITH tj
   ENDFOR
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
统计最密5生素数

SELECT 1
   USE D:\最密三生素数文件\素数表.DBF ALIAS 素数表
   SELECT 2
   USE D:\最密三生素数文件\最密5生素数逆.DBF ALIAS 密5中逆
   kssj=SECONDS()
   SELECT 1
   GO 1
   For i=1 to 50876228
    @ 5,12 say i
   SELECT 1 &&打开盛放素数的表
   jl=recno()
   PUBLIC A
   A=素数
   SKIP 1
   PUBLIC B
   B=素数
   SKIP 1
   PUBLIC C
   C=素数
   SKIP 1
   PUBLIC D
   D=素数
   SKIP 1
   PUBLIC E
   E=素数
   cj1=B-A &&求前后两个素数的差值
   cj2=C-B &&求前后两个素数的差值
   cj3=D-C &&求前后两个素数的差值
   cj4=E-D &&求前后两个素数的差值
     IF cj1=4 AND cj2=2 AND cj3=4 AND cj4=2
     SELECT 2
     APPEND BLANK  &&增加一条空记录
     REPLACE nwss WITH A+6 &&将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,"运行时间提示")
用素数表制作最密5生素数逆(0,4,2,4,2)

SELECT 1
   USE D:\最密三生素数文件\素数表.DBF ALIAS 素数表
   SELECT 2
   USE D:\最密三生素数文件\最密6生素数.DBF ALIAS 密6中
   kssj=SECONDS()
   SELECT 1
   GO 1
   For i=1 to 50876228
    @ 5,12 say i
   SELECT 1 &&打开盛放素数的表
   jl=recno()
   PUBLIC A
   A=素数
   SKIP 1
   PUBLIC B
   B=素数
   SKIP 1
   PUBLIC C
   C=素数
   SKIP 1
   PUBLIC D
   D=素数
   SKIP 1
   PUBLIC E
   E=素数
   SKIP 1
   PUBLIC F
   F=素数
   cj1=B-A &&求前后两个素数的差值
   cj2=C-B &&求前后两个素数的差值
   cj3=D-C &&求前后两个素数的差值
   cj4=E-D &&求前后两个素数的差值
   cj5=F-E &&求前后两个素数的差值
     IF cj1=4 AND cj2=2 AND cj3=4 AND cj4=2 AND cj5=4
     SELECT 2
     APPEND BLANK  &&增加一条空记录
     REPLACE lss WITH A+8 &&将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,"运行时间提示")