超级素数问题(4)

作者在 2006-10-10 20:06:00 发布以下内容

***
***  统计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

VFP程序 | 阅读 3791 次
文章评论,共0条
游客请输入验证码
浏览119096次