|

楼主 |
发表于 2023-6-28 12:59
|
显示全部楼层
CLEAR ALL
SELECT 1
USE d:\孪生素数\素数表万.DBF ALIAS 万素
SELECT 2
USE d:\孪生素数\孪素19.DBF ALIAS 孪素19基表
SELECT 3
USE d:\孪生素数\孪生素数表.DBF ALIAS 孪素表
kssj=SECONDS() &&取出开始时间
For i=1 to 10
@ 3,6 say i
SELECT 孪素19基表
go 1
SELECT 孪素19基表
for m=1 to 378675
k=孪19
N=k+9699690*(i-1)
Kf=INT(SQRT(N)) &&求出奇数的开方根
SELECT 万素 &&打开素数表
LOCATE FOR 素数>=kf &&根据开方根,查找最大素数
DO CASE
CASE EOF() &&如果超出素数表最后一条记录
GO BOTTOM
CASE 素数>kf &&如果找到的素数比开方根值大,就将记录指针退回一个,保证使用的最大素数在开方根内
SKIP -1
ENDCASE
SSS=RECNO() &&记录下最大素数(开方根内)记录位置(如果找到的素数等于开方根,则直接读取这个素数的位置)
GO 9 &&从第一条记录开始读取素数(3)
FOR j=1 TO SSS-8 &&内循环开始。这个循环实质上是从小到大顺序,依次读取素数。循环值是记录序号
qmz=MOD(N,素数) &&以读取的素数为条件,对被判断值求模 【qmz=MOD(bpz,素数)】文本号括着的原在此行语句体部分.特别需要注意的是在用标点符号时,一定改成英文模式及半角状态。
IF qmz=0 OR qmz=2
EXIT
ENDIF
SKIP &&素数表指针向下移动一个
ENDFOR
IF j>sss-8
SELECT 孪素表 &&打开盛放素数的表
APPEND BLANK &&增加一条空记录
REPLACE 孪素 WITH N &&将N值付给素数
ENDIF
SELECT 2
skip
ENDFOR
ENDFOR
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
这是筛选10个周期中的孪生素数程序(周期值9699690). |
|