***
*** 统计X以内的超级素数的个数 By Tiger5392 (C)Copyright 2006-2006
***
*** 语法
*** ?超级素数(nExpr)
*** nExpr:正整数,表示要计算从1到该数之间的超级素数
*** 返回:超级素数的个数
***
*** 相关知识
*** 超级素数:一个素数,去掉个位是素数;去掉个位和十位也是素数;
*** 直到保留到最高位依然是素数,这个素数称为超级素数。
***
PARAMETERS nNum
PRIVATE nNum,tStartTime,nTotal,nCnt,lIsSupperPrime
CLEAR
tStartTime=SECONDS()
nTotal=0 &&超级素数的个数的初始值
FOR nCnt=1 TO nNum &&nNum以内的数都需要进行判断
IF IsSupperPrime(nCnt)
nTotal=nTotal+1
* ?STR(nTotal,5)+STR(nCnt,10)+STR(SECONDS()-tStartTime,10,3)
ENDIF
ENDFOR
?"Total Time:"+STR(SECONDS()-tStartTime,10,3)+" Seconds"
?"Total Prime Numbers:"+STR(nTotal,5)+"(1-"+ALLTRIM(STR(nNum))+")"
RETURN nTotal
*** EOF(): CSP By Tiger5392 (C)Copyright 2006-2006
FUNCTION IsSupperPrime
PARAMETERS nNum
PRIVATE nNum,nNum1,cNum,nCnt,lIsSupperPrime
IF SUBSTR(ALLTRIM(STR(nNum)),1,1)$"14689"
RETURN .F. &&若数值最高位不为2,3,5,7之一,则判定该数不为超级素数
ENDIF
FOR nCnt=0 TO LEN(ALLTRIM(STR(nNum)))-2
nNum1=INT(nNum/10^nCnt) &&将原数去掉从个位开始的nCnt位
cNum=SUBSTR(ALLTRIM(STR(nNum1)),2) &&提取除高位外的其余部分
IF "0"$cNum OR "2"$cNum OR "4"$cNum OR "5"$cNum OR "6"$cNum OR "8"$cNum
RETURN .F. &&要判断的数除高位外其余数值中含有0,2,4,5,6,8之一,则判定不为素数
ENDIF
IF IsPrime(nNum1)=.F.
RETURN .F. &&去掉nCnt位后不为素数,则判定不为超级素数
ENDIF
ENDFOR
RETURN .T.
*** EOF(): IsSupperPrime By Tiger5392 (C)Copyright 2006-2006
FUNCTION IsPrime
PARAMETERS nNumber
PRIVATE nNumber,nCnt
FOR nCnt=2 TO SQRT(nNumber)
IF nNumber%nCnt=0
RETURN .F.
ENDIF
ENDFOR
RETURN .T.
*** EOF(): IsPrime By Tiger5392 (C)Copyright 2006-2006