[原创]请教vfbpgyfk先生一个编程问题

白新岭

截止2022年11月02日周三20:34分浏览量78164，回复1763，热度111度
截止2022年11月03日周四21:34分浏览量78274，回复1764，热度111度

白新岭

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就是\(C_{2(2)}\)了
6*∏P(P-4)/(P-2)^2=6*0.3968803638672010=2.3812821832032100

独木星空谁

截止2022年11月03日周四21:34分浏览量78274，回复1764，热度111度
截止2022年11月04日周五06:31分浏览量78319，回复1767，热度112度

独舟星海

截止2022年11月03日周四21:34分浏览量78274，回复1764，热度111度
截止2022年11月04日周五14:36分浏览量78350，回复1768，热度112度

独木星空谁

截止2022年11月03日周四21:34分浏览量78274，回复1764，热度111度
截止2022年11月04日周五23:06分浏览量78382，回复1769，热度112度

独舟星海

截止2022年11月04日周五23:06分浏览量78382，回复1769，热度112度
截止2022年11月05日周六12:39分浏览量78452，回复1770，热度112度

独木星空谁

截止2022年11月04日周五23:06分浏览量78382，回复1769，热度112度
截止2022年11月05日周六20:06分浏览量78482，回复1771，热度112度

独舟星海

截止2022年11月05日周六20:06分浏览量78482，回复1771，热度112度
截止2022年11月06日周日11:11分浏览量78532，回复1772，热度112度

独木星空谁

截止2022年11月05日周六20:06分浏览量78482，回复1771，热度112度
截止2022年11月06日周日23:27分浏览量78592，回复1773，热度112度