SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
k=4
SELECT 1
   GO 3
    s=1.0000000000000000
       FOR  i=1 TO 50876320  &&又遇到以前类似事情，外循环变量用了j与下面的内循环变量一致，所以提示0不能做除数。
        @ 5,12 say i
        SELECT 1
        p=素数
        jl=recno()
        s=s*p^(k-1)*(p-k)/(p-1)^k
         SELECT 1
            GO jl+1
     ENDFOR
?s
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
好久不发编程了，这是计算孪中合成数的系数极限值。

SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
k=4
SELECT 1
   GO 3
    s=1.0000000000000000
       FOR  i=1 TO 50876320  &&又遇到以前类似事情，外循环变量用了j与下面的内循环变量一致，所以提示0不能做除数。
        @ 5,12 say i
        SELECT 1
        p=素数
        jl=recno()
        s=s*p*(p-k)/(p-2)^2
         SELECT 1
            GO jl+1
     ENDFOR
?s
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
孪中合成数的数量公式中的系数（一个式子的极限值）:0.3968803638672010
运行时间:  32分46.79秒   确定。
那个值*6就是C2(2)了
6*∏P(P-4)/(P-2)^2=6*0.3968803638672010=2.3812821832032100

SELECT 1
USE D:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
s=1.0000000000000000
SELECT 1
  GO 5
  A=素数
  s=s*A^7*(A-8)/(A-1)^8
  SKIP 1
        FOR  j=1 TO 50876310
        SELECT 1
        B=素数
        s=s*A^7*(A-8)/(A-1)^8
        SELECT 1  
        SKIP
        ENDFOR

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

因为VFP版块的排序限制，即便是新发的帖子，它也不会出现在第一页。

SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
s=3.000000000000000000
SELECT 1
  GO 3
  A=素数
  s=s*A*(A-3)/(A-1)/(A-2)
  SKIP 1
        FOR  j=1 TO 50876310
        SELECT 1
        B=素数
        s=s*B*(B-3)/(B-1)/(B-2)
        SELECT 1  
        SKIP
        ENDFOR

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

SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
s=2.000000000000000000
SELECT 1
  GO 2
  A=素数
  s=s*A*(A-2)/(A-1)^2
  SKIP 1
        FOR  j=1 TO 50876310
        SELECT 1
        B=素数
        s=s*B*(B-2)/(B-1)^2
        SELECT 1  
        SKIP
        ENDFOR

?s
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
这是求二生素数的公共系数。

SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
s=5.000000000000000000
SELECT 1
  GO 4
  A=素数
  s=s*A*(A-5)/(A-2)/(A-3)
  SKIP 1
        FOR  j=1 TO 50876310
        SELECT 1
        B=素数
        s=s*B*(A-5)/(A-2)/(A-3)
        SELECT 1  
        SKIP
        ENDFOR

?s
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
计算孪中+最密3生素数的中项合成数的数量公式前的系数（公共系数）

SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
s=4.500000000000000000
SELECT 1
  GO 3
  A=素数
  s=s*A^2*(A-3)/(A-1)^3
  SKIP 1
        FOR  j=1 TO 50876310
        SELECT 1
        B=素数
        s=s*B^2*(B-3)/(B-1)^3
        SELECT 1  
        SKIP
        ENDFOR

?s
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
最密3生素数（P，P+2，P+6）的系数

SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
s=8.7500000000000000000
SELECT 1
  GO 5
  A=素数
  s=s*A*(A-7)/(A-3)/(A-4)
  SKIP 1
        FOR  j=1 TO 50876310
        SELECT 1
        B=素数
        s=s*B*(B-7)/(B-3)/(B-4)
        SELECT 1  
        SKIP
        ENDFOR

?s
=MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
最密3生素数的中项与最密4生素数的中项和，合成数的数量公式中的公共系数。

MODIFY PROJECT e:\三角递增法\三角项目.pjx
USE 四生素数表中 AGAIN IN 0
SELECT 四生素数表中
BROWSE LAST
?RECCOUNT()
GO bottom
GO 1000
GO 200
COUNT ALL FOR 四中<1000000 TO tj
?tj
GO 1
COPY TO e:\三角递增法\四生素数中项.xls NEXT 165 TYPE XL5
USE
MODIFY PROJECT f:\最密三生素数\正三生素数.pjx
USE 三生素数d24 AGAIN IN 0
SELECT 三生素数d24
BROWSE LAST
GO bottom
COUNT ALL FOR 三中24<1000000 TO tj
?tj
GO 1392
GO 1
COPY TO f:\最密三生素数\三中表.xls NEXT 1392 TYPE XL5
USE

今天的实际操作。