|

楼主 |
发表于 2010-9-13 10:55
|
显示全部楼层
[原创]请教vfbpgyfk先生一个编程问题
[这个贴子最后由白新岭在 2010/09/13 11:00am 第 1 次编辑]
以下是182楼修改后的程序,用来求对称8生素数群,首先需要制作表:k8生成初始值29,里边存放着可能产生8生素数群的代表元素,共计2*4*6*(17-8)*(19-8)*(23-8)*(29-8)=1496880个元素,分段周期为:2*3*5*7*11*13*17*19*23*29=6469693230.如果筛选10^12内的对称8生素数群,需要循环154.5668个周期,即不到155个分段周期。总处理数据量=1496880*154.5668=231367951,即2.3亿个数据。
CLEAR ALL
SELECT 1
USE E:\用vfp求k生素数数量\最密8生素数群.DBF ALIAS 最密8生素数
SELECT 2
USE E:\用vfp求k生素数数量\素数表2.DBF ALIAS 素数表
SELECT 3
USE E:\用vfp求k生素数数量\k8生成初始值29.DBF ALIAS 生成初始值29
*INPUT "请输入预先值 K= " TO yxk
*INPUT "请输入步长值 bcz= " TO bcz
*INPUT "请输入初始值 csz= " TO csz
*INPUT "请输入外循环起始值 xks= " TO xks
*INPUT "请输入外循环终结值 zds= " TO zds
GO top &&***************这条语句是将“k8生成初始值31”表的指针移到“41001”处**************
DO WHILE NOT EOF() &&*********既然是调用表中信息,而且是后部分,就应该使用“DO WHILE NOT EOF()”语句,由它来判断是否到文件尾**********
@ 2,5 say RECNO() &&显示记录号
scy=生成元29 &&计算被判断值
For i=1 to 2
bpz=scy+i*6469693230
Kf=INT(SQRT(bpz)) &&求出被判断值的开方根
SELECT 素数表 &&打开素数表
LOCATE FOR 素数>=kf &&根据开方根,查找最大素数
DO CASE
CASE EOF() &&如果超出素数表最后一条记录
GO BOTTOM
CASE 素数>kf &&如果找到的素数比开方根值大,就将记录指针退回一个,保证使用的最大素数在开方根内
SKIP -1
ENDCASE
SSS=RECNO() &&记录下最大素数(开方根内)记录位置(如果找到的素数等于开方根,则直接读取这个素数的位置)
GO 10 &&从第二条记录开始读取素数(37)
FOR j=10 TO SSS &&内循环开始。这个循环实质上是从小到大顺序,依次读取素数。循环值是记录序号
qmz=MOD(bpz,素数) &&以读取的素数为条件,对被判断值求模 【qmz=MOD(bpz,素数)】文本号括着的原在此行语句体部分.特别需要注意的是在用标点符号时,一定改成英文模式及半角状态。
IF qmz=0 OR qmz=2 OR qmz=6 OR qmz=12 OR qmz=14 OR qmz=20 OR qmz=24 OR qmz=26
EXIT
ENDIF
SKIP &&素数表指针向下移动一个
ENDFOR
IF j>sss
SELECT 最密8生素数 &&打开保存求解结果的信息表
APPEND BLANK &&增加一条空记录
REPLACE 被判值 WITH bpz &&将被判断值保存到被判断值字段中(被判断值数据栏)
REPLACE 开方根 WITH kf &&将开方根保存到开方根字段中(开方根值数据栏)
REPLACE 外环值 WITH RECNO()
ENDIF
ENDFOR
SELECT 生成初始值29
SKIP
ENDDO
CLEAR ALL
文字 |
|