SELECT 等差3中升序
BROWSE LAST
COUNT ALL FOR dcsz<5000000 TO tj
?tj
GO 9888
USE
USE 等差三中和统计表 AGAIN IN 0
SELECT 等差三中和统计表
BROWSE LAST
COUNT ALL FOR dcszh<5000000 TO tj
GO 1499900
USE
USE 等差三中和统计表 AGAIN IN 0
SELECT 等差三中和统计表
BROWSE LAST
SELECT * FROM 等差三中和统计表 WHERE tj2=0
COUNT ALL FOR dcszh<7000000 TO tj
?tj
GO 13119
COUNT ALL FOR dcszh 5000000 TO tj
COUNT ALL FOR dcszh<5000000 TO tj
?tj
GO 11859
GO 1
COPY TO d:\三生素数中项\等差3生素数中项五百万前无合成数.xls NEXT 11859 TYPE XL5
USE
USE 等差3中升序 AGAIN IN 0
SELECT 等差3中升序
BROWSE LAST
COUNT ALL FOR dcsz<10000000 TO tj
?tj
GO 17194
USE
USE 等差三中和统计表 AGAIN IN 0
SELECT 等差三中和统计表
BROWSE LAST
GO bottom
USE

SELECT 1
 USE d:\三生素数中项\等差3中升序.DBF ALIAS 等差3中表
 SELECT 2
 USE d:\三生素数中项\等差三中和统计表.DBF ALIAS 等差三中和统计
 kssj=SECONDS()
 SELECT 1
 GO 9889 &&定位在第一个参与运算值上,本次运行交叉部分，外循坏取前段，内循环取后段，一次产生两组
 FOR i=1 TO 7306
   @ 3,6 say i
   jl=recno()
   PUBLIC A
   A=dcsz
      SELECT 1
      GO 9889
      FOR j=1  TO 7306
      PUBLIC B
      B=dcsz
      n=A+B
      c=n-INT((n-1)/30)*30
      d=INT((n-1)/30)*9
      IF c<21
      SELECT  2
      GO   INT(c/3)+d
      REPLACE tj2 WITH tj2+1
      ELSE
      SELECT  2
      GO 3+INT((c-10)/3)+d
      REPLACE tj2 WITH tj2+1
      ENDIF
      SELECT 1
      skip
      ENDFOR
      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:\三生素数中项\等差3中升序.DBF ALIAS 等差3中表
 SELECT 2
 USE d:\三生素数中项\等差三中和统计表.DBF ALIAS 等差三中和统计
 kssj=SECONDS()
 SELECT 1
 GO 1 &&定位在第一个参与运算值上,本次运行交叉部分，外循坏取前段，内循环取后段，一次产生两组
 FOR i=1 TO 17194
   @ 9,6 say i
   jl=recno()
   PUBLIC A
   A=dcsz
      SELECT 1
      GO 17195
      FOR j=1  TO 3902
      PUBLIC B
      B=dcsz
      n=A+B
      c=n-INT((n-1)/30)*30
      d=INT((n-1)/30)*9
      IF c<21
      SELECT  2
      GO   INT(c/3)+d
      REPLACE tj2 WITH tj2+2 &&本次运行（a+b）的平方交叉部分
      ELSE
      SELECT  2
      GO 3+INT((c-10)/3)+d
      REPLACE tj2 WITH tj2+2 &&本次运行（a+b）的平方交叉部分
      ENDIF
      SELECT 1
      skip
      ENDFOR
      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:\三生素数中项\等差3中元素.DBF ALIAS 等差元素
SELECT 2
USE d:\三生素数中项\等差三中和统计表.DBF ALIAS 等差三中和统计
kssj=SECONDS()
FOR i=666670 TO 1000003
   @ 3,6 say i
      SELECT 1
      GO 1
      FOR j=1  TO  9
        SELECT 1
        PUBLIC a
         a=dcys
         m=30*(i-1)+a
         SELECT  2
         APPEND BLANK   
         REPLACE dcszh WITH m
         REPLACE tj2 WITH 0
         SELECT 1
         skip
     ENDFOR
ENDFOR
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")

SELECT 1
 USE d:\三生素数中项\等差3中升序.DBF ALIAS 等差3中表
 SELECT 2
 USE d:\三生素数中项\等差三中和统计表.DBF ALIAS 等差三中和统计
 kssj=SECONDS()
 SELECT 1
 GO 1 &&定位在第一个参与运算值上,本次运行交叉部分，外循坏取前段，内循环取后段，一次产生两组
 FOR i=1 TO 17194
   @ 9,6 say i
   jl=recno()
   PUBLIC A
   A=dcsz
      SELECT 1
      GO 21097
      FOR j=1  TO 2526
      PUBLIC B
      B=dcsz
      n=A+B
      c=n-INT((n-1)/30)*30
      d=INT((n-1)/30)*9
      IF c<21
      SELECT  2
      GO   INT(c/3)+d
      REPLACE tj2 WITH tj2+2 &&本次运行（a+b）的平方交叉部分
      ELSE
      SELECT  2
      GO 3+INT((c-10)/3)+d
      REPLACE tj2 WITH tj2+2 &&本次运行（a+b）的平方交叉部分
      ENDIF
      SELECT 1
      skip
      ENDFOR
      SELECT 1
      GO jl+1
 ENDFOR
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
计算a∈（0,1千万），c∈（1千3百万，1千5百万），之前已经计算了a^2,2ab,b∈(1千万，1千3百万）

SELECT 1
USE D:\素数式系数\四生成元至23.DBF  ALIAS 四生元参照表
SELECT 2
USE D:\素数式系数\素数13万.DBF ALIAS 素数表十三万
SELECT 3
USE D:\素数式系数\最密四生素数.DBF ALIAS 四生素数
kssj=SECONDS()                      &&取出开始时间
    FOR i=1 TO 1
     @22,20 SAY i
     SELECT 1
     GO 1
     FOR k=1 TO 700245
     ly23=四中23+4
     bpz=ly23+(i-1)*223092870 &&给被判断值置数，以循环值做变量。
     Kf=INT(SQRT(bpz))                   &&求出被判断值的开方根
        SELECT 2                     &&打开素数表
        LOCATE FOR 素数>=kf                 &&根据开方根,查找最大素数
        DO CASE
        CASE EOF()                          &&如果超出素数表最后一条记录
        GO BOTTOM
        CASE 素数>kf                        &&如果找到的素数比开方根值大,就将记录指针退回一个,保证使用的最大素数在开方根内
        SKIP -1
        ENDCASE
        SSS=RECNO()                         &&记录下最大素数（开方根内）记录位置（如果找到的素数等于开方根，则直接读取这个素数的位置）
        GO 10
                                 &&从第二条记录开始读取素数(3)
        FOR j=1 TO SSS-9
        ys=MOD(bpz,素数)                    &&内循环开始。这个循环实质上是从小到大顺序,依次读取素数。循环值是记录序号
        IF ys=0 OR ys=2 OR ys=6 OR ys=8 && OR ys=12 OR ys=16 OR ys=MOD(36,素数) OR ys=MOD(38,素数) OR ys=MOD(120,素数) OR ys=MOD(122,素数) OR ys=MOD(126,素数) OR ys=MOD(128,素数)
        EXIT                                &&因为符合条件,则做完相应工作后跳出内循环
        ENDIF
        SKIP                               &&素数表指针向下移动一个
        ENDFOR
        IF j>sss-9
        SELECT 3              &&打开保存求解结果的信息表
        APPEND BLANK                        &&增加一条空记录
        REPLACE 四素 WITH bpz             &&将被判断值保存到被判断值字段中(被判断值数据栏)
        ENDIF
        SELECT 1
        SKIP
        ENDFOR
ENDFOR
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
用生成元及素数表制作最密四生素数