|
|
SELECT 1
USE D:\素数式系数\素数表万.dbf ALIAS 素数表
SELECT 2
USE D:\素数式系数\k生素数系数表26的.dbf ALIAS 系数表
kssj=SECONDS() &&取出开始时间
FOR n=10 TO 1000
@ 5,12 say n
s=1.0000000000000000
IF MOD(6*n,13)=0
s=s*(13-1)/(13-2)
ENDIF
SELECT 1
GO 3
zhi=.f.
DO WHILE NOT zhi
SELECT 1
P=素数
IF 6*n>=P
zhi=.f.
ELSE
zhi=.t.
ENDIF
DO case
CASE MOD(6*n,P)=0
s=s*(P-2)/(P-4)
CASE MOD(6*n,P)=MOD(26,P)
s=s*(P-3)/(P-4)
CASE MOD(6*n,P)=MOD(-26,P)
s=s*(P-3)/(P-4)
otherwise
s=s*1
ENDCASE
SELECT 1
SKIP
ENDDO
SELECT 2
APPEND BLANK &&增加一条空记录
REPLACE k生 WITH 4 &&将B值付给素数式
REPLACE 二kd WITH 6*n &&将B值付给素数式
REPLACE 系数 WITH 5.07366549990956*s &&将B值付给素数式
ENDFOR
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示") |
|